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FIGURES D'ÉQUILIBRE 

D'UNE MASSE FLUIDE 



CHAPITRE PREMIER 



THÉORÈMES GÉNÉRAUX SUR LE POTENTIEL NEWTONIEN 



Objet du cours, — Nous allons étudier la figure d'équilibre 
d'une niasse Uuide animée d'un mouvement de rotation : les seules 
forces auxquelles est soumis le système sont des forces intérieures 
dues h l'attraction newtonienne : deux points s'attirent en raison 
directe de leurs masses et en raison inverse du carré de leur 
distance. Nous allons d'abord rappeler quelques résultats connus 
sur le potentiel newtonien. 

Définition et propriétés du potentiel. — Un point de masse 
m, soumis à l'attraction de tout le reste du système, est soumis 
à une force dérivant d'un potentiel V, c'est-à-dire que les com- 
posantes, suivant les axes de coordonnées, de la force agissant sur 
ce point, sont les dérivées partielles d'une môme fonction 

dV dV OV 

m -r — , m -r — , m 



— — > //* -r — j lit. -^ — • 
^x ùj/ dz 

Si les masses attirantes sont discontinues, le potentiel a pour 
expression 

en appelant ;*, la distance du point de masse m au point de 
masse /w^. 

PoiNCAUÉ. — Figures d'équilibre. i 
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THEOREMES GENERAUX SUR LE POTENTIEL NEWTON lEN 

Si ces masses sont continues, le potentiel a pour expression 




p' étant la densité de l'élément de volume d'z\ r étant sa distance * 

au point /;^, et l'intégrale étant étendue à tout le volume attirant. 
On sait que l'intégrale a un sens. 

Si la masse attirante est répandue sur une surface, avec une 
densité superficielle \k\ le potentiel a pour expression 



ix'(h' 




Enfin, si c'est un volume infiniment mince d'épaisseur s', alors 
et le potentiel est 






Pour un volume, en dehors des masses attirantes, V est une 
fonction continue, ainsi que ses dérivées premières, dans tout 
l'espace, et on sait que l'on a 

iVV cVV D-V 



AV = -^-^ +-^-v- +- Tir = o- 



A rintérieur, en un point de densité o, on a 

AY = — 4'^?. 

Les dérivées secondes sont discontinues en une surface de sépa- 
ration de deux milieux diflerents. 

Pour une surface, V est continue dans tout l'espace, mais ses 
dérivées ne le sont pas sur la surfiice elle-même. À 

Sur la normale en un point A d'une surface S (fig. i), prenons 

deux longueurs 

AB = dn,.y AC = d/ii. 
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Si V^ est le potentiel en A, le potentiel en B et C est 

On a craillcurs v\g, , 

dV i>V 

+-xr-=— 4-;^' 




0/2,. 0/// 

jjL étant la densité de la masse attirante. 

Si Ton désigne par /, m, n les cosinus directeurs de la normale, 
supposée dirigée vers l'extérieur, par 

les valeurs des dérivées partielles à l'intérieur et à l'extérieur de 
la surface, on a 

- S-==^ (-17).+ '" (-17).+ " (if).- 

Théorème de Gauss, — Gauss a démontré que, en dehors 
des masses attirantes, V n'a ni maximum ni minimum, puisque 
l'on a 4V= o. 

D'une façon générale, si AV est positif ou nul, la fonction V 
n'a pas de maximum; si AV est négatif ou nul, il n'y a pas de 
minimum. Dans le cas qui nous occupe, AV est nul ou négatif: V 
n'a pas de minimum et n'a pas de maximum en dehors des 
masses agissantes. 

Donc, si sur une surface fermée ne renfermant pas de masses 
attirantes, V est compris entre deux nombres ^ et /f, il en 
résulte qu'en tout point situé à l'intérieur de la surface, on aura 
ff<Y<h. 
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4 TÛEOREMES GENERAUX SUR LE POTENTIEL NEWTOMEN 

Théorème de Green, — Sî U et V sont deux fonctions con- 
tinues h rintéricur d'un volume T limité par une surface S, on a 
l'égalité 



(.) / U41rf^= i UAVrf-c 




J. 



I / 5U 5V dU tW OU dV \ 
"^ / \ dj: dx "^ dy dy "*" d; d^ j ' ' 



Dans celte égalité, on peut permuter U et V, et en retranchant 
les deux égalités obtenues, on a 



(^) / (u^-^1^)^' - / (U^V-VAU)rf.. 





Cas particuliers de l'égalité (i) : faisons U == i. 



/^ 



Faisons : 




(3) I ^ rf„ I AVrfT. 



u=v 



(4) / V -1^ rf^= / VAVrfT 




D/i 




ÔV\2 /dV' 






)V(-f){V(-f)']'''. 



Si V est une fonction satisfaisant à l'équation de Laplace, 
AV=o, 



(=) /^^*-/[(-S)V(f )V(4^)>. 



le second membre étant positif, le premier Test aussi. 

Si U et V sont des potentiels, sur une sphère de rayon suffi- 
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samment grand, U et V sont de l'ordre de — , leurs dérivées 

premières de Tordre de -^ . On peut alors appliquer la for- 
mule (i) à l'espace extérieur à une surface S et intérieur à une 
sphère de rayon très grand. Si ce rayon est suffisamment grand, 
l'intégrale prise sur la surface de la sphère cp est négligeable : 

c'est-h-dire est de Tordre de -rr- ; mais on a la relation 



''^/^+ ; u-^'''^" 



car les éléments de Tintégrale se détruisent deux à deux. 

Quant aux intégrales des seconds membres, leur somme donne 
des intégrales étendues h l'espace entier, et on a, en ajoutant les 
deux égalités obtenues : 

les intégrales étant étendues à tout l'espace. 
De même, la formule (2) donnera 



(UAV— VAU)rfT = o; 

et, en se rappelant que 

— AV = 4^? 
_AV, = 47rp„ 

on a 



(7) / (V.?-VpJrfT=o 




l'intégrale étant étendue à tout l'espace. 
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Cette formule est analogue à celle que Ton trouve dans la 
théorie de rélectrlcité, 



2!/^"^ '~2j"^'^ "^ ^' 



Dans la formule (^), faisons 
on en conclut la formule 




p — c,d\)d'z = o. 



Travail de Fattraction. — Supposons un système de masses 
oi\ m" ^ ..., etc. attirées par un autre /«',, w"^, ...., etc. 

Soit Vj, V',,... le potentiel aux points m\ m'\...j e\.c., dû aux 
masses attirantes w'^, in'\^ ..., etc. 

V, V"..., le potentiel aux points m\^ m'\,,,^ dû aux points 



m' ^ m" , 



Les masses attirées se déplaçant, on a comme travail 
en appelant Sx, Sy, S::, le déplacement du point de masse m ; 



Si je pose 



j aurai 






ê=:ôH. 



Si les masses attirantes se déplacent aussi, mais seules, le 
potentiel V, deviendra 



\ 



THEOREME DE GREEN 

et on a 

5^ 



;ii=V/«o'\^. 



Si les deux déplacements ont lieu simultanément, on a 
Supposons que les masses attirées forment un volume T, on a 




11= I ?V,rfT, 



cette intégrale pouvant d'ailleurs s'étendre à Tespace entier, 
puisque p est nul en dehors de T. On a 



oII= I V.oorfT 




étendue de même à l'espace entier. 

Si la masse attirante est la même que la masse attirée, on 
a V = V., 

211= / Xd^dx, 
ce qui peut s'écrire en tenant compte de (8) 

c . 

ôH = I pSVcÎT, 



ou encore 



5„=. / ^^+JV d.. 



Posons maintenant 



VV= I ^di, 
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rintégrnie étant toujours étendue à tout Tespace; West Ténergie 
du f^ystcme. 
On a 

oW = t. 

On peut d'ailleurs écrire, puisque AV = — ^r.z. 



w = -^Jy^yd. 



ft, en tenant compte de (6), 



(9) 



-=w/[(^)V(f)V(|^)>. 



rintégralc étant toujours étendue à Tespace entier. 

Conditions d'équilibre. — Ceci étant rappelé, nous allons com- 
mencer Tétudc du problème que nous nous sommes proposé. 

Considérons une masse fluide, isolée de toute influence exté- 
rieure, tournant autour d'un axe fixe que nous prendrons comme 
axe Ozj d'un mouvement uniforme de vitesse angulaire to; prenons 
les axes Ox et Oy, invariablement liés à la masse fluido. 

Soity^ la pression en un point (a:, y, z); p ne dépend que de 
r, //, Zy la force s'exerçant sur l'élément de volume di a pour 
composantes 

X^^rf., Y=-^rfT, Z^^d-:. 
ô.r iSy UG 

Mcrivons les conditions d'équilibre relatif, en remarquant que 
raccéléraiion de Coriolis est nulle. Si on appelle X, Y, Z les 
composantes suivant les axes de coordonnées de la force agis- 
sant sur la molécule (.r, y, z), on a 

V "^^^ I s 



» * 



et, en posant 
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u=v+-^-(x^ + r). 



on a 



c> dU 

— 



0.r ' 0^- ' 



d'où Ton tire 



-:r — dx-^—^clu 4--J—dz==p (- — rfx+ -r — dt/ -+~ —— dz] , 

dp=Qd\J. 

Cette relation montre que p est uniquement fonction de p ; U csl 
par suite, aussi uniquement fonction de p et, par suite, de p. 
On arrive alors au théorème suivant: 

Théorème. — Quand une masse fluide, animée d'un mouve- 
ment de rotation autour d'un axe fixe^ est en équilibre relatifs les 
surfaces de niveau sont les surfaces d'égale pression^ et aussi 
d^ égale densité, 

A la surface de la masse, on a /;= o. Donc U est une cons- 
tante. La surface libre est donc une surface de niveau. 

S'il n'y a pas de rotation, on a U == V. 

Dans le cas général on a AV = — 4^p. AV est donc simple- 
ment fonction de U. 

o 1 r , , OU ÙU dU , , , 

Remarquons que la iorce totale -r — , rr — > -:r— est en gênerai 

normale aux surfaces de niveau, U = C. 

Ces conditions d'équilibre sont nécessaires. 

Une condition nécessaire et suffisante d'équilibre est que le 
travail résultant d'un déplacement virtuel soit nul. Ce travail 
comprend le travail de l'attraction, plus le travail dû à la force 
centrifuge. 
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Le premier a pour expression oW, que nous avons déjà 
étudié, le second a pour expression 



pc^T (co^.roa: + w"y2y) 




,/.oi^!if!±^ 






et, si J est le moment d'inertie par rapport à O-S, le travail dû à 
la force centrifuge est 



J. 



La condition d'équilibre est donc 



2W+i:^oJ=o 



pour tout déplacement compatible avec les liaisons. 
Cette condition n'implique pas nécessairement que 



AV + — J 

2 



soit maximum ; il faut encore d'autres conditions ; mais je dis 
que si cette fonction est maximum, l'équilibre est stable. 
Soit T la force vive relative, on a l'égalité 



_^w+^j)=e. 



O)* 



Posons W -|- — J =1], et supposons que pour le maximum 

U soit égale à U^. 
• On a 



T-U=T,-U„. 



' * 



RELATION ENTRE LA MASSE, LE VOLUME ET LA VITESSE n 

Si on donne maintenant à U une valeur U < U^ et à T^ une 
valeur a, on aura 

T-U = T,-U, 

T = a-(U.-U). 

La force vive dans les instants suivants sera inférieure à Ty, 
Téquilibre est donc stable. 

Relation entre la masse^ le volume et la vitesse de rotation. 
— Une condition nécessaire de l'équilibre est, qu'en tout point 
de la surface libre, la force soit dirigée vers Tintérieur, autrement 

une partie se détacherait. Il faut donc -- — <o et par suite 




La formule (3) transforme cette inégalité en la suivante : 




AUrfT<o. 



Donc il faut que Ton ait 




AVrfT+ 1 2(x}hk<0 



c'est-a-dire 






En désignant le volume par T et la masse par M, il vient 

— 471M + 2(0^1 <o. 

Dans le cas particulier où la densité est constante et égale à p, 
on a 

2tO^<4îîp. 
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12 THEOREMES GENERAUX SUR LE POTENTIEL NEWTONIEN 

Dans Tapplication de cette formule, il faut se rappeler que les 
unités sont telles que la force d'attraction newtonienne s'ex- 
prime par la formule 

^ mm' 



CHAPITRE II 



MASSE HOMOGENE FLUIDE 



LA MASSE EST SANS MOUVEMENT DE ROTATION 

Simplîûcation des formules générales. — Etudions le cas 
particulier d'une masse homogène fluide et sans mouvement de 
rotation. 

En appliquant les théorèmes généraux, on trouve que la sur- 
face libre doit être une surface équipotentielle V = V^. 

Si on prend la densité de la masse pour unité de densité, le 
potentiel en un point a pour expression 




L'énergie est 




Yd-z 



V est continu dans tout l'espace ainsi qu'à la surface libre ; il est 
nul à l'infini ; h l'extérieur il n'a ni maximum ni minimum et est 
compris entre o et V^j ; à l'intérieur il est supérieur à V^ et a un 
maximum. 



L'intégrale j r — ch étendue h une surface fermée entou- 



i4 MASSE HOMOGÈNE FLUIDE 

rant le corps a pour valeur i AV^t, étendue au volume limité 

par cette surface, et on a I AV^t = — 4"^ / prfT = — 4''^^ï" 

Ce théorème est encore vrai pour la surface limite, et il s'ap- 
plique que le corps soit creux ou plein, et même si Ton a affaire 
h une distribution superficielle. 

Comme enfin nous supposons le corps homogène et de den- 
sité 1, on a 



Wd'z = — I ^Tz^d'z == — 4'^T 



en désignant par T son volume, et par S une surface qui le 
renferme complètement, en particulier sa propre surface. 

Si on porte sur la normale en M h la surface extérieure un 
segment MM' = 6?/ic(fig. 2), le potentiel en M'a pour expression : 

Or, on a 

Fi g. '2. y^f <C Vj, ; 





on en conclut que 



<o. 



0/2 



e 



On a vu que cette condition était nécessaire pour l'équilibre : 
ici elle est toujours remplie. 

Soit Vq le potentiel de la masse fluide sur sa surface, et suppo- 
sons maintenant une couche électrique de masse M, en équilibre 
à la surface du corps, et produisant dans Tespace un potentiel 
V'tel que sur la surface et à son intérieur le potentiel soit V'=Vq. 

A l'extérieur, on a V =V. En effet, sur la surface V' = Vo et 
à l'infini V'= o. Or V et V sont deux fonctions satisfaisant h 
l'équation de Laplace, ayant les mêmes valeurs sur la surface et 
sur une sphère de rayon très grand. Ces deux fonctions coïnci- 
dent donc : c'est le principe de Dirichlet. 



THÉORIE DE LIAPOVNOFF i5 

On a, d'après des formules connues, 



ô/i 



On a d'ailleurs 




d'7= — 4^T. 



0/2 



Donc la masse M répandue h la surface a pour valeur T et réci- 
proquement si 

M = T, V' = V. 

T , 

Le rapport -^r^ est ce que Ton appelle la capacité électrique 

du corps. 

Cette notion va nous être utile pour résoudre le problème 
proposé. 

Théorie de Liapounoif, — Parmi les figures d'équilibre, il y 
a une solution évidente qXii est la sphère : mais peut-il y avoir 
d'autres figures d'équilibre ? Nous n'en savons rien : mais si Ton 
se propose de chercher les figures d'équilibre stable, le pro- 
blème est différent. 

Le maximum absolu de W est bien atteint pour la sphère, 
mais a priori nous ne pouvons affirmer que W n'ait pas d'autres 
maximums relatifs. Cela a été démontré par Liapounoff; la 
démonstration est assez longue et on va procéder par étapes. 

L'énergie a un maximum, — Le potentiel d'un volume 
quelconque par rapport à un point quelconque, est inférieur au 
potentiel d'une sphère de même volume par rapport à son centre. 

Du point M, comme centre, je décris la sphère S ayant le môme 
volume que le corps S (fig. 3). 

Le corps S est partagé en deux volumes : le volume K inté- 
rieur à S, et le volume A extérieur. 

Soit B le volume de S extérieur à S, 



J^ 
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On a par b^'potbêse B = A. 

Si R est le rayon de la sphère, on a 




Potentiel de A < 



A 



Potentiel de B 



U 
B 



A 



> 



Fig. 1. 



R R 

D'où 

Potentiel de B > Potentiel 
de A 

Potentiel de S = Potentiel 
de K -f- Potentiel de A 

Potentiel de i! = Potentiel 
de K 4- Potentiel de B. 



Donc 



Potentiel de I > Potentiel de S. 
Si V est le potentiel au point M, et R le rayon de la sphère I, 
on a donc 



V < 27:R-. 



L'énergie qui est égale à 




\ih 



'A 



est donc inférieure b 




et il vient 



W<-R^T 



en appelant T le volume de la masse. 
On a d'ailleurs 

T=-|-^R% 



d'après la définition de R. 
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Donc 

L'énergie W due à un volume T est donc limitée supérieure- 
ment : W a donc un maximum. Je dis que ce maximum n\^st 
atteint que si la surface extérieure est une sphère. 

Nous allons d*abord démontrer le théorème suivant : 

La capacité électrostatique a un minimum. — De tous les 
corps ayant même volume, la sphère est celui qui a la moindre 
capacité électrostatique. 

En effet, la formule ^9) (p. 8) donne 

rintégrale étant étendue h tout Tespace. 

Supposons que sur un corps conducteur de volume T il y ait 
une masse électrique T en équilibre. Si V est le potentiel corres- 
pondant, sur la surface du corps on a 

CV = T 

par définition de la capacité C du conducteur. 

A rintérieur du corps, V est constant; aTextérieur, V satisfait 
à l'équation AV = o et est nul h Tinfini. 

L'énergie électrique est, d'après une formule connue de Télec- 
trîcité : 

1 ! T^ 

\V = — CV'= — -7T-- 

2 2 L 

Elle est d'ailleurs égale à 



hfim'H^ÏHW-^ 



8 



l'intégrale étant étendue, soit a tout l'espace, soit au volume 
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la 
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extérieur seulement, puis({ue V est constant à rinlérîeur du 
A'olume. 

Nous allons voir quVlle est inférieure à Ténergie AV due an 
potentiel newtonien du corps supposé homogène et de den- 
sité égale à l'unité. 

Soit V== V+U ce potentiel : je pose 



\V=- 




i[^^^]-[^^]+[^^]!-. 



rintégralc étant étendue à tout l'espace. 




m'^mv- 



dv ou . OV ÙU 



*>!/ ^!/ 



Oc Oc 



d'Z 




La première intégrale est W. 

La deuxième se décompose en deux : l'une prise à l'intérieur 

OV OV 



du corps, qui est nulle, puisque 



OV 



= G. 



D.i- 0// ùr 

L'autre, prise à l'extérieur, se transforme par la formule de 



Green en 




\ — — an 



— i VAUrfT. 



La seconde de ces deux intégrales est nulle puisque AU est nul 
à l'extérieur. 

La première est également nulle : en elFet, sur la surface, V 
est constant et égal à Y',,. 
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ï9 



L'intéffrale est donc 




or 



ao" = 4'«^ 1 , 



d//. 



(h = 4t:T ; 



0/?, 



celle intégrale est donc nulle. 
On a par suite 




dV OU ùV i)U dV ùU . , 

«7=0. 



d.r 0.f dy 0// dz àz 



La troisième n'est pas nulle, car U n'est pas nul à l'intérieur 
du corps, puisque V y est constant, et que V y varie. 
La troisième intégrale est alors positive, et on a bien 

On a donc, en réunissant les deux résultats précédents : 

«R^T>W> ^' 



d'où 

T 
C> 



27:R-' 
et enfin 
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C a donc un minimum. Je dis que ce minimum ne peut être 
atteint que pour la sphère. 

Le minimum est atteint pour la sphère. — Supposons qu'un 
corps conducteur ait la forme pour laquelle ce minimum est 
atteint : on a 



'AV = 



^=ir/[(^)'-(^)'+(^)"J- 



"/ 




-r — désignant la dérivée suivant la normale a la surface équl- 
potentielle passant par le point (.r, y, ^), 



2W = 



Déformons le conducteur de façon que l'élément de surface 

(h devienne da' (fig. 4). Appelons 
^ la projection sur la normale à la 
surface du déplacement rfo- compté 
avec son signe. 

L'augmentation du volume est 
alors 




l'ij?. -.. 



L'augmentation de l'intégrale est 




C- 



- m^ïA^y^m. 




X,ch 



Ceci suppose (jue le potentiel reste constant; mais, par hypo- 
thèse, la figure est en équilibre. Donc il y a un minimum pour 
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— - ; la variation de l'intégrale due à la variation du potentiel 

est donc nulle, puisque cette intégrale est minimum. 
Or, si on appelle [à la densité superficielle, on a 

= — 47C[X, 



Ô/l 



On a donc 



^2 = 1 4^[^%«<y 




= 4^ I ;x-Çrf7. 



Or, il faut que dC=o \ donc 

\k%d(7 == o 

et comme on doit avoir en même temps /ÇrfT = o, quel que 

soit !^, il faut que l'on ait [à = C''', 

On a d'ailleurs M = uiS en appelant S la surface du conduc- 
teur : on a donc 

— — - = 47:^1, 

Supposons maintenant que le conducteur se déforme en restant 
semblable à lui-même : la capacité est proportionnelle à la 
racine cubique du volume, et on a 



d'où on conclut 



dC 

c 




I 

3 


T 


S' 
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Tous les corps dont la capacité est minimum sont donc tels 
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que cette relation existe entre leur capacité, leur volume et leur 
surface. 

Mais alors, pour tous les corps de même volume, la capacité est 
minimum si la surface est minimum ; or, parmi tous les corps 
de même volume, la sphère est celui qui a la surface minimum. 
C'est donc celui dont la capacité électrostatique est mini- 
mum. 

Nous arrivons maintenant au théorème que nous nous étions 
proposé de démontrer, 

La sphère est la seule ûgare d'équilibre. — Je dis que W est 
maximum quand le corps de volume T a la forme d'une sphère. 

Soit T le volume correspondant, donnons-lui un accroisse- 
ment (fîg. 4) 



Comme on a 



W= 




dW= I î^l±^^d. 




V dp dz. 



Le corps étant supposé de densité égale à i, dp est égala — i, 
o ou I ; dp n'est différent de zéro que dans le voisinage de la 
surface, qui est une surface équîpotentielle, le corps étant sup- 
posé en équilibre ; c'est-h-dire que V est constant. 



rfW=V, / dpdx 



LA SPHÈRE EST LA SEULE FIGURE D'ÉQUILIBRE xà3 

Ci rintégrale est égale à dT, On a donc 



0^ 



On peut établir une autre relation entre rfW et dT. 
Supposons que le corps varie en restant semblable à lui-même : 

l'énergie varie en restant proportionnelle à la puissance -"- du 

volume. On a donc 

^W __5_ dl 

\y ~ 6 ï ' 

Eliminant (ifW et dT entre ces deux dernières relations, oi» 

I 

trouve 

W = 4-VJ. 

D'après la définition de W, on voit que Ton peut appeler 
le potentiel moyen. Il est égal à 

Or AV est aussi donné par ï*intégrdle 




S)-(4^)"]* 



étendue a l'espace entier ♦ 

La partie étendue a Teepace intérieur «est, d-nprës la formulé 
de Green (4) -: 

ôV I I 
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Sar la surface, V est constant. L'intégrale a donc pour 
valeur 

^^ rfcr — rf- I VAVrfT, 





d/i 87c 



ce que Ton transforme immédiatement, puisque — AV =t 4'î^o 
= 47r, en 



VAVrfT 




o o 

La partie de l'intégrale étendue au volume intérieur est donc 
a la partie étendue au volume extérieur, dans le rapport cons- 
tant -p-, en supposant évidemment le corps en équilibre. 

Mais l'énergie de ce corps est en raison inverse de sa capa- 
cité électrostatique. 

Nous voyons donc bien que la sphère est le volume pour lequel 
W est maximum. 



LA MASSE EST AMMEE d'uN MOUVEMENT DE ROTATION. 

Formulea générales, — Supposons maintenant le corps, 
homogène, soumis à une rotation w autour d'un axe fixe. 

Appelons J le moment d'inertie autour de cet axe. 
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Si le corps change de forme, de façon que la projection du 
déplacement de l'élément de surface rfa sur la normale soit Ç, 
J varie et on a 



di= I r^d^Çv^ + î/). 



^2 

Nous avons appelé U la fonction \-\ (-^^H-y^)' 

W représente toujours l'énergie due au potentiel new- 
tonien. 
On a 

dVf+ — dJ= I nY+::L^fa:' + y')\d7 



= /^[^+T-(^'' + ^')] 



= 1 UÇrfcr. 



A la surface d'équilibre U est constant et égal à U^. 
On a donc 



co^ 




rfW + -=^ rfj = U, I i;dŒ = VJT. 



D'ailleurs si le corpsse déforme en restant semblable a lui-même, 

5 
W et J varient comme la puissance -^ de T. 



On a donc 



rfW + - 


w* 

2 


di 


w+ 


W^ 


J 



5 dJ 



3 ï ' 



et par suite 



lnT = W + — J=tJ. 
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On a 






AU — 2(0^ + AV, 




2tL)' 4^^> 




2tL)^ 4'^- 



On a déjà vu (page ii) qu'une condition nécessaire pour 
l'équilibre est que 

c'est-a-dire que AU soit négatif. 

Donc U, à rintérieur du corps, ne peut avoir de minimum, et 
comme sur la surface il est constant et égal îi U^, à l'intérieur on 
aU>U,. 

Si AU est nul, U est constant et égal à U^. 

Enfin si AU était positif, U serait dans tout le corps inférieur 

Considérons maintenant l'intégrale 



(O^ 



U£f7=2WH J. 



Suivant que AU est négatif, nul ou positif, le premier membre 
est supérieur, égal ou inférieur à 



On a donc dans les trois cas : 

5 



O \ 2 / 2 



W >a)2J 



AU = o ~ (w 4- -^ j) -= aW 4- -^ J W = w'^J 
AU>o -|-(w + -^j)>2W4--^J W<w^J. 
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Limite de la, vitesse de rotation, — Supposons maintenant 
que co varie d'une façon continue. La figure se déforme d'une 
manière continue, et on a 

d\N+ — di+ cojrfco = 4- Tef U,. 

Lia figure étant supposée en équilibre, W -| J est maximum 

ou minimum, en tous cas 

/2 



rfVV + -^rfJ = o. 

2 



Il reste 



Donc 



tiiidiù — 


5 ''""^ 


7 


>o. 



Ug croit qaand <<> augmente. 

V étant le potentiel newtonien, on a sur la surface 



to' 



V = U,-— (.•r^+./). 

Si donc la surface rencontre Taxe de rotation, U^ est le poten- 
tiel newtonien au pôle. Si 10 croît, ce potentiel croît également. 
Uq ne peut d'ailleurs dépasser 27wR^ (p^^g^ ^6), R étant le rayon de 
la sphère ayant même volume que la masse fluide. 

En divisant membre à membre, on a 

iùidtù d\\ 



w+— J . ^» 



Si on fait croître w indéfiniment, to^ finira par dépasser -. 
W sera inférieur à (o^J ; 



> 



— J+W ^ 
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et par suite 

2 doy d\]^ 



< 



3 0) ^ U, 

Si oj croît indéfiniment, U^ devra croître indéfiniment ; mais 
\.\ ne peut dépasser 27rR^ ; il est donc nécessaire, ou bien que to 
cesse de croître, ou bien que la surface d'équilibre ne rencontre 
plus Taxe ; on a alors une figure annulaire. 

Nous verrons que parmi les figures d'équilibre, il y a des 
ellipsoïdes de révolution ou des ellipsoïdes a trois axes inégaux. 
Pour les premiers, to est inférieur à 4^ X 0,112. Pour les 
seconds, (o < 4?^ X 0,098. 

Il est aussi possible qu'il existe une succession de figures 
d'équilibre, pour lesquelles w passerait par une infinité de maxi- 
mums ou de minimums. Ici, le raisonnement ne s'appliquerait 
plus et lù pourrait augmenter indéfiniment. 

JPermanence de F axe de rotation. — Jusqu'ici nous ne nous 
sommes préoccupés que d'un mouvement de rotation uni- 
forme. 

On pourrait se demander s'il peut y avoir des figures d'équi- 
libre relatif pour une masse fluide animée d'un mouvement de 
rotation non uniforme. 

Nous allons prouver que c'est impossible. 

Nous pouvons d'abord supposer le centre de gravité du corps 
immobile; puisqu'il n'y a aucune force extérieure, son mouve- 
ment est un mouvement rectiligne et uniforme^ 

L'équilibre ne sera pas troublé si on ajoute une nouvelle liai- 
son au corps, si on le solidifie. Alors le mouvement du solide 
autour de son centre de gravité est un mouvement h la Poin- 
sot. 

D'après le principe de Dalembert, le travail virtuel de la force 
totale, pour tout déplacement compatible avec les liaisons, est 
donc nul. Dans le cas qui nous occupe, il n'y a qu'une liaison, 
c'est l'incompressibilité du fluide, s'il s'agît d'un liquide. Si l'on 
a affaire à un gaz en équilibre, l'équilibre subsiste a fortiori en 
supposant le corps incompressible. 

Si Sj:, 5y, Se représentent les déplacements virtuels d'un point 
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.r, y^ Zy la condition d'incompressibilité s'exprime par la rela- 
tion 

rf ^ , rf ^ , d ^ 

ox -\ — oij -| — - — oc = o . 



-N^ 



dx dy "^ dz 

Il peut y avoir trois sortes de déplacements virtuels : 
i^ Les déplacements de toute la masse, considérés comme les 
déplacements d'un corps solide. L'application du principe de 
Dalembert montrerait que le corps doit avoir un mouvement à 
la Poinsot autour du centre de gravité, comme on l'a déjà vu. 

2^ Il y a les déformations du corps. i, ^ /* 

3** On peut supposer que la masse soit animée d'un mouvement ^' /^ 
tel que les molécules se déplaceraient sur des surfaces à densité 
constantes. 

Les surfaces d'égale densité étant aussi des surfaces équipo- 
tentielles, on aurait 

-r^ ox +-V- ^y H — ^ o^ = o- 

^x dy oz 

Les surfaces d'égale densité restant les mêmes, il en est de 
même de la forme extérieure. Tout se réduit donc à des dépla- 
cements internes dont nous allons étudier le travail virtuel. 

Soit un point M de masse w, de coordonnées x^ //, z-, animé 
d'une vitesse x^^y^z', l'accélération du mouvement étant x'\ y" <i^" > 
La force d'inertie a pour composantes mx'\ /«//', mz!' , 

Appelons co^, (o.^, (O3 les composantes de la rotation instantanée 
suivant les axes liés aux corps Ox, 0?/, 0^. 

D'après les formules connues, on a 

/ x' = iù^y — (0,::, 



tOjC WjJT, 



L'accélération a pour composantes 



y" = [iù^z' — 0),.i:') + {l^y\z — iù\x), 
z" = {i0.^x' — iùj/) + {t0\T — Oj\7j). 
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L'équilibre relatif ayant lieu, il y a équilibre entre le travail de 
l'attraclion et la force d'inertie mx", my" ^ mz" . 

Cette force se compose de deux parties : Tune est ce que Ton 

trouverait si Ton supposait le mouvement 
uniforme; l'autre est due à l'accélération 
du mouvement de rotation. 

A un instant /, l'axe de rotation est OP 
((Oj^, (Ojj, (O3) ; à l'instant i+rf/ l'axe instan- 
tané de rotation est OP' (fig. 3). Le vec- 

PP' 

teur OQ équipollent au segment —7 — 

tend, lorsque dt tend vers zéro, vers 
une position limite, et ses projections sur les trois axes sont 

Si on prend cet axe OH lui-même comme axe Oz 

{ù[=z 0)^ = o. 

La force d'inertie due à l'accélération angulaire se réduit ii 
togy, — (03?, o. 

Appliquons maintenant le principe de Dalembert aux déplace- 
ments virtuels que nous considérons. 

Le travail de l'attraction est nul, puisque la forme du corps n'a 




9 
(jj" 



pas changé. Le travail de la force centrifuge est égal à — oJ, en 

posant 0)^ = wj + ^1 + ^s- 

Ce travail est nul puisque J ne varie pas. 
Le travail de la troisième partie est 



^'3 /. m{ijlx — xZy), 



Puisqu'il y a équilibre, il doit être nul. Or, on peut toujours 
choisir un déplacement virtuel tel que 



2^m{7j^j; — xhj) 



ne soît pas nul; il suffît de supposer un courant autour de Oz, 
La somme aurait pour mesure l'aire du plan des xot/, limitée 
par la projection du courant sur ce plan. 
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Il faut donc que (O3 soit nul, c'est-à-dire que Ton ait un mou- 
vement de rotation uniforme. C. Q. F. D. 

11 faut de plus que Taxe de rotation soit un des axes principaux 
d'Inertie ; cela résulte de l'étude du mouvement d'un corps 
solide. 



•i 



Stabilité de l'équilibre,/— Jusqu'ici on ne s'est pas préoccupé 
des conditions de stabilité de l'équilibre. 

Lejeune-Dirichlet a démontré que la condition nécessaire et 






suffisante de la stabilité de l'équilibre est que W + J soit 

maximum. 

Mais il faut faire la distinction entre la stabilité temporaire et 
la stabilité séculaire. Lord Kelvin qui le premier a fait cette dis- 
tinction, appelle stabilité séculaire celle qui a lieu même en 
tenant compte du frottement; tandis que la stabilité temporaire, 
qui subsiste quand on ne tient pas compte du frottement, ne 
subsiste plus si on en tient compte. Lord Kelvin a démontré que : 

La condition nécessaire et suffisante de la stabilité séculaire 
de l'équilibre relatif est la-condilion de Dirichlet. 

Evidemment cette condition est toujours suffisante ; mais sup- 

posons que W H J ne soit pas maximum, il pourra, d'après 

le théorème de L. Kelvin, y avoir équilibre stable, à condition 
qu'il n'y ait pas de frottement. Mais s'il y a un frottement, si petit 
qu'il soit, l'équilibre ne sera plus stable. 

On ne peut appliquer immédiatement le théorème de L. Kelvin, 
car il suppose que tout mouvement détermine un frottement, ce 
qui n'est pas le cas pour cette masse fluide : en effet, si la masse 
isolée dans l'espace se déplace d'un bloc, a la façon d'un corps 
solide, il n'y aura pas de frottement. 

Solide équivalent, — Avant d'aller plus loin, il est nécessaire 
de définir ce que nous aj^pellerons le corps solide équivalent 
à la masse fluide. C'est un solide où, à l'instant considéré, 
les molécules ont la même position que dans le système fluide. 

La vitesse de son centre de gravité" est la même que pour le 
fluide. 



f (j/tM ^/r«^//>i»iv.O- h ' t'-h^u .•>'--.^^ 




. j/^- 
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Les axes principaux d'inertie ayant évidemment la même posi- 
tion, les trois moments de rotation autour de ces axes sont les 
mêmes que pour la masse fluide. 

Son mouvement est donc bien défini a Tinstant considéré ; 
mais il faut remarquer que le solide équivalent a l'instant /, n'est 
pas le solide équivalent a l'instant t' , 

Théorème. — La force vive du fluide est équivalente à la force 
vive du solide équivalent^ augmentée de la force vive du fluide 
dans le mouvement relatif par rapport à des axes invariablement 
liés au solide équivalent. 

En effet, appelons Ç, y|, J^, les composantes selon ces trois axes 
de la vitesse absolue de la molécule m ; Ç', tj', Ç' les composantes 
de la vitesse de la molécule correspondante du solide équiva- 
lent ; ^\ r/', 'Q' les composantes de la vitesse relative de cette 
molécule par rapport au solide équivalent. On a 

1 = 1' + 5", 

ç = ç' + !;". 

Les trois forces vives considérées ont pour expression : 






et il faut démontrer que 

T = T' + T". 

On a 









CONDITION DE STABILITÉ DE LEJEUNE-DIRICllLET 

Je vais démontrer que 



2 /H(i'$" + fS" + Ç'C') = O. 



Cette expression représente le travail d'un système de force S", 
lel que la force agissant sur la molécule m ait pour composantes ^ 
mW /wr/', mQ'j mais ce système est la différence de deux autres S.\ç^ O 
et S', dont les forces sont ni^y m'r\y mX^^ et /wi', mr/, mXl ; ces deux ^->À 
systèmes de forces sont équivalents. La résultante générale est la v:^ '-^, 
même et les moments résultants des deux systèmes sont les . "^ 
mêmes par rapport aux trois axes, par la définition même du Vs^^ '^ 
solide équivalent. ^ ^ '* 

On a en effet T $i 

puisque la vitesse du centre de gravité du solide équivalent est la ^ 
même que celle de la masse fluide. 
De même on a 

puisque les moments de rotation sont les mêmes autour des trois 
axes de coordonnées. Donc S = S' et S'' est nul; son travail Test 
aussi. 

Par suite 



^m(ir+Vr/' + C^")=o. 



Condition de stabilité de Lejeune-Diricblet. — Revenons 
maintenant a notre masse fluide. 
On a Tégalité 

T — W = C\ 

D'où 

T'+T" — W = 0% 

Cette égalité suppose qu'il n'y a pas de frottement : autre- 
ment T'-f-T'^ — W ira constamment en diminuant. 

PoiNCABÉ. — Figures d'équilibre. 3 
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Si on pose 

T' + T" — W = £, 

—7- est nul s'il n'y a pas de frottement. 

Mais alors c'est qu'il n'y a pas de mouvement relatif, et par 
suite T'' = o. 

Si T">o, ^<o. 

Dans ce mouvement de rotation le théorème des aires s'applique, 
et on peut considérer comme des données les trois moments de 
rotation. 

Je vais démontrer que la condition nécessaire et suffisante de 
stabilité est que T' — W soit minimum. 

Pour une position d'équilibre relatif, on aura 

X' T' 

^ ^ 

W = Wo 

T" = o. 

Si l'on part de cette position d'équilibre, et que l'on fasse varier 
très peu les éléments des mouvements des molécules, mais de 
façon que les moments de rotation ne soient pas altérés (c'est 
une restriction sur laquelle nous reviendrons), la masse fluide 
ne s'écartera pas beaucoup de la position d'équilibre. 

En effet, soient T' et W, les valeurs de la force vive et de 
l'énergie, on a 

T' — w > T'„ - W,. 

On a fl fortiori 

T + T' — ^N>T, — ^y,^ 

T-fT'' — W = e est d'ailleurs voisin de T^— W^. 

Cette quantité ne peut jamais croître, donc elle restera toujours 
voisine de Tq — Wq. 

Les valeurs des variables ne différeront que très peu de leurs 
valeurs pour T^ et Wo, donc il y a équilibre stable : la condi- 
tion est suffisante. 
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Elle est aussi nécessaîre. En effet, si T^ — W*o n'est pas mîni- 
mum, je puis choisir T' et W de manière que 



et même que 



T' — W<T'„ — W„ 



T'+T''— W<T' — Wo, 



T" étant suffisamment petit. 

Mais alors T'-f- ^" — W ne peut que décroître ; T'' lui, tendra 
vers zéro; T' et W ne pourront alors reprendre leurs valeurs pri- 
mitives. 

Il y a une restriction à lever. Supposons que Ton écarte le 
corps de sa position d'équilibre sans conserver la valeur des 
moments de rotation autour des axes principaux ; on sait alors 
que le solide équivalent aurait un mouvement où les moments, de 
rotation auraient des valeurs peu différentes de celles qu'ils 
avaient précédemment et l'on est ramené au premier cas. 



L'axe de rotation est le petit axe de Vellipsoïde d'inertie, 
— Je vais démontrer maintenant que l'équilibre n'est stable que 
si l'axe de rotation est le 
plus petit axe de l'ellip- 
soïde d'inertie relatif h 
la masse fluide. 

Soient (fig. 6) Oxyz un 
tricdre fixe de coordon- 
nées ; OA,OB, OC, les axes 
principaux d'inertie; w la 
vitesse instantanée de rota- 
tion dont les composantes 
sont (1)^, a>2, a>3 suivant les 
trois axes fixes,/?, y,/' sui- 
vant les axes OA, OB, 
OC ; on a Fig. 6. 




.2 



Soit OR l'axe instantané de rotation, J le moment d'inertie 
par rapport à cet axe. 
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On a 

Soient [JLj, [i.g, [Xg, les moments de rotation par rapport aux trois 
axes de coordonnées. Le segment OM, dont les projections sont 
(l^i> \^29 l^s)» ®^* ^^^ dans l'espace (théorème des aires), et l'on a 

Les projections de OM sur les axes OABC sont A/>, B^r, Cr, 
et l'on a encore 

JJL^=Ay + By+cv^ 

La condition de stabilité, comme nous l'avons dit, est que 
T' — W soit minimum. Or, si on change l'orientation sans 
changer la forme du corps, W ne variera pas. Donc il faudra que 
T' soit minimum. OM est fixe dans l'espace. 

Si on pose 

p v/A = X, 

q V/B"= Y, 

/• y/c = Z. 

On a 

T' = — (X^ + Y^ + Z^) 

[A^ = (AX^ +BY^ + CZ2) = C^ 

Donc, le point X, Y, Z doit être sur un ellipsoïde donné et sur 
une sphère de rayon minimum. Donc Taxe instantané de rota- 
tion est le plus petit axe de l'ellipsoïde donné et par suite, de 
l'ellipsoïde d'inertie. 

Si A est la plus grande des quantités A, B, C, l'axe de rota- 
tion sera l'axe des x» Alors on a/j = to, y = o, /•=o, 

A = J, 



2J • 
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Il faut donc que 

2 



(0 W — T' = W — 



I* 



2J 

soit maximum ; donc que 



5(w_4)=« 



Il suffit aussi que i .^ 

soit maximum (p. lo), condition qui n'est pas nécessaire. 
Donc 

.2 



s/w + -^j) = o. 



Les deux conditions donnent : 



(-^ ) 8w + — ôj = o. ^' • / '-r ^'^~ ^"^ 

On doit donc prévoir que si la condition (2) est remplie, la 
première le sera également. 

Supposons qu'il en soit autrement, c'est-à-dire que l'on puisse 
avoir en même temps 

w + 4j<w.+4j. •' 

2 2 



w — -.->w„ — 



2J " ° aJ • 



Comme on a u, = (uJg, on conclurait facilement 



"2 2 J 



'0 






in' 

'4 



a ij^- 'l-HKK'h 



(I 



,1* «JJ, I .1/.' 0, 

(4 .!„)*. 0. 



Ml l|m (1*1 itltollHlt* I litliU lu l'<Vl|«i'0(|uo uViit pua vruio ; lu 

(l«>M^I(i()m mMulUtiKt ii't>hU'uliit« |Mt» Ih |tnMui6i-(>, 



^V ^^-'^ ^inrtf^ — /^ 






o 



^ >y 



>, 
v^ 
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FONCTIONS SFHÉRIQUES 



Expression de AV dans un système triple orthogonal. — 
Soient 

Ç(^, y, -s) = y, 

trois familles de surfaces, constituant un système triple ortho- 
gonal. 

Considérons le petit parallélépipède rectangle dont les faces 
sont les 6 surfaces. 

Les arêtes de ce parallélépipède ont pour longueurs 

ad^, hdy\^ cd^, 

La diagonale est 

ds^ = a^d^' + iVr;^ + 6•Vî;^ 

Supposons une fonction V (Ç, vi, 2^) et cherchons à déterminer 
la valeur de T'intégrale 



/ 



AVrfT 
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étendue au volume du parallélépipède; elle est égale à 




étendue a sa surface. 

Le volume étant înfinîment petit, la première intégrale a pour 
expression 

La seconde est la somme de trois intégrales étendues chacune 
aux deux faces opposées du volume. 
Pour la face \ =^\^^ ou 2i 

Pour la face opposée 5 = ?o -h^ly on a 

et par suite pour l'ensemble des deux faces, on a 

On a donc 

ûiYabcd^ d'r\ rfî^ 

jrij jyV ^ hc ^Y ^ d ca dV . d ab dV 1 



et finalement on a 

AV 



~ abc Zj%^ \"â" di / 



Application aux coordonnées polaires, — Dans ce système 

on a 

X =^v sin 8 cos o, 

7^ =/• sin sin cp, 

^ == /• cos 9, 
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et, en posant cos 8 = ja, on a 

a: = /• y I — [jl^ cos 'f , 
y =ry/i_jx^ sin cp, 

Z = /-[JL. 

L'élément linéaire dans ce système de coordonnées a pour 
expression 

ds^ = dx^ + dtf + rf^"^ = dr^ + /-Vô' + /•' sîn^Qtf'f '^ 

= rf;.2 + ;.2 _^1^_|_ /-(l — a2)t;C0^ 

I [A'* ^ *^ ^ * 

Au lieu de Ç, y;, ?^, nous avons /•, [jl, '^ et dans ce système 

a= I 

h = 



9 

r 



c = r\/i — [JL% 
abc = r^. 



On calcule alors facilement 






I d^V 



Polynômes sphériques. Fonctions sphériques. — On appelle 
polynôme sphérique un polynôme homogène en x^ y, z satisfaisant 
à l'équation de Laplace, A P = o, 

b il est de degré n en jt, t/, z il contient ^^ -^ coelti- 

cients arbitraires. Mais il y a un certain nombre de relations 
entre ces coefficients puisque l'on a AP„ = o. Cette relation est 

de degré n — 2. Il y a donc n ' — • relations entre les coef- 
ficients. 
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Il y a donc 

(/i+ i)(/i-f- 2) n[n — i) 



2/1+ I 



polynômes sphérîques, indépendants du degré /i. 

On peut exprimer a;, y^ z en coordonnées polaires : on aura 
alors 

p ,.n Y 

Y„ étant uniquement une fonction de y et de [jl, que nous 
"appelons fonction sphérique d'ordre n. 
On a évidemment la relation 

(i) r^^^n?n. 

Considérons maintenant l'égalité (2) démontrée page 4 

Supposons que la surface S soit une sphère de rayon i , que U 
et V soient deux polynômes sphériques d'ordre m et /i, le second 
membre sera nul, et on aura 



p,. ^f\d<' = o. 



Doue, ea tenant compte de l'égalité (i), on a 




-">/ 



{m — n) I P„, P„ d<T = o. 



Si m est différent de n, on a 




PjR Pn d<7 = o, 



poirNOMES sPOÉRiquES 43 

rintégrale étant étendue à une sphère de rayon i ; et on peut 
encore écrire ceci : 



[J: 



Si m = 71, ce résultat nVst plus vrai. 

Fonctions spbèriques fondamentales. — On a vu que Ton a 



a; = /• y I — tx^ cos cp = r y] 



' — r — 



y = r y I — \k^ sin © = /• y i — 



L* 

7,1 



Si on remplace x q\ ij par ces valeurs dans P„, on aura 






Xj étant une fonction de [jl seulement 5 ce sera un polynôme 

en [JL et en ^ i — [JU. Les polynômes où/? est pair, ne contiendront 

y I — \k^ qu'à des puissances paires. Les polynômes où p est 

impair, ne contiendront y i — [x^ qu'à des puissances impaires. 
Les premiers seront donc des polynômes en [jl, les seconds 

seront des polynômes en \k multipliés par y i — [jl^. 
Substituons la valeur de P dans Téquation A P„ = o. 
On obtient 

- I o. 

I — a'' J 

En annulant le coefficient de /"-^ eO>?, on obtient Téquation 
à laquelle doit satisfaire Xfj. C'est 

jLr(._t.^)^1- .(.+.-)(.^--0+/.- X;; = o. 

rf[x L^ ^^ tfix J '— [^ . 
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L'intégrale générale de cette équation est une fonction trans- 
cendante, mais il y a une intégrale particulière qui est de l'une 
des deux formes exigées par la définition de Xj, et qui est 
déterminée à un facteur constant près. 

On voit que XJJ et X^^ satisfont à la même équation. On a 
donc 

A étant une constante quelconque. On a ainsi la forme des 
2/? -{- I fonctions sphériques cherchées. 

<«v^ • *^f^ t/ • , ^^n s • • • • î ^^n ' * * • 1 

On peut écrire autrement 2/1+1 fonctions indépendantes en 

posant 

XJ , X^cosp, XJcos2cp ,..., X"COSW'^, 

Xi sin cp, X^ sin 2'^ , . . . , X" sin fVf, 

Deux fonctions sphériques différent donc soit par /i, soit par p^ 
soit que n et p soient égaux, mais que l'une contienne le sinus, 
l'autre le cosinus de p'f. 

On a déjà démontré que 




est nul si mr=f=.n\ nous allons achever la démonstration et démon- 
trer que l'intégrale 

est nulle si Y et Y' sont deux fonctions difiérentes. 

Il suffit de le démontrer dans le cas de deux fonctions de 
même ordre. 



On a alors 




+1 

XJ XJ d]k I 00s p'^ cos y!p d'^ 

1 J ^ 
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OU bien 



X^ Xj d]k I sîn jy^ sin q'^ (f'f , 

-1 .7 

X«^ ^'[J^ / smjrfcosp'^do, 

-1 .7 




Dans les trois cas, la seconde intégrale est nulle. 
Enfin il est évident que 




est différent de zéro. 

Détermination des polynômes X{|. — Proposons-nous main- 
tenant de déterminer Xf^y et pour cela, commençons par le cas 
où /; = o. 

L'équation est alors 

C'est l'équation qui définit les polynômes de Legendre, et 
on a 

" dix'' ^ ^ ^ 

Xj (Cos 6) est donc une fonction sphérique où 9 représente 
l'angle d'une direction avec l'axe O^. 

Par raison de symétrie, Xj, (Cos e), où z représente l'angle de 
la direction considérée avec une autre direction quelconque, mais 
fixe, sera aussi une fonction sphérique ; mais on a 

cos e = cos 9 cos 0'+ sin sin Q' cos (-^ — '^'j 
et cela, quels que soient 9 et 9'. 



46 FONCTIONS SPHÉlilQUES 

Donc 

X^ I cos 6 cos 6' + sin 9 sin 8' cos {f — '/) 1 

sera encore une fonction sphérique. 
Si on pose 

sin ft' e-'î' = 2$, 

sin 8'e*?'= 271, 

d'où on conclut, en ajoutant membre à membre après avoir mul- 
tiplié ces deux équations respectivement par e*? et e~ *?, 

sin 6' cos (cp — o') = $e'? +7^6"'? 

et en multipliant membre à membre, 

sin^9'=4ivi; 

d'où 

cos^ô' = I — 4i'i ; 

la fonction 

X« [v/i — 4i7i [X— \/i— [x'-^ (?e'? +r.e-^>) ] 

est donc encore une fonction sphérique quels que soient ^ et r». 
En particulier, si on fait 

"1= o, 



X„[fx + v/T^^e>] 



est une fonction sphérique. 

Si je développe suivant la formule de Taylor, j'ai un poly- 
nôme : 



Xn [[X+v/l— [X^e»>]= Xn([x)+6^>V/^ 



Y" ^n 



d^ 



■ I ■ ■ i^ w 1.^ — I — ■■ ■ ■■ I ■ I » ■ ' i -^ I ^ ■ I I ■ I ■■-■!. . .^ ■■ . 

1,2.., JJ dltJ' '" I.2.../Î d[L^ 
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Chacun des termes de la somme est une fonction sphérique. 
On a donc défini les fonctions X^ à un facteur près. 



x; = (,-,-)-5-i^ 



Propriétés des fonctions spbériques. — Soit U une fonction 
de r, et P„ un polynôme sphérique : 

Proposons-nous de calculer A (U I\). On a 

,^ dU dP„ 



mais 



A(UP„)=P„AU + .^-^^-^+UAP,., 



' dar dr ' dx dr 



r 

d^u 2 cm 



MUP,.) = P«[-^ + -^J 

2 rfU / dP,. dP„ dP,. \ 

■+" ~ "rfT r "dJ~ +^ ^ "^^ "ôT j 

et, parce que P" est une fonction homogène, 

_p r rfnj , 2(/i+i) rfu-1 

Proposons-nous maintenant de calculer 

^(uy„)=a[^xp,.] 

puisque P„ = /•"¥„. 
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p'oscTioss fit^iii^niQVEs 



Oôtiche êptiêVOÏdAÎe. — Noun nlloni* npplîqucr ce qui précède à 
l'élude tl'utie uirtftfte fluide nphéroïdrtle, c'est-a-dire d\me masse 
lluldê dont lit flguie d'ét(uilil)ro orI peu diflerente d'une sphère, 
uil, d'une utnttiîMe pluR pri^otse, dout lu surface est comprise entre 
deu^ sinTrtocR RpluMîquea de rayons voisins. Une couche sphéroï- 
dide seiu Uhe couche comprise cuire deux surfaces sphéroïdales. 

Action tt^nM tonthÉ ëphétoMAl^. — Soît (fig. 7) r/^rêlément 
de ftuH'rtce d'une couche de rayon 1^ € IVpaisseur de la couche 

Rphi^ix^Vdale, rft rclcment de volume 
de la couche sphéroîdale ; on a 

On peut snpposct que la masse 
<^ontenwe dans d'z est sur rélément 
/i(Ty et e«t représentée par ixdfs^ 

Soit un point M dont les coordon^ 

nées C4ïrtésiennes sont .?', î/, r^ les 
coordonnées sphériques 

7\ tJL = COsO, Ci, 







>ioil un point M' dont les éléments analogues sont ,r', t/\ r'; 
7 ' == I , Vi.' =<y>s h'. «'. 

La di^tam^e MM e^t égale à i + 7 - — p. 7- cos ^OM, OM' . 

<]os "' = cos fOM, OM " = cos h cos h' + sin h sin 6 cos ;t — :t'i . 

Développement en Tsérie de-^^.-. — Le pot<Mifiol en IVI, du à 
la couche sphéroïdale, esl 



V=. 




t' esl une fonction de u el de cp : ;• es! phîs polil qn(* i. 
Or 



MM''' = I — srcos *'-i- ;■- =1^ / 1 — /v'. ' I - ;v ',, 



rrrr' I VC^^ 

MM' V 



I — fC 



-7- 
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On a 



/ \ J_ • ï I 3 

( I — /•e'Y)"' 2 = I H /e'v-^ 1— /•V''+... 

\ / 2 2.4 . 

. 1.3.5... (in — i) . . 
H TTr ''"^"'"^ + • • • 

2.4.0... 2/2 

C'est une série absolument convergente pour /•< i. 
De même 

/ .\ L .1 . , . 1.3... (2/1 — i) .• . 

\ / 2 2.4... 2/1 

Le produit de ces deux séries est une série encore conver- 
gente, si /• < I ; le coefficient du terme en /•** est 

. 1.3... f2/i — . 1.3... (2/1 — 3) I , ,,. 

2.4 .. 2/1 2.4... (2/1 2) 2 

, 1.3... (2/1 — 5) 1.3 ,, ... , , 1.3.3... (2/1 — i) 

2.4... (2/1 4} 2.4 2.4.0... 2/1. 

Le coefficient du terme en é*'^ égale celui du terme en c"^"*. 

On a donc un polynôme à coefficients réels du cosinus de 
l'angle y et de ses multiples jus^[u'à n ; comme les coefficients de 
ce polynôme sont tous positifs, ce polynôme atteindra sa valeur 
maximum pour y== o. Mais alors : 



MM' I— /• 

Donc les coefficients du développement sont inférieurs ou égaux 
h l'unité. 

On peut donc écrire 



IM' Zj " ' 



MM' 
A„ étant une fonction de (x, de », de a' et de o'. 

Développement en série du potentiel d'une couche spbéroï- 
dale. — Le potentiel V sera alors égal à 



2,'-" / A„eW«r. 



PoiNCARÉ. — Figures d'équilibre. 
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Les termes de V sont moindres que /•" I e'rfo". 
L'erreur commise quand on prend pour V la somme des n pre- 
miers termes, est moindre que 

La série sera donc convergente si r<i. 

^--^ peut être développé suivant les puissances croissantes 
des variables x^ y, z, 




mM 



y/(y _ ^)^+ (y _ yf + (^/ _ zf 



^x^yz'^ I 



Kipq ^ài! 



et, si on revient aux coordonnées polaires, on aura un développe- 
ment suivant les puissances de /•". 

Si j'appelle V„ l'ensemble des termes de degré n dans V, j'au- 
rai 

Or V satisfait a l'équation de Laplace : on a donc 

AVo + AV, + ... + AVn + ...=o. 

Cela exige que chacun des polynômes V„ soit un polynôme 
sphérique; et, si on remplace dans V„, x^ y, z^ en fonction de r,[jL 
et çp, on aura un résultat de la forme /•" Y^j où Y„ est une fonction 
sphérique. On a donc 



=i-^.- 



Cette formule est valable si r est plus petit que i. 
Supposons maintenant r supérieur à i, et considérons le 
point M^ image du point M, tel que OM^ X OM =: i (fig. 8} ; 
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DI 



j'ai OM = i\ OM^ =— - . Si M' est un point de la sphère, on a 
les égalités 

MM' OM OM' v^OMxOM' 



M'M, OM' OM. -^ÔM^OM 
Or le potentiel en Mj est 



-M' _ . /O.M 



/ t'dit / t'dt 

M,M ~'' / MM' ■ 



On 


a donc 




Or 


le potentiel en M,, 


, point 


situé 


à une distance 


' du 
r 


centre, est 






Zj /•» 






M 



Fig. 8. 



Le potentiel en M est donc 






Supposons maintenant que M et Mj se rapprochent de la 
surface de la sphère en restant inverses l'un de Tautre. Les 
valeurs du potentiel M et M. sont voisines, mais les dérivées 
ont des valeurs difTérentes de part et d'autre de la surface et 
on a 

dV dV, 



mais 



d/-. 



V. = Vr, 
dV 



r dV ,,n dr 



= — V/-'— 



ÔV 
d/- 



,s 
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On a donc lorsque /• tend vers i 

Dans les mêmes conditions, puisque -^r — ■ est fini, on ii 

11 m -^:^ — ^,r — /•') = o; 
Or ^ ' 

enfin 

lim V (i — /•^)^o. 

Si on ajoute ces quatre égalités, on voit que l'on a 



Posons 



liml s/*-— — t-V l=47^£. 



Or 



Si r est inférieur a i, 4> peut se représenter par une série de 
la forme 



<ï>=^(2/e+i)/-Y„; 



<ï> est donc une fonction satisfaisant îi l'équation de Laplace, 
puisque c'est une somme de polynômes sphériques. 

Si r tend vers i, ^ tend vers ^ize. Si on démontre que ce 
développement est encore valable pour /*== i, on arrive au théo- 
rème de Laplace : 

Théorème. — Une fonction quelconque de deux variables est 
développable en une somme de fonctions sphériques. 

Potentiel en un point de la couche sphèroïdale. — Soit e la 
fonction donnée sur une sphère de rayon i. 
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Elle définit a Tintérieur de la sphère un potentiel V, prenant 
sur la sphère la valeur ^tzz. 

Ce potentiel est développable en une série entière 



(0 2^A„' 



Jfl 



convergente si r est inférieur à i ; et, si on regarde r comme une 
variable complexe, cette série définit une fonction ^ convergente 
h l'intérieur du cercle de rayon i, et qui sera obtenue en faisant 
dans V, r^pe*?. 

Dans V, faisons r = e*»; on a une fonction 4>' définie sur le 
cercle de rayon i, sauf peut-être en un certain nombre de points 
singuliers. Si la série (i) converge pour p = i, elle représentera 
4>' d'après le théorème d'Abel. 

Je dis qu'il en est ainsi dans le cas où ^' est développable en 
série de Fourier, 

Nous avons alors 

Je dis que Ton a 

A;,=tAn . B,; = o. 

En effet, d'après le théorème de Gauchy, la suite des valeurs 
^' définit à l'intérieur du cercle de rayon i une fonction qui 
n'est autre que la fonction V (p^*-). 

Cette fonction est développable suivant la série (i), et on a 
d'autre part 

les deux intégrales étant prises le long du cercle de rayon i . 
Mais la première intégrale se réduit à A',», la seconde à B'„. On 
a donc 

B;=o. 

Donc, si ^' est développable en série de Fourier, la série sera 
encore convergente pour'r=;e'? et en particulier pour r= i* 
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Tout revient donc à démontrer que 4>' est développable en 
série de Fourier, c'est-à-dire que <ï>' est une fonction de o admet- 
tant une dérivée, sauf en quelques points singuliers, et qu'en ces 

points singuliers l'intégrale / 1 <Ï>M rf:) a un sens. 

Nous allons maintenant démontrer que e est développable en 
une somme de fonctions sphériques. Le potentiel en un point M 
de coordonnées (r, 6, '^) est 



MM 



s' étant la densité de la couche sphérique, au point (6', »'), 
On a 

dV / t'd^' dMW 



dV__ / t'd^' 
dr ~ I MM'' 



dr 



l e' (ces Y — '■) dv' 
~ MM^^ 

Comme cela a été démontré, nous avons 



:,itt=^(^.<o)= 



Supposons que OM soit une direction fixe que nous allons 
prendre comme axe O^ : 

e = o, [x=i, 

MM' =1 — 2/' cos Y + /'^ ^ I — ajjiV + r^. 

L'élément de surface est di' ==rf[jLW'f', et on a pour valeur en M 
de la fonction 4> 




rfj/.' 



_.l (l— 27y+/-2) 



2 
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car e' seul dépend de ^', Si on pose 

F ([.')= / edi, 
on a 

<ï>=(l — /-) I 

t/-i (i — 27-[À' + r2) 



2^ 

2 



Supposons que la sphère soit partagée en un nombre fini de 
polygones sphériques dont les côtés soient des arcs de courbes 
analytiques, et que dans chacun de ces polygones e' soit une fonc- 
tion continue de ^' et de ip '. 

Dans ces conditions, F([jl') sera une fonction analytique de [x', 

excepté pour les valeurs de ^I correspondantes aux sommets des 

polygones, ou aux parallèles de la sphère tangents aux côtés des 

F (aO 
polygones. La fonction ' ^^ sera alors intégrable entre — i et 

Si on pose r=:e*^, on aura 

MM'= 27- — f/'H — ^j — ^' =2/*| cosij/ — ^'\ 

et si on pose ensuite 

cosij; — \f''= Z"^» 
on aura 

^"' F (fxO d^! 




V 

-1 

Nous démontrerons plus loin que ^ est une fonction de ^ 
admettant une dérivée, et est par suite développablê en une série 
de Fourier. 

Enfin, si nous développons suivant les puissances de /*, nous 
aurons une série convergente pour|7*|<i. Elle sera sûrement 
convergente même si 1 7* | = i ,"dans le cas où s né comprend qu'un 
nombre fini de termes. 
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Voyons quelques exemples qui nous serviront par la suite. 

Si on fait e' = , on aura 










si on fait e' == —^ , on aura 



On pourra prendre pour fonctions sphériqucs correspondantes 
Y.=-î-ouY„ = — . 

De môme, si e' = —^, on a en appliquant la formule de la page 52 



47:e' = Y.+3Y, + 5Y,+ ... 
Dans le cas qui nous occupe 



Yj, = Y^ = . . . = o 

et 

^-=2/' = 2e'K 

Théorème de Laplace. — On a toujours le droit d'écrire 

a et i étant des constantes quelconques; 4> sera alors de la forme 

2a + 2ie'^ + <ï>^ 
et si on pose 



, ,v I — r^ — isintL . * 
, 2/rr v^A 



on aura 



4.,=-/l(^) 




THÉORÈME DE LAPLACE 5; 

Ceci posé, pour établîr que 4> est développable en série de 
Fourier, il faut démontrer : i° que --tt- existe et est fini sauf en 
un nombre limité de points singuliers de la fonction 4> ; 2** que 

1^1 rf^ est fini en ces points singuliers. Il y a une difficulté, 

car Z s'annule si [jl' = cos i{/. Il faut donc que F ([jl') s'annule pour 
jjL' = ebs^j^ et alors on doit étudier vers quelle limite tend la 
fonction lorsque [jl tend vers cos ^. 

I** — 7î— existe. En effet on a 

d^ ,.,. . rfvi / ^^d\J' ,3 . , / V,du! 





d^ 4 / Z« ' 2 ' ^ } Z^ 

— 1 

Or, je puis disposer de a et i, de façon que Ton ait 

Fj (cos ij;) = o, 

F( (cos if) = o. 

Si donc cos if — \k' est un infiniment petit du premier ordre, 
F(jjl'), qui est analytique au voisinage de cos if, sera un infini- 
ment petit du second ordre, et, sous le signe d'intégration, les 

3 5 
dénominateurs sont respectivement d*ordre — et — ; le numé- 

rateur est d'ordre 2. L'intégrale est donc finie, sauf si cos if est 
une valeur singulière de F([a'). Il faut qu'en ces points singuliers, 
où i|; == if^, l'intégrale 




soit finie. 
Supposons 

a = F (cos if^ 
b = 0, 
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On a alors 



14>M.{-<2 / lF(cos;}/)|rf.}+ / T,d-l I liil^l^i£2îilirf^, 



Le second membre doît rester fini. La première intégrale est 
finie. Il faut démontrer qu'il en est de même de la seconde. 

C'est une intégrale double qui, aux points singuliers, doit être 
un infiniment petit d'ordre inférieur à 2, et aux lignes singulières 
d'ordre inférieur à i . 

Or, nous avons ici une ligne singulière, la ligne [x' = cos ^, 

3 
mais pour [/.' = cos ^^ le dénominateur est d'ordre — > le numéra- 

teur est d'ordre i, sauf aux points singuliers de F(ul'). Donc 

l'élément de l'intégrale devient infini de l'ordre — , l'intégrale 

a un sens. 

Aux points singuliers F([jl') =^ F(cos ^), la fonction est d'ordre 

3 
— , l'intégrale a encore un sens. 

La fonction <t> remplit donc les conditions nécessaires pour 
être développable en série de Fourier. 

La formule de Laplace est valable pour /•=: i, et Ton a 



L-^{^n+i)\\. 



Le théorème de Laplace est ainsi démontré. 



CHAPITRE IV 
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Préliminaires. — Supposons le mouvement de rotation de la 
masse suffisamment lent : la surface différera peu de celle d'une 
sphère, le carré de la vitesse de rotation o)^ sera considéré 
comme un Infiniment petit du premier ordre. 

On a déjà vu que si le mouvement de rotation est nul, les 
surfaces équipotentlelles sont également des surfaces d'égale 
densité. Le potentiel V^ en un point est uniquement fonction 
de /• et Tequatlon de Laplace se réduit à 

(0 ^^o = -^ + -^ = -4«?. 

La masse contenue entre deux sphères Infiniment voisines de 
rayon r et r -j- dr est 

et la masse totale contenue à l'intérieur de la sphère de rayon /• 
est 



Tzr^^dr, 



Si on appelle D la densité moyenne de cette sphère 
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D est une fonction de r dont la dérivée D' est négative dans 
l'hypothèse d'un équilibre stable. 

En difTérentiant cette dernière équation, on a 

o 
En supprimant le facteur ^Tzr^dr^ on a 

(2) 3p = 3D + /D'. 

En prenant encore une fois la dérivée, on obtient la rela- 
tion 

(3) 3p' = 4D'+7D" 

qui nous sera utile par la suite. 

Faisons encore la remarque suivante. On a évidemment : 

3D^3D + /D'^o. 

La première inégalité est vérifiée puisque D' est négatif; la 
seconde Test, car p est positif. On a donc 

On a • - =0 si la masse est homogène. 

On a -r=-- = — 3 si toute la masse est concentrée au centre. 

Le potentiel de la sphère de rayon /• en un point de sa surface 
est 

o 
La force qui s'exerce sur une molécule de la surface est 

Puisque p et V^ sont uniquement fonction de /•, on peut les sup- 
poser fonctions l'un de l'autre, et, en vertu de la relation (i), 



DEVELOPPEMENT DU POTEr\TlEL EN SERIE 6l 

AVj, est fonction de V^, et on peut écrire A V^^ =/*(¥,,) ; en diffé- 
rentiant 



df{N,) rf(AV,) dr - 4irp^ 

TT/D 



^V, d\, dY, 4 



rf/- '6 

En tenant compte delà relation (3), on a 

4D' + /D" 



(5) /' (Vo) = 



/•D 



Développement du potentiel en série, — Supposons main- 
tenant un mouvement de rotation lent. La surface se déformera, 
mais peu. Le potentiel en un point d'une sphère sera dévelop- 
pable en une somme de fonctions sphériques, et Ton aura 

Y représentant une fonction sphérique et H étant un coeffi- 
cient qui est une fonction de /• de Tordre de w*, puisque V dif- 
fère peu de Vq. 

L'équation de Laplace peut s'écrire 

(6) AV = AV. +2 A (UY) = - 4^:? . 

On a d'ailleurs, d'après une formule démontrée plus haut 

(P- 47) : 

en supposant que n soit l'ordre de la fonction sphérique Y. 
Si on pose 

U=V+ — (.r^+y^), 



la condition d'équilibre est, comme on Ta démontré plus haut, 
que les surfaces U = constante coïncident avec les surfaces 
d'égale densité. 
On a donc 

AV = — 4^^? = fonction de U, 



62 MASSE FLUIDE HÉTÉROGÈNE. PROBLÈME DE CLAIRAVT 

(0^ étant très petit, V difTère peu de V^ ; 

AV diffère peu de AV^. 

La relation entre AV et V diffère peu de la relation entre AYy 
et Vq, et on peut par suite écrire 

AV=/-(U) + o(U), 

c3 étant une fonction dont la valeur est constamment de l'ordre 
de tù^. 

Si l'on développe le second membre suivant la formule de 
Taylor, on a 

AV = /-(V.) + (U- V,) /'(¥,)+.. . 

+ ?(Vo)+(u-\g'f'{v.)+... 

En vertu des hypothèses faites, on peut négliger les termes 
infiniment petits du deuxième ordre et écrire 

AV =/-(V,) + (U - V,) f (V.) + ? (V„) 
= AV. + (U-V,)/'(V„) + ?(V„). 

Enfin, en tenant compte de (6), on peut écrire 

2)a (HY) = (U_V,) /'(V„) + 'f (V„). 
On a d'ailleurs 



iù^ 



(x^-^-y^) est lui-même décomposablc en une somme de 
fonctions sphériques 

w*/'* w'/'^ x^ H- ;/^ — 2z^ 



3 ' 6 r^ 
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Le développement se réduit h la somme des deux fonctions 
sphériques Y = i et Y = ^^ affectées respective- 

ment des coefficients C, =- — r. — etC, = — - — • 

^ o -* 6 

Les autres coefficients sont nuls. 

On peut dans ces conditions écrire 

U-V,=2(H + C)Y, 

2)a (HY)==^(II + C)Y/'{V„) + »(V.). 

Détermination des coefûcients du développement. — Pour 
trouver les valeurs de II, on peut égaler terme à terme les coef- 
cients d'une même fonction Y. 

A(HY) 
contient Y en facteur ; on a donc 

A(UV)=Y(IH-C)/'(V.) 

Si on prend Y= i, on a 

(8) Mi=(H+^,..):^:+«:+,(vj, 

H est donc déterminé par une équation linéaire du second 
ordre. 

Si on prend 

a;* + y* — 2z^ 



Y = 



/•^ 



H est déterminé par l'équation : 

wV*\-, /•D"+4D' 



(9) A(HY)=(h+^)y 



/•D 



« 9 
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Pour les autres, on a simplement Téquation 



(lo) A(IIY) = IIY 



/•D 



Ces dernières équations ont une intégrale évidente, H==o; 
nous allons démontrer que Toii ne doit pas en considérer d'autres. 

Dans le cas où' Y est une fonction du premier ordre, on peut 
trouver la valeur correspondante de II. En efiet, il y a trois fonc- 
tions indépendantes du premier ordre : 

X y z 



j, > /• ' /• 



Considérons Tune d'elles, -^ par exemple, on a vu que Ton 

avait 

AV„=/-(V,). 

Si je prends les dérivées des deux membres par rapport a ;r, 
j ai : 



^i^hriy.) 



ùV. 



\ dr ■ iix)~' ^ "' dr • dx' 

.[dV, X ] dy„ .T rW + ^lù' 
L dr ' r J dr ' r rD 

X 

On voit que si Y est égal à —, l'équation 

A(HY) = HY -"" + ^^' 
^ ' ru 

est satisfaite pour H = a — j^ , Il en est de même si \ = -^ 

* dr r 



Àf 



ou Y^ = — . 



LES SURFACES DE NIVEAU SONT DES ELLIPSOÏDES 6 ) 

On a alors pour les termes du premier ordre 



dr \ 



r 



rt, i, c étant des constantes quelconques. — Nous verrons pro- 
chainement qu'on peut les supposer nulles, à condition de 
prendre le centre de gravité de la masse comme origine des 
coordonnées. 

Les surfaces de niveau sont des ellipsoïdes. — Considérons 
une surface de niveau primitivement sphcrique, lorsque la rota- 
tion est nulle. 

Pendant la rotation (fig. g), 
le point M est venu en M' sur la 
nouvelle surface de niveau. Le 
déplacement normal est 

î; = M]Vrcos(/-, MM'). 

Avant la déformation, on avait 
en. M Fig. 9. 

Après la déformation, la densité en M' est encore o.jLa dcn- ^J^ ^ 

AV 

site en M est 0^= — - — • , en appelant V le potentiel en M, et 




on a la relation 



AV = AV,-j; 



dr 
On a de plus 



rfAV„ rfAV, d\\ 



dr dV, dr 



et, en se servant des relations (4) et (5), 

AV = AV„ + -|i:(4D' + /-D"). 

PoiNCARÉ. — Figures d'équilibre, :> 



• # 
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Nous avons trailleiirs déjà écrit 

av=av.+Vaiiy. 

On a donc 

• ^ * 47:(4]yH-/-D") Zj 

Revenons maintenant aux équations déterminant les coeffi- 
cionlsll. 

Jj'augmenlatlon de volume entre la figure d'équilibre au repos, 
et la ligure d'équilibre dans le mouvement de rotation, est nulle. 
On a donc 



Îjrf3"=0, 



ou 





Ceci peut s'écrire 




\ Cp (/•) I YrfT= O. 



(Chaque terme s'annule séparément, sauf toutefois le premier, 
pour lequel Y est égal à Tunité. On a donc, en désignant le coef- 
cient de Y^ = i, par 11^, 

A(ig=o, 

mais IIq est uniquement fonction de R si bien que 

l /intégrale générale de cette équation est 

IL = A+— , 



,. 



A et B étant deux constantes. 



A PLA TISSEMEyr G; 

Mais remarquons que V = V^ -f- I II Y; à rinfini, V doit être 
nul comme V^. Donc A est nul. De môme B est nul ; autremeul 
le potentiel au centre serait infini, ce qui est absurde. 

I/équation (8) qui définit II se réduit dans ce cas à 






/"'(V„)-^ + '-?(V„) = o. 

C.ette équation définit o (VJ. 

Il nous reste encore à voir que pour les autres valeurs de Y, les 
seules valeurs possibles de II sont les valeurs II = o. Si nous 
admettons ce résultat, nous voyons que î^ se réduit à un seul 
terme, et que Ton a 

'\ AÏIY 

^ 4- 4iy + /jr -*^^ 

On a donc 



■> I •> •> 



Ceci montre que les surfaces de niveau sont des ellipsoïdes. 
En effet, on a 

{x' +u'+ -') = ('• + ^y = '" + 2/< + K' ; 

mais nous pouvons négliger ÎJ", et il reste 

^' -hf-hz'^ ,' + 2/- •'-' + ?/- a. ^ ,^ ^^.^ _ 

C'est Téquation d'un ellipsoïde de révolution autour de Taxe 0.3. 

Aplatissement. — L'aplatissement est le rapport de la 
différence entre le rayon équatorial et le rayon polaire, au 
rayon moyen, que nous pouvons supposer égal à r. 

Pour avoir ces rayons nous pouvons calculer Ç pour les valeurs 
^.r = y = Oy z = /•) et pour (y = c = o, .r :=: /•), on trouve les 
deux valeurs respectives 



-= — ar'y (/•) 



V "2 f 
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La diflereiice est 

et Taplatlssement est 

e = '-f-^ = ira (r) 

_ 3/< _ 3;: 

en désignant par Y la fonction spliérique 

■r^' + Z/" — a 3- 
, 

On a vu que 

3 XAIIY 



4^: /D'' + 4D'' 
La somme se réduit ici à un seul terme. L'équation (9) donne 



A (IIY) = (11 + Jl^) Y/' (V„) == (11 + 4I) Y ii^l+iP:^ 

Donc, 

4.,.^DeY = (lH-41)Y, 

On a d'ailleurs aussi 

A/-Y = o, 

]'iiisque Y est une fonction sphérique du second ordre. 
Ou a aussi 

-?- AIIY = [rD" + 4D') e/Y. 
41: 

L'équation devient donc 

Acz-DY = ^/■D"+ 4D') e/-Y. 
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On a démontré que si P est un polynôme sphérlque d'ordre /î, 



on a 



(UP)=p[u'' + ^(^ + i)-^]. 



Si nous posons pour /i = 2 

? = /•% 
On aura Téquation 



en développant les calculs on trouve, 

e"Y) + 2e'D' + -^ eD' + — p'D = o. 

C'est une équation linéaire du second ordre que Ton peut 
ramener a une équation du premier par le changement de 
variables. 

re' 
, re" re'^ . d 

V= — H — > 

e e e 



2 '/ 



= /V+^^ — *1- 



L^équation devient 

(12) D(/V + *^' + 5-/i)+2rD'(i+yi.) = o. 

C'est l'équation de Claivaut. 

Limites de r^, — Nous avons vu (p. 60) que Ton avait 

— 3<^-^<o. 

Dans l'équation (12), substituons ces valeurs extrêmes de rl>\ 
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et considérons les deux paraboles (fig. lo), où l'on a porté en 
abscisses r, et en ordonnées i-r' 



.■ + v,' + 5,-6(. 




Le point dont les coordonnées sont /y/ et /, ne peut se trouver 

qu'entre les deux paraboles (daps lu région hachurée sur la figure.) 

Je vais maintenant démontrer que r, est compris entre o et 3. 

r, ne peut être négatif. En ellet, s'il est négalif pour une 

)■„, on pourra trouver deux nombres positifsa et ^ tels que 



l'o 



— ?<-1< — ^, 
i et ^ pourront même satisfaire ii l'inégalité 
0<ï <a <5 <p. 
Je dis que c, ne peut dépasser — k. 



LIMITES DE r, 7 1 

En effet, si r, atteignait le valeur — a pour une valeur r infé- 
rieure à i\^ --—-serait négatif, /'—t^ serait négatif. 

Le point dont les coordonnées sont /r/' et vj = a ne se trou- 
verait pas alors dans la région hachurée. 

De même r, ne peut dépasser — ^3. En effet, si pour une valeur 

de /• inférieure à /• on avait /• < /• , -7^ serait positif, — r-^ serait 

^ dr cir 

positif aussi et le point représentatif serait dans une région 

du plan où le point (rr/, t^) ne peut se trouver. 

Donc si pour une valeur de j\^ t, est négatif, il est, pour toute 

valeur de /'inférieur à r^, compris entre deux nombres — a et — [4, 

on a donc 



re' 






— — 


-< 


— a; 


e 






c' 




a 




<- 




e 




/• 



lîn intégrant de i\ à /• </o, on a 



I,^.>L(^)-, 



^0 ^ 



e>A/— «. 



Or e croîtrait plus que lorsque /• tend vers zéro, ce qui 

est impossible pour /• suffisamment petit, car e serait très 
grand. Donc v) ne peut être négatif. 

Je dis maintenant que 7; est inférieur à 3. Le raisonnement est 
le même. 

En effet, supposons que pour r =: r^ on ait 7, = t^^ > 3, je dis 
que pour /• inférieur à i\^ on aurait encore •/; > 3. 

En effet, Tinégalité ne cesserait que si */; devenait inférieur 
à 3, ou si y; devenait négatif en passant par Tinfini. Cette dernière 
hypothèse est à écarter d'après ce que Ton a vu plus haut. 11 
reste donc à montrer que r, ne peut devenir égal a 3. En effet, 
pour Tj = 3, r/ serait positif, 1^ serait positif, ce qui est absurde, 
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le point (/t/, r) ne serait pas compris entre les deux courbes. 
On a donc toujours 

T, <3. 
On a d'ailleurs 

/•7/-|-7i^ — 'f\ — 6<o pourr, >;i. 
A fortiori on a 

/•r/ 4- /i^ — h — 6 — 5 (r^ — 3) < o, 

/V+(^.— 3)'<o 

et en multipliant par dr qui est négatif si r décroît de /^ à zéro, 
on a • 

dr , dr^ 

r "T TNr> ^» 



c'est-h-dire 



d"L \rd T-<o, 



T ' . " T ' 



/• 71—3 r, vio — 3 

Le second membre est une quantité finie ; le premier est une 
somme de deux termes, dont le premier seulement devient infini 
lorsque r tend vers o, l'inégalité est donc impossible. 

On a donc bien 



0<71<3 

Les limites peuvent-elles être atteintes? 
Si on a ■ - = o, l'équation se réduit à 

/V + 7^2+571 = 0, 

et admet la solution 

7i =0. 



C. Q. F. D. 
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Si— p-- = — 3, réquation se réduit à 

qui admet la solution r, = 3. 

Le premier cas correspond h un point intérieur de la masse 
supposée homogène; le second correspond aux points extérieurs 
à la masse, et en particulier au cas où toute la masse est concen- 
trée au centre. 

Je dis que dans ces cas particuliers, les seules intégrales 
admissibles sont les solutions vj = o et Tj = 3. 

En effet : 

Supposons - = o. L'équation se réduit à 

dr dr^ 

dont la solution est 

L/+4-L —^ = 0\ 



Pour r = o, on aurait 



7, + 5 



donc vi = — 5 ; or pour r quelconque, r^ est toujours positif : 
donc la constante d'intégration ne peut être que — oo , et l'inté- 
grale se réduit à 7\ = o. Il en est donc de même pour y, = 3, 
dans le cas de la masse concentrée au centre. A l'extérieur de la 
masse, on a l'équation 

/v/ + 7f — 7i — 6 = o, 
dont l'intégrale générale est 

La constante d'intégration ne peut être que o et on a y, = 3. 
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Ce sont les seuls cas oii Tj puisse atteindre ses limites. En effet, 
supposons d'abord que pour une valeur i\^ on ait r, = o, pour /• 
différent de zéro, r, est supérieur à zéro, donc 7/ [r^ = o, et en 
portant ces valeurs dans l'équation de Clairaut, il vient 

D' = o, 
mais comme 

/•D'+3D = 3o, 

on a D= p. La densité en ce point est donc égale à la den- 
sité moyenne, et comme par hypothèse la densité ne peut 
diminuer vers le centre, on a un noyau homogène. 

Supposons de même que pour /• = r^ on ait Tj = 3, de 
même que précédemment il faudra que y,' = o ; donc 



D 



p = o en ce point, et pour r > r^ on a encore p = o ; p ne peut 
augmenter puisqu'il y a équilibre, et par suite r, = 3 pour /• > 7^. 

11 reste à voir ce qui se passe pour r < r^. Je dis que Ton a 
encore Tt == 3. S'il en était autrement, r^ serait plus petit que 3, 
r^ (t, — 3) serait une quantité qui diminuerait avec r, sa dérivée 
serait positive. On aurait donc 

/•r/4-5(T. — 3)>o. 
Or, l'équation de Clairaut peut s'écrire 

rV + 5 (t. - 3) + (r, - 3)^ + [ ^ + e] f I + vi] =: o. 

Les deux deriiiers termes étant essentiellement positifs, il en 
résulte que l'on devra avoir 

7V + 5(t.— 3)<o 

ce qui est contraire à ce que l'on a établi plus haut. 
Pour /• < /-y, on a donc encore t^ = 3, et par suite 



r^+el=o. 



FORME DES COURBES ISTÉGRALES ^^ 

Donc p = o pour /• < 7*^. 

La niasse de la planète est donc entièrement concentrée au 
•centre. 

Forme des courbes intégrales, — Etudions maintenant la 
iorme des courbes représentées par Téquation 

2/D' 

'V+V + 5t,H — — (i 4-r,) =o. 

Traçons les axes de coordonnées Or, et Or, 
Il suffit évidemment d'étudier les courbes comprises dans les 
régions où Ton a 

o<Y, <3. 

Traçons la droite D d'équation r» = 3. 

Par un point du plan passe en général une courbe et une seule, 
il moins que t/ ne soit pas une fonction holomorphe. 
Or, on a 

2/D' 



y/ 



r 



7/ n'est pas holomorphe si r == o ; la droite /• = o est une courbe 
intégrale de l'équation ; il y a sur cette droite deux points remar- 
•cjuables Tj = o et 7| = — 5, pour lesquels y/ est indéterminé 

il condition que — p-— tende en effet vers zéro avec r. Or, on a 



D —^ VD 



Le second membre est essentiellement négatif, la densité d'une 
«ouche étant inférieure à la densité moyenne de la masse inté- 
rieure à cette couche, puisque par hypothèse il y a équilibre. 

Lorsque r tend vers zéro, p tend vers une limite p^ et l'on a 

?<ï^<?o. 
■^ i< -^ — \<o; 
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— • — I tend vers zéro avec r, donc -Jr i tend au 



zéro et par suite —r^ tend vers zéro avei 
Aux deux points 



les deux coui'bcs intégrales ne sont pas déûnîes d'une façoi 
unique. 

Faisons mtiintennnt r suOisamment grand, alors 



/■D' 






L'équation devient 

Klle s'intègre iiniuédiatenient, et son intégrale générale esl 



= C". 



Les courbes représentatives sont toutes asymptotes ii la droilo 

r, =3, mais elles ne sont pas nécessairement asymptotes à la droite 

/■D' 
;■= o, car on a -j^ :^ — ;i pour 

j- = o : le raisonnement ne s"ap- 
pliijuc plus. 

Considérons le rectangle OABC 
(fig. Il) formé par les droites 

étant une valeur quelconque. Par 
un point P, pris sur le côté OA, 
passe une courbe intégrale et une 
seule. Pour cette valeur 




■IV 



r>o 



la portion de la courbe qui pénètre dans le rectangle est doi 



FORME DES COURBES ISTÉGRALES 



i y 



P' 




dirigée vers le côté AB, quel que soit le point P. Une courbe 

intégrale ne peut donc sortir du rectangle par le côté OA, car si 

elle ressortait en un point P', elle aurait une forme analogue à 

la courbe PMP' (fig. 12), ce qui est impos- 

dr 
sible puisque -^ — ne peut être nul dans le 

rectangle. 

Donc cette courbe ne peut sortir du rec- ^ 

tangle que par un point M de la droite AB. — ô" 
De même, si l'on considère la courbe inté- yw. 12. 

grale passant par un point Q du segment 
GB, la tangente en ce point a comme coefficient angulaire 

dr r 



rfri 



4K-^] 



quantité qui est négative puisque Ton a 

2/D' 



D 



— 6. 



La branche de courbe qui pénètre dans le rectangle est donc 
dirigée vers l'arc AB, le même raisonnement que précédemment 
montre que cette courbe ne peut sortir du rectangle que par un 
point N de la droite AB. Evidemment, d'ailleurs, le point M est 
à gauche du point N puisque deux courbes intégrales ne peuvent 
pas se couper. 

Les points M ont donc une limite M' lorsque le point P se rap- 
proche de G, de même les points N ont une limite N' lorsque le 
point Q se rapproche de G. 

Le lieu des points M' sera une courbe complètement comprise 
entre les droites t^z=^o^t\ =3. Il en est de même pour les lieux 
des points N'. 

Est-il possible que M' et N' coïncident ? Supposons, ce qui 
sera démontré un. peu plus loin que cela arrive. Alors il y 
aura une courbe intégrale et une seule qui passe par le point O, 
Gette courbe sera évidemment la courbe intégrale convenable, 
puisqu'elle restera constamment comprise entre les droites y^ = o, 
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7j =1= 3. lien résulte que la fonction r, (/•) est parfaitement déter- 

mince, est donc parfaitement déterminé ; e est alors, en 

e 

désignant par E la base des logarithmes népériens 







La constante e^ sera connue si on se sert de Tune des équations 
précédentes, par exemple 



e/*' 



■>=f ("+4^)- 



Une seule courbe est acceptable. — Démontrons maintenant 
qu'il y a une seule courbe intégrale acceptable. Nous allons, pour 
cela, démontrer qu'il y en a une et une seule qui passe par le point 
/•izzizo, v) = o, et que toutes les autres passent par le point /• = o, 
7, = — 5. 

Démontrons d'abord qu'il y a au moins une courbe passant 
par le point vj = o, et une passant par y^ = — 5. 

Nous avons déjà vu quel était le signe de rr/ pour diverses 
valeurs de vj ; r/i' est négatif si T^ est inférieur à — 5, ou supé- 
rieur a + 3 et positif si vj est conipris entre o et — 2. 

Dans les autres intervalles, les deux signes sont possibles. 

Pour 7) = — 5, /"Yj' peut être nul ou négatif, mais non positif. 

Pour vj = — 2, /t/ ne peut être que positif ou nul. 

/•D' 
Enfin, nous avons vu que — — — est une fonction tendant vers o 

avec r, de façon, qu'étant donné un nombre positif a voisin de 
zéro, on peut déterminer un nombre r\ tel que l'inégalité 

entraine 

o>— — > — a. 

L'équation de Clairaut peut s'écrire 
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Ji et vêlant des nombres positifs, tendant Tun vers zéro, raiitre 
vers 5 lorsque a tend vers zéro. 

Ceci va nous donner des indications sur le signe de rr/, lorsque r 
est inférieur à i\. On a le tableau suivant: 



•^1 



/•// 



00 



5 — Y o j3. -f-^ 



•> 



1 + 



Il y a deux intervalles pour lequel le signe de /t/ est douteux, 
mais ces intervalles sont moindres que ceux précédemment 
déterminés. 

Ceci posé, je vais démontrer que pour /'</',, r, est nécessaire- 
ment inférieur a ^. Supposons que pour une valeur /'2</'i on ait 
T, = Tj^> |3. Je dis que si r < i\ on aura encore fi > p. En effet, 
l'inégalité ne pourrait cesser que si t'== ^, mais*alors dans l'in- 
tervalle, r/sera positif, ce qui est contraire à la règle des signes. 

Or l'équation peut s'écrire : 

,V+(yi-?)(r. + T)+2(-^ + a) (. +r.) = o; 



par hypothèse 



/•D' 



D 



-> — a, 



on a donc 



'V+(-0 — ?)(^. + ï)<o 



A fortiori on a 



'V + (^ — ?)'<o 



. I 



et, si on Intègre de /•, à r<.r^, 



+ — <0' 



L 


I 


— 




I 




>- 


-L 


I 
'1 






I 




f^ 




? 


■'.. 




? 


L 


I 


' + 




I 




<L 


I 






I 






U 




^ 


-'n 




? 
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Or 7j — P serait positif : on aurait a fortiori 

L — <L h 



Le second membre est fixe, le premier peut grandir au delà de 
toute limite, Thypothëse faite est donc absurde. Par suite 

Tj a pour limite zéro lorsque /• tend vers zéro. 
Je dis maintenant que l'équation de Clairaut admet au moins 
une intégrale tendant vers — 5. 

En effet, étant donnée une valeur t,j comprise entre — y et — 5 

— 5<7)^<— y, 

considérons une courbe intégrale passant par le point (rj,r,J, je 
dis que si r^ < z*^, r,^ <-y. Autrement y/ serait négatif pour Tj > — • Y, 
ce qui est contraire a la règle des signes. 
On a donc pour /• suffisamment petit 

— 5<r,<— y; 

mais y pourra être choisi aussi voisin de 5 que Ton voudra ; il en 
résulte quer, tend vers — 5. 

Soient Yj^ et y^j les deux intégrales dont nous venons de démon- 
trer Texistence. Il leur correspond des intégrales e^ et e^ de 
Téquation linéaire du second ordre 

\e'D + 2e'D'4- Ad' + — e'D = o. 
L'intégrale générale de cette équation sera 

A et B étant deux constantes. L'intégrale générale de Téqua- 
tion de Clairault sera 

re^ kre^ -f- ^re[ 



'f\ 



Ae, + Be, 



Ar,A+ By,, 

A-^ +B 



RELATION ENTRE V APLATISSEMENT, LA PESANTEUR, ETC. 8l 

Or, quand r tend vers zéro, r,^ tend vers — 5, Tj^ tend vers zéro. 
Il en résulte que tjj — r.^ pourra être Inférieur h — 4- On aura 
donc 

re/ /•(?' 



'- ^< — 4. 



c c 






^1 ^0 '• 



• Vf : 



L i 



^i<^o''" ; r, " 



«■ 



ç 

c'est -à-dîre que la limite de —^ est zéro quand /• tend vers 

zéro. Donc quand r tend vers zéro, la limite de r, est la même que « , 
celle de r.^ tant que B est diflFérent de zéro. Si B = o la limite 
est rjjj. Il en résulte qu'il y a une intégrale et une seule de l'équa- 
tion de Clairaut tendant vers zéro avec r : c'est celle qui corres- 
pond a B = o. 

Relation entre V aplatissement^ la pesanteur et la force 
centrifuge à l'équateur, — A l'extérieur de la planète, comme 
cela a déjà été remarqué, on ne peut parler de surfaces d'égale 
densité, mais on peut encore parler de surfaces de niveau. 

Pour ces surfaces, on a toujours l'équation 

AV=AV, + A(IIY)=o; 

mais AVq = — 4^p est nul à l'extérieur de la masse fluide. On a 
donc 

A(IIY) =0 
A(HY) =y| if +-^ ^,_ /?(/i+i) jjl ^ 

Faisons n= 2, on a 

IF+— IF — -4-H = o. 

L'intégrale générale de cette équation est 

H=A/-2+B/-3; 

PoiNCARÉ. — Figures d'équilibre. G 
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mais à rinfini II doit être nul, il faut donc que A = o. On a donc 

II = Br-3 ; 

si M représente la masse fluide, Y le volume enfermé par la sur- 
face de niveau considérée, on a 

en posant L = — - — . 

47w 



L'équation 



ei 



devient 



D'où 



r-C _ 9 r B o)',-^ 1 

9r B coVI 

47: L C/-- ^ 6C J • 



4- L C/- ^ 6C J 



On a donc 






re' +^e — -f-. 




i5 ,3 
87.C "■'• 




iSco^ 



StiD • 



Comme on a posé 



on a finalement 



re'=er^ 



i5 W-' 



<-^.+ 2)=-8^.-^. 
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Cette relation établie pour Textcrleur de la niasse, est encore 
valable à la surface. On peut lui donner une autre forme. 

La force centrifuge à Téquateur est égale à wV ; la pesanteur 

à la surface est -5- tc/'D. Si on désigne par tp, le rapport de la 

première h la seconde de ces quantités, on a 

Or ■t\ est compris entre o et 3 ; il en résulte la double inégalité 



4^2 

Tj == 3 seulement si la planète est homogène. L'aplatissement 
ne peut atteindre la limite -^que si la planète est homogène. 

Détermination complète des coefâcients. — Nous allons 
maintenant démontrer que les équations 

AHY:=-P^+/P'^ IIY 

/•D 

qui déterminent les coefficients II du développement 



v=v,+2hy, 



n'admettent que la solution H = o. 

En effet, soit n Tordre de la fonction sphérique, posons 

H = Z/«D. 
L'équation devient 

A(Z/'»DY) = ^^ Z/»Y ; 

/•"Y est un polynôme sphérlque ; en appliquant la formule 

AfPU) = 



(pu)=ru'' + -^^^y^ u^ 1p, 



# » 
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on obtient finalement 



/• ' ' /• ' 



ce qui, développé, s'écrit : 

Z"D + aZ' (d' + -Îi±-L D W .^^i=^ ZD' = o . 

Si l'on fait le même changement de variables 

/•Z' 

-z- = ^' 

on a 

rt + £2 ^ (2/1 -h 1 )£ H j^ (e + /i — I ) = o . 

Cette équation se discute de la même manière que l'équation de 
Clairaut ; — — — est compris entre les limites o et — 6. 
On est donc amené à considérer les deux équations : 

re'H- £^ +(2/i-|- i) £ = 0, 

/V + £'^ + (2/1 — 5)£ — 6 (/i — i) =o. 

On construit les courbes en portant £ en abscisses et /V en 
ordonnées (fig. i3) ; on obtient deux paraboles, et comme pré- 
cédemment on démontre que £ ne peut être compris qu'entre 
o et 3. 

A l'extérieur de la planète on a A (HY) =o, ce qui peut s'écrire 

r r 

L'intégrale générale de cette équation est 

II devant s'annuler à l'infini, on a a = o, et par suite 

C 
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Vlaiss doit être compris entre o et 3 ; donc b doit ètro nul. 




Fîg. i3. 

H^o est donc la seule solution possible, à l'extérieur de la 
masse. Il en est donc de même à l'intérieur. 

Ce résultat ne subsiste pas si n = i ou i, mais nous avons 
déjà examine ces cas. 

Détermina f ion précise de l'aplatissement. — Nous avons 
trouvé 

i(UY)=4ii,a'i;, 



Nous avons donc Téquation 

i(HY) = 



- rap'i'Y, , 
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et, en se rappelant la formule 

A(HY) =y[h"+-^- <^±ilH], 

dans laquelle /i == 2, on arrive pour déterminer II à l'équation 
linéaire avec second membre : 

2ir 6H 4 , 

L'intégrale générale de l'équation sans second membre 
étant 

où a et p sont des constantes d'intégration, on déterminera l'in- 
tégrale de l'équation avec second membre, en regardant a et ^ 
comme des fonctions de /•. La méthode de variation des cons- 
tantes donne 

(i) aV2+pV-« = o, 

par suite 

H'=z=2ar — 3?/-'^ 

H'' = 2aV— 3pV-'* + 2a+ i2p/- 

et, portant ces valeurs dans l'équation primitive, on a 
(2) 2aV — 3|3V-^ = -|- TrepV, 

Les équations (i) et (2) déterminent a' et ^' 



— 5 



a 



/ 



47t , 



i5 ^^' 



P'=_4j.pv^, 



i5 

On a donc 



ec^'dr— ^ ^ I e^'r'^dr, 
' i5 f ' 
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Ce^le fonction comprend deux constantes arbitraires a et h.. 
Comment faut-il les choisir ? Soit r, le rayon de la planète ; pour r 
supérieur h r^, p = o, f ' = o, par suite les intégrales se réduisent 
à des constantes 




et 




Pour /'=oo, Il doit s'annuler. Donc le coefficient de /• es'l 
nul ; on a donc a = Tj. 

De même pour r = o^II doit être fini; donc le coefficient de 
/• '^ doit être nul, si /• = o, on a donc Z> = o. 

II est ainsi complètement déterminé. 

Portons cette valeur de H dans l'équation 



ei 



•^D=--^ [11+ 



0)2/2 -, 



47^^ ^ 



11 vient 




5ï> / ^?'*-w / '^''"^''-^-w4 




Que devient cette expression à la surface pour /• = r,?D devient 
D, la densité moyenne de la planète, et on a 



e. 



>~ 5D,/' / ' ' 8- D, 



(3) 



6,=^ 




9 



cp étant comme on Ta déjà indiqué le rapport de la force cen- 
trifuge à l'équateur, à la pesanteur. 



88 MASSE FLUIDE HÉTÉROGÈNE, PROBLÈME DE CLAIRAUT 

Moments d'inertie de l'ellipsoïde, — Calculons Tîntégrale 
/ Y^rfo-, Y étant la fonction 



r' 



On a 

X = /• y I — i*} cos y, 

y = /• y/i — [JL^ sin cp, 



Z = /-[A, 



On a donc a calculer 




ce qui donne 



(4) 



+1 /»+l 

I I 3jJL^ I d^kd^f = 27Z j I I 6|X*+9[JL* |rf|X 

— 1 



Calculons les moments d'inertie A, B, G de la masse 



On a 



Or on a 



A V = — 4^p , 

-4TrA=jAV(y^ + ^VT, 

— 4-B =j'^y{z^-i.x^)d1, 

— 4itC = fà\{x^ + y^)di: 

— 27t(A + B + C) = j'^\{x^-\-y-+z'')d'z. 

AV = AV, + A(HY) 
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et de plus 

rA(HY)/VT = o; 



en effet, on a d'z = r'^drd\kd^ = drdfj. 
Il vient 



— ' .^ffl- ;.V,..Yrfcr = I f{r)dv 



h 



ce qui est nul, car le second facteur est nul. On a donc. 

27t( A H- B + C) = rAV,/-VT 

et comme AV, ne dépend que de r, on a 

27t(A + B + C) = 47: C^yydr. 
Or 

' ^^0 — ' dr' +^' dr ~ dr L dr J' 

mais on a 

dN, 47^/D 



donc 




'"-JF ^^'"^'^''' 



Si on intègre par parties et que l'on remarque que A = B, on 
voit que Ton a 



m P 
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En négligeant des intiaiment petits on a A == C ; et par suite 



= |L[,o,_ 




Calculons maintenant C — A : on a 

4-(C — A) =.Ç\V{z^—x'-)d-, 

4r.{C — h)=flY(z*—y^)dr. 

Donc 

8-(C— A) =—f^Yr^Ydr: 

2 I 2 *> 



en désignant par Y la fonction sphérique 
On a, comme on a vu précédemment : 



r' 



Stt 



[C — k)-^—f\\^r^^ch — fl[Wi)r^\d':. 



La première intégrale est nulle, car c'est le produit de deux 
intégrales : 



flY,r'drJ 



YrfT. 



dont la seconde est nulle. L'intégrale 



A(HY)/nWcr= / ^(^^) , ;.2YVt 



/ A(HY) ,^ / ^.,, 



et comme on a 



(4) / V^i. = ^ 



D 



î 
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il vient 



c^ fr^ A\ i6k I A IIY) , , 

87:(C — A)= ^ I ■ ' ^ r'dr 



167c 14 / • ï 



91 




64-- 



Donc on a 



C— A = ^ / e^'r'dr 




i 



quantité positive. 

On a donc en se reportant à la formule (3) page 87 

. ? _ 9 (C-A) 

On aurait pu faire une hypothèse différente de celle de Clai- 
raut, par exemple, supposer un noyau solide non homogène, 
recouvert d'une masse fluide; la relation établie plus haut sub- 
sisterait encore. 



Comparaison de la théorie et des observations. — Posons 

C — A 

— = J, on a 



2 ~ /JD, L * 




.r, 

2 / r'Ddr 






» 



92 
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e. est mesuré par la géodésie, on trouve — rr—p î '-^ <^st mesuré 
* * ^ 290,5 * 

par la physique, on trouve — - ; J est mesuré par la préces- 

sion des équlnoxes, on trouve 0,0032^53 = ^ ^ on a donc 

00 J j o I 

Tégalité 



1 ^'i 



Revenons maintenant h l'équation de Clairaut 

/•7)'D + ('/)'+ 5yi)D + 2rD'(i +'/i) = o ; 



on a 






_ 5;*D r 



/V /-(ir!--/!) D^ 



10 D D J 



On peut donc écrire l'équation de Clairaut 



^[^VTT.o]=-^[.+,-^-i] 



J'intègre de o à /•, et j'ai 



4/1+V1 L 2 10 J 




1 + 



2 1 
2 10 J 



On peut appliquer le théorème de la moyenne, puisque 5/'^ D 
est une fonction positive ; Ç désignant une valeur comprise 
entre o pt 3, et correspondant à une valeur de /* comprise entre 
o et f\, j'ai 






Dt//\ 



COMPARAISON DE LA THÉORIE ET DES OBSERVATIONS y 5 

Voyons les variations de 

' + — — 



^ "L^K; 



+ 5 

lorsque 5 varie de o ii 3 on trouve : 

l = o K=i 

? = -TT- K = I ,00075 maximum 

i = — K = 1 , 000a 

2 

\ z= 0,53 K = I 

b = -^i = 0,544 K, = 0,99954 

i = I K = 0,989 

; = 2 K = o,9 

Ç = 3 K = 0,8. 

On voit que K varie très peu et reste compris entre 1,00075, 
et o, 99954. K est donc très voisin de i ; et si dans la formule ou 
remplace v), par sa valeur 0,544? on trouve 

= 7n7 \/* + ''il 



f\'D^ 5K 

= o>497-j^>^>49663. 

Précédemment on a trouvé pour la même valeur 

0,49145. 

Comment peut-on expliquer cette différence? 
On peut admettre qu'il y a des erreurs sur les quantités 
qui entrent dans le calcul : c'est-à-dire des erreurs sur J et e^. 

La première équation donnerait alors comme valeur exacte : 



'■/'-'' 



«P 8J 



,1D, - = «49'45 - -5- Se, + -^ 

= 0,49145 — 3o5,3i 8(?^ + 155,26 oj. 



* # 
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La seconde équation donnerait : 

_ cLK(r'\ydr / ^ ^^ 

2 J — (7)^+2)rfe 

J =0,49700 120,10 0^^. 

Supposons 5J = o; pour que ces formules soient comparables, 
il faudrait que 5^^ soit égal à — 0,000027 et la valeur de e^ serait 

alors TT-FT^' 

295,85 

Si on suppose que Terreur est commise sur J, on a 

oJ= 0,000034. 

J'aurai alors pour valeur -rr . 

'^ 002,18 

Cela est-il admissible ? Remarquons d'abord que le J qui 

(] ^ 

figure dans cette formule n'est qu'approximativement — — — . 

Cela revient à remplacer J par 

J 

'^ 

2 

L'erreur commise est 3J=-^ J^ : on trouve 

o 

8J = 0,000006. 

Ce n'est pas de là que provient l'erreur. 

D'autre part les mesures de nutation sont assez précises. Si 
on désigne par N la constante de la nutation, on a 

N = — ^J [5,36542]; 

I + [A 

en désignant par [jl le rapport — de la masse lunaire, à la masse 



COMPARAISON DE LA THEORIE ET DES OBSERVATIONS 9$ 

terrestre ; le nombre entre parenthèse étant non le coefficient 
mais le logarithme de ce nombre. 

Soit p la précession des équinoxes : c'est la somme de deux 
termes 

P^p'+p", 

celui dû à l'attraction de la lune 

p' = 866i35 — £ — J = 34,38 ; 

I — }— [A. 

celui dû à l'attraction du soleil 

p'^ = 48^1, o5 J= 15,95. 
Posons 



on aura 





V- 




s : 




I 


+ i^ 




5N 


0£ 


1 


ûj 
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e 
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J ' 




-4 ' 
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■^~1 




+/'• 


5J 
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+p 


J 


't 





P est mieux connu que N, on peut donc supposer 5/; = o ; et 

gj 

résolvant les deux équations où Ton regarde— r- comme égal 
o,oooo34 X 3o5,3i =o,oio3, on trouve : 






s 

8N 



= OjOID 



= o,oo5. 



N 

Cela revient à prendre pour N la valeur 9'' 17 au lieu de9''2i 

et pour [ji la valeur -^ — ^ au lieu de-^. Ces chiffres paraissent 

difficiles à admettre. 

On peut alors supposer que l'erreur provient de la somme des 
erreurs commises sur chacun des termes. 
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Une autre explication serait que l'hypothèse de Clairaut n'est 
approchée qu*aux termes du second ordre près. M. Callandreau, 
qui a fait le calcul, a montré que considérer les termes du second 
ordre reviendrait à creuser l'ellipsoïde de 9 mètres seulement. 
C'est négligeable. 

On pourrait alors se demander jusqu'à quel point l'hypothèse 
de Clairaut serait exacte : par exemple, on peut supposer la terre 
complètement solide sauf la mer. Au moment de la solidifica- 
tion, les couches obéissaient à la loi de Clairaut, mais depuis ce 
moment, la vitesse de rotation aura pu varier (par exemple, les 
marées diminuent cette vitesse de rotation) ; et de plus l'aplatis- 
sement des ellipsoïdes internes aura pu se modifier par suite des 
tassements qui ont produit les continents et les montagnes. 

Enfin on pourrait supposer qu'une portion intérieure de la 
terre est restée fluide : les deux parties fluide et solide, agiraient 
alors séparément dans le phénomène de la précession : mais 
l'hypothèse de la fluidité intérieure de la terre n'est guère vrai- 
semblable. 

Jusqu'à présent nous avons considéré les planètes comme 
fluides et nous avons montré quelle devait être la distribution 
des densités à l'intérieur. 

Supposons maintenant que l'on ait une masse solide, de forme 
quelconque, où la distribution des densités est quelconque, mais 
entièrement recouverte d'une masse fluide de faible épaisseur, 
soumise à la pesanteur et à la force centrifuge due à la rotation. 

La surface extérieure sera une surface de niveau, sur laquelle 
on aura 

mais à l'intérieur les surfaces équipotentielles ne sont pas néces- 
sairement les surfaces d'égale densité. 






CHAPITRE V 
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Théorème de Stokes. — Supposons que Von ail déterminé la 
forme de la surface extérieure 5, la masse M et la rotation w, 
on pourra déterminer en tout point extérieur la valeur du poten^ 
tiel V. 

Il y a une infinité de lois de distribution de la matière dans le 
noyau, qui peuvent donner une même surface extérieure, mais 
quelle que soit cette distribution, le potentiel V, h l'extérieur, est 
complètement déterminé par les éléments S, M et o). 

En effet la fonction V satisfait à l'extérieur de la masse à l'équa- 
tion AV = o. A l'infini V = o ; sur la surface on a 



V+-î^(x^+y^) = U = K, 



et enfin l'égalité 



acr= — 47ïM. 



dn 



Le problème de la détermination de V est complètement déter- 
miné. En effet, soit V une autre solution. Sur la surface, on a 

V— V'=U — U' = K — K'=constante. 

On a d'ailleurs A (V— V) =o. 

PoiKCARÉ. — Figures d'équilibre. 7 
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Donc V — V est le potentiel d'une charge électrique répartie 
sur la surface, puisque c'est une fonction constante sur cette 
surface. 

La charge a pour expression 

d/e 

l'intégrale étant étendue à la surface extérieure ; mais elle est 
nulle puisque l'on a par hypothèse 




d/e 



(17= — 4'ï^M, 



à/i 



V — V^. le potentiel d'une charge nulle, est donc nul. C'est ce que 
Ton voulait démontrer. 

L'expérience du pendule qui peut faire connaître les valeurs 
de g à diverses altitudes ne peut donc rien nous apprendre sur la 
distribution des masses intérieures, que ne nous fassent connaître 
les éléments considérés dans le théorème de Stokes. De même, 
nous ne pouvons rien apprendre de plus par les perturbations 
produites sur le mouvement des autres astres. 

Développement de ^^, en série de polynômes. — Soient 

O l'origine des coordonnées (fig. 14), M' un point attirant de 
coordonnées x\ y\ z' , de densité p', M {x^ y, z) un point attiré. 
En prenant les coordonnées polaires on a : 

X = /-y I — [j? cos (p, 
y ==; /-y/ 1 — [jL^ sin ç, 

z = /-[A, 




[ x' = r'\/i — ^'^ coscp', 
I y^==r'\/i — [ji'^ sin y', 



ig. 14. \ ^ — / [i. , 



^ :a 
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MM'^ = r^ — 2 r/' cos y + /'^ 



= a;'^y' + z' — 2{xx' -^yy'-^ zz') -\-x'' -h y'' + z 

I I 



«/2 



MM' 


. / 2 /-'cos Y , /-'^ 




'•'— ^'+y'+A 


■ 


,.'■!— ^c'î+yî+s/î, 


osY = 


/• /•' /• /•' r r' 



Développons -TjTjrtrp- suivant les puissances croissantes de ; 

le coefficient de — -^ sera un polynôme U^ [x, y, z ; x\ y\ z') . 
On doit avoir les deux relations 

la première étant obtenue en considérant x, y, z comme varia- 
bles, x'y'z\ comme fixes; c'est Tinverse dans la seconde. 



^ V MM' / ~" 
\MM7 



52 ,^4^ 52 ..JL^ 52 ï 



MM' MM' MM' 

— = G, 



do;^ ^n^ ^z 





d^ 


I 
MM' 




^f 


1 


d^ 


I 
MM' 



^, I 



MM' MM' MM' 

H r-75 — =0 



d.r'2 ' ô?/'^ ' d^'^ 

On sait que l'on peut écrire 

I ^7-' . ^ / /•' ^ " 



\/ 



1—2 --COSYH j 



^=' + p.-+--- + p«^T 



où P„ est un polynôme en cos y de la forme 

AgCos^Y + AgCos^'^y 4" ' 
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si on y remplace cos y par sa valeur, on aura 

Le polynôme Q sera un polynôme homogène et de degré n 
par rapport h chacun des systèmes de variables [x^ y^ z^ /*), 
(x' ^ y' f z' y r'). Il ne contiendra d'ailleurs /• et /•' qu'à des puissances 
paires, on pourra donc remplacer r^ et /•'^ par leurs valeurs 

;j;2 _|_ y2 _|_ ^2 ^t ^/2 J^ ytl ^ Jl . 

on aura donc : 



,./2 Za 






r' 



cosy+ 






MM 

Q' représente un polynôme homogène et de degré n par rapport 
à chacun des systèmes de variables x, y, z; x\ y\ s', chacun des 
deux ensembles y entrant de la même façon que Tautre, 
On a d'ailleurs 



Donc on a 



A'-2-=o; 



et, comme r ne dépend pas de x\ ?/', z\ on a : 

A'(Q')=a. 

Q' est donc un polynôme sphérique d'ordre n par rapport aux 
variables x' ^ y\ z' et par suite par rapport aux variables [x^ y, £), 
Nous savons qu'il y a 2/?+ i polynômes sphériques d'ordre/;, 
de forme 
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OÙ q prend les valeurs 

— P, — 7^+1, — i,o, I, 2, . . . p — i,p. 

On peut donc écrire 

formule où P est un polynôme fondamental en .r, y, r, Q un 
polynôme en {x' , y', z') ; Q' est d'ailleurs fonction de 



cos 



Y = a[x' + y/i — 1^ V^ — I^'' cos ('f — cp'j 



Q' ne dépend donc que de la différence (o — '^') : il en résulte 
que le coefficient de P^^ne peut être que 



;^' 



et on pourra écrire : 






MM' 



A^^ est un coefficient numérique que Ton peut déterminer faci- 
lement. 

Application. — Supposons que /• soit suffisamment grand par 
rapport au rayon de la planète. On pourra développer le poten- 
tiel suivant les puissances croissantes de et on aura 

'^ r 

B.„ P. 



V \^ , P9 ^ PQ 



D'après ce que Ton a démontré plus haut, la connaissance de 
S,Met(o suffit pour déterminer V, donc ces éléments permet- 
lent de déterminer B^, ^ ; or, on a 



\^ ^PQ^ PQ J ^/p/ 



— ' MM' — /. ..2B-+-1 f P^i>.-?"''» 
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il vient donc 



= A^ / p'P', 



B.,=A^ I pn..,^^ 



P';,, q est l'un des polynômes sphériques fondamentaux d'ordre p. 
On doit donc savoir calculer 



P'PV^. 



Considérons d'abord le polynôme d'ordre o : P^, ^ = i ; on a 



p'd:T' = M. 



Les polynômes d'ordre i sont x\ y\ z' : on connaît donc la valeur 
des intégrales 



1 pVrfr, 1 py^T, 



p z d'z. 



Ce sont, h un facteur près, les coordonnées du centre de gra- 
vité de la masse. Supposons que la surface extérieure possède 
un centre de symétrie, et prenons l'origine des coordonnées en 
ce centre ; les intégrales sont nulles puisque ce sont les coeffi- 
cients des termes en —3-, -^, -^, dans le développement du 

potentiel et que les surfaces équipotentielles doivent admettre 
l'origine comme centre de symétrie. Le centre de gravité de la 
planète coïncide avec le centre de symétrie. 
Les polynômes de second ordre sont 

^'y\ y'z\ ^'-^7 
Si / x' y' ^'d'z est nulle, ainsi que les deux autres intégrales 
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analogues, les axes de coordonnées sont les axes principaux 
d'inertie du système. Si / (o;'^ — y'^) prfw est connu, on connaî- 
tra la différence entre les moments principaux d'inertie. 

On a trouvé une relation entre e^, j\y D^ et w ; les deux premiers 
éléments sont donnés par la connaissance de S ; D^ est connu 
puisque Ton connaît M et S ; cp est connu puisque Ton connaît co ; 
C — A est aussi connu. On doit en conclure que la relation déduite 
de l'équation de Clairaut est vraie, quand même la distribution 
intérieure de la matière ne serait pas conforme à l'hypothèse 
faite pour l'établir. 

Relation entre l'aplatissement, la pesanteur et la force 
centrifuge à Féquatear. — Proposons-nous de calculer les pre- 
miers termes du développement de ^...^ • 

On a d'après la formule du binôme 

I I r I / 2/-' , /•'* \ 

+ _^_cosY— ^) + • • -J' 
négligeant les termes de l'ordre de — j-î ï^ous avons 



I / 



. rcosY , î r ' . ^ /> 2 1 

MM' - ,-+-7^ + 7^L— T + T'-"*=°^*^J' 

IM' /• ^ /•» L ^^ 2 I 



MM 

3 



2/ 



- {.Tx' + yy' + zz'Y 



= -7 + TT {^^'+yy' + -') 

-[(a:^ + y^-^z')ix'' + y''+z'')-i(a:x'-hyy' + zz'y]. 



2/ 

Le coefficient du terme- — ^ peut s'écrire 



3^ 

— 6 {xx' y y' -]-xx' zz' + y y' zz') , 



-4 [x'-f) [x'^-y'^) -- [x' + f—iz^) [x'-^ + y--^z'^) 
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On a alors 




+Y,tA^''''' 



+ 3(.r^ - 



V— I P'^' — 

MM'~ 



^î^ / (^--y-)p'./T 



I (X^ + y^ - 2C ^) 

4 '•' 



/ (y^ + y^_25'0pW-r 



+-iS^ Av?'''-' 



SI M est la masse, a, j3, y les coordonnées du centre de gravité, 
A, B, C les axes principaux d'inertie, on a 

V= -^ + "^ + ^f+-^ +-4^[3(-'-3.^) (B-A) 



+(x^ +y^ — 2^^) (_ A — B + 2C) -f-; 



-js^ A'^>''^^-i 



Supposons que la surface soit de révolution, que l'origine des 
coordonnées soit placée au centre, que les axes de coordonnées 
soient les axes principaux d'inertie : alors 

B = A 

et 

Calculons maintenant la valeur que Ton déduit pour V de la 
forme extérieure de la surface, que nous supposons un ellipsoïde 
d'aplatissement c,. Ce potentiel est de la forme 

V=-^+IIYH-SQZ, 
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OÙ Y est la fonction 

et Z est une fonction sphérique différente de Y. 

H et Q doivent être d'une forme particulière, puisque lu poten- 
tiel doit s'annuler a l'infini, nous avons déjà vu que Ton avait 



^.3 ' 



Hy et q^ sont des constantes; on a donc 

M , HJ y y„Z 
r ~^ ,.3 -1-^,.» + 

A la surface libre on a 



2 9 •>V 



11 vient donc 






Posons 

avait 



^ 



On a vu (p, 6b) que l'on avai 



1 

' 1 



/•.Y ' 



de plus on a 
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Tous les termes étant petits, sauf le premier, on peut y négli- 
ger e\ et on a 

M M M / e.Y 



r / eM 

r 



■■'(■+^) 



"(-^) 



L'équation U = C*® se réduit donc à 

M Me. Y «oV? «oMY , H„Y Sy.Z _ 

r. 3r, "^ 3 "^ 6 "^ /? "^/•ï+' ~ 

Comme par hypothèse nous avons un ellipsoïde : ^^ = o ; donc 

V se réduit à 

M H„Y 



/• /•* 



et on conclut, en comparant les deux valeurs de V, que 

C — A 



H. 



On a donc l'équation 



M /Me, Ç-A_\ «Vî , i^v-Vi _^,, 

77""V'377 1?^)^i 3 ' 6 ~^ ' 

on en conclut que le coefficient de Y est nul, et que l'on a 



Me, C— A , wV? 



On a d'ailleurs 



3/-, 2/-Ï 6 



M 



&0 ,.2 ' 

'1 



M 



'1 



Il vient donc 



C — A 



0. il tDOi 



2r] "«" 3 6 ' 

rî ~ 3 r a/' 
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et, en supposant le rî^yon de la terre égal à l'unité, on a 



y 



C_A=^(e,-L) 



I^esanteur en un point de la surface. — Les composantes 

dU dU ôU 
de la pesanteur sont -^ — ,-^— , -r— , et on a 
* ojc oy oz 

^=\/(ïF(fRïr- 

La composante suivant le rayon OM est, s étant le cosinus de 
l'angle de la verticale avec le rayon OM, 

dU 

srcosz : 



àr 



i? ^^'^«' > 



comme e est petit, g cos e est sensiblement égal à g^ aux quan- 
tités du second ordre près, on a donc 

dU 

M ,, (C — A) Y 2coV wVY 



ty 



3 



Si maintenant on fait /• = i -|- Ç, 



7T='— 2^' 



on trouve 



2 03 



g=M{i-^^)+g,(e,-^)Y 



.2 ..2Y 



2(1)"' (O 



tr or I 



3 3 



é'o 



io8 
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C'est la formule de Clairaut quî donne la pesanteur en un 
point, en fonction de la latitude de ce point. 

Influence de F altitude, — Si on s'élève en ballon, la cor- 
rection est facile à faire, on n'a, dans la formule, qu'a rem- 
placer r par i + Ç+ /*, au lieu de le remplacer par i + Ç, mais 
ce n'est pas dans ce cas que l'on a à faire la correction d'alti- 
tude : le plus fréquemment cette correction se présente pour un 

lieu situé sur un plateau, et il faut 
tenir compte de l'attraction due à 
ce plateau. 

Il faudra donc tenir compte des 
aspérités positives ou négatives; il 
faudra d'abord choisir une surface 
définissant le niveau général du 
plateau : cette surface sera la sphère 
décrite sur OM comme diamètre 
(fig. i5), O étant le centre du sphé- 
roïde, M le lieu de l'observation. 
On calcule l'attraction due au vo- 
lume compris entre les deux sphères, 
l'erreur est peu considérable et d'autant moins que l'attraction 
des montagnes est très faible. On prend, en général, comme 
densité de la masse la densité moyenne des roches à la surface, 

c'est-à-dire D,, en représentant par D^ la densité moyenne 

du sphéroïde. 

Correction de Bouguer, — Supposons une couche sphériquc 
homogène d'épaisseur h et de densité p ; le volume de cette 
couche sera 

et sa masse 

47rA/7p. 

La correction apportée à la mesure de g sera ^Tzphrl et on a 

éT ~ 4 






■ô '^^''•' 



D,,- 
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D'un point M situé à une certaine distance au-dessus de la 
couche, menons le cône tangent à la sphère (fig. 16) ; le cercle 
de contact partage la couche en deux portions, et nous savons 
que les attractions de chacune de 
ces deux portions sur le point M 
sont égales. 

Appliquons ici ce théorème. 
Nous avons supposé que la surface 
de la planète était limitée à la 
sphère décrite sur OM comme 
diamètre ; si donc on suppose 
que l'attraction de la masse MRS 
se compose de l'attraction d'un 
certain nombre de couches homo- 
gènes et minces, la partie de ces 
couches limitée a la sphère (OM) 
aura sur M la même action que la 
moitié de la couche totale. Il faut donc calculer l'attraction sur M 
d'une couche entourant le sphéroïde de densité p et d'épaisseur h 
égale a l'altitude du point M. La correction à apporter à la 
pesanteur sera la moitié de cette attraction. 

L'épaisseur h de la couche peut se représenter par un dévelop- 
pement de la forme 




Fig. iG. 



le potentiel dû à cette couche, d'après une formule connue sera 

4t:AY 



1 



(2/1+ i)r 



,n +1 



l'accroissement de g sera 



ôg= 



IF' 



Zj ^ an + I ' 



car on fait ensuite /' = i dans la formule, et on a 



2K. 

C»0 



D^ Zj 2/1+ I 
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Si on se borne au terme n = o 

o«^ 3p 



si on en prend la moitié, c'est la correction de Bouguer. 

La correction totale comprend donc deux termes : i° celui 
qui provient de l'augmentation d'altitude sans tenir compte du 
plateau; 2° celui où Ton tient compte de la masse du plateau, on 
a finalement la formule 

S«^ ah 3 pA 



^0 '\ 2 D,r, 



Si on applique cette formule on trouve une valeur trop faible : 
la valeur fournie par l'observation est, en général, sur les conti- 
nents, inférieure à la valeur calculée ; elle lui est au contraire 
supérieure sur les océans et dans les îles isolées. 

M. Faye a ainsi trouvé qu'il vaut mieux ne pas tenir compte 
de la correction de Bouguer, ce qui correspond à des masses qui 
devraient se trouver sous les océans, et h des vides qui existe- 
raient dans les continents. 

A la surface des mers on a U = U^. La surface U = U^ est une 
surface qui, se prolongeant sous les continents, doit donner une 
surface de révolution, mais il n'en est pas ainsi. Cette surface 
s'appelle le géoïde. 

Pesanteur en un point de la surface du géoïde. — Considé- 
rons une sphère de rayon voisin du rayon de la terre ; soit Ç l'al- 
titude du continent au-dessus de cette sphère, on peut écrire 

Ç' la distance à l'ellipsoïde, peut s'écrire 

Ç — Ç' est la distance de la sphère à l'ellipsoïde au lieu consi- 
déré ; 
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Soit h Taltitude da point au-dessus du géoïde, 



"=1 



kY 



Î7 — h est la surélévation du géoïde au-dessus de Tellipsoïde : 
on a 



X' - h ="2^{ll -q-k)^, 



Ç — Ç' se réduit à deux termes : celui qui correspond h Y = i , 
et celui qui correspond a Y == 



,.2 



A la surface du géoïde U = ^^ == constante ; en un point d'al- 
titude &, on a 






àr 



On a d'ailleurs 



dU 



(r 



Donc 



_=.+2(H-A-)Y 

x = ,+y(ii_/,_v)Y. 

Les coefficients Y sont connus : ils correspondent au développe- 
ment de 



(O^ 



2^ 



i^' + f); 



si on rétablit l'homogénéité, on a 



I ^ S<H — /: — y ) Y 



O- 1. i.« + 1 



et comme on a : 






àr ' 
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on en conclut : 

I , s (n-)-.) (H- A— y) Y 



cr 7-2 ' ,.n + î 



Les opérations géodésiques font connaître A, les opérations 
de nivellement font connaître ^, les observations pendulaires 
font connaître g". Théoriquement, deux séries d'observations suf- 
firaient pour connaître Ç, A etg^; mais pratiquement les observa- 
tions d'aucune des trois espèces ne sont assez nombreuses, ni 
assez précises pour qu'on puisse se passer d'une série d'observa- 
tions. 

Jusqu'à présent Ton peut dire que l'on a trouvé la surface du 
géoïde peu différente de celle d'un ellipsoïde de révolution. 

Nous allons, dans ce qui suit, montrer qu'il peut exister des 
figures d'équilibre autres que l'ellipsoïde de révolution. 
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FONCTIONS DE LAME 



Coordonnées elliptiques, — Considérons les qiiadrlqueshomo- 
focales 

' -^ ' —1=0. 



Par un point de l'espace passent trois de ces surfaces : en efTet, 
x^ y^ z étant donnés, on a pour déterminer X^ une équation du 
troisième degré dont les trois racines sont séparées par les 
nombres a^, i^, c^. 

Il en résulte que Tune des racines, la plus grande, supérieure 
il a^y correspond h un ellipsoïde; la racine moyenne, comprise 
entre a^ etb^, définit un hyperboloïde à une nappe; la troisième, 
comprise entre b'^ et c^, définit un hyperboloïde à deux nappes. 

Soit p2, [JI.2, V" ces racines, on a donc 

p^ > a^ > [A^ > i^ > v^ > 6'^ 

Réciproquement si p, «jl, v sont donnés, on a trois surfaces se 
coupant en huit points, placés symétriquement par rapport aux 
divers plans de coordonnées. 

Si on ne considère que les points situés dans un trièdre trirec- 
tangle donné, les trois nombres p, a, v définissent un point et un 
seul. 

On peut écrire 

PoiKCARÉ. — Figures d'équilibre. 8 
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Si on multiplie par X^ — a^, et qu'on fasse ensuite \^ = a\ 



on a 



~ [a^ — U') {a' — c^) 



on obtiendra donc 



V la' — b^] id' — c^) 



et de même pour 3/ et ^, 



=v 



_ {^^ — p) {^' — b') [y' — b-) 

^— V (^b'—a^)(b'~c'} 

^ — V (t^ _ o») (c^ — 6^) 

Le signe qu'il faut prendre devant le radical correspond au 
trièdre des coordonnées que l'on considère. 
Soitfix) un polynôme, et posons 

R=/"(?^)> 

N=/-(v=), 

RMN sera un polynôme symétrique des trois quantités p^, [x^, v"^ 
raciïies de Téquation 



et par suite une fonction entière de o:^, y^, z^. 
Si on pose 

R = /-(P^)V/P"=^^ 
OU bien 

OU encore 



m m^ 
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et que Ton prenne pour M et N les mûmes fonctions de [jl et 
de V, RMN sera un polynôme en x'^^ y^,^^, multiplié par ^, par^r?/ 
ou par xyz. 

Supposons que R soit un polynôme entier en p* ; on peut le 
décomposer en facteurs du premier degré et écrire 

R=n(p«-Aî). 

Alors 

RMN = n (p= — >.f) [y>? — ÀJ) (v^ — Xî) 

RMN sera donc, dans tous les cas possibles, un produit de fac- 
teurs de la forme (i) multiplié par un certain nombre de quan- 
tités x^ y, z. 

Élément linéaire en coordonnées elliptiques, — Les sur- 
faces p =C*% [jL=C^*^, V = C*® forment un système triple ortho- 
gonal et l'élément linéaire peut s'écrire : 

ds' = oî'dp^ + p^diL" + rdy\ 

a, [3, y étant des fonctions de p, j^., v que nous allons déterminer; 
on a 

ds^ = dx^ + dy^ + dz\ 



dx pdp 1 ^du. 


v^v 


X p^ — cr [x^ — a^ 


1 » "^ 


dy pdp iKd[x 


v^/v 


y ?'-b^' ' \^'-b'' 


' v^ b' ' 


dz pdp 1 [JL<f[JL 


"/(h 


z p2 c' ' [Ji-^ c' 


V c' 


On a donc 





rf.«+rf,'+^.^=pvp^[^^+^.^,+^^] 



I 2^ 2 r__ _f^ I u^ I f__n 



^ "' |_(v-i_rt^j^ ^(«/•î—/.^)^ ^(vi—c-^j^J 



N 



1 16 FONCTIONS DE LAMÉ 

et ensuite 

g' (^»-p»)(.v»-p') 

p. -(p^_«^)(p^_Z,-^)(p«_c>) 

Je vais poser 

Q« = (i.'_p^)(p^_v»)(v'-[x»), 

(p'-«')(p^-&')(p-^-c-') 

p^ 



u 

4 

^ -z , 



A* et C^ sont positifs, B^ est négatif; donc A et C sont réels, 
B est purement imaginaire. 
Avec ces notations on a 



a 



p= 



V 



Q 

Bv/v^— p» ' 

Q 



cv/p^-i.' ' 

yv^ — p'est purement imaginaire, donc fi est réel. 

Equation de Laplace. — Voyons ce que devient Téquation 
de Laplace en coordonnées elliptiques. Nous avons déjà établi 
(p. 4^)j que si on désigne l'élément linéaire par 

ds^ = aVp2 + p'di»} + fdy\ 
on a 






L'équation AV = o devient donc : 



m^^h" 
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On a 
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El = ^ v/__J^lziz!__ 

a BC V (v'^ — p^) (?'•■ — [A-^) 

_Av/:~r(a^-v^), 



BC 
et Téquation de Laplacc est alors : 



S^[è(^'->')^]=- 



A seul dépend de p. Ou a donc 



2j B^ d^h-Ô^J^"' 



en multipliant par ABC on obtient l'équation 



W I[A"f- + A^f ](.■-., = «. 

Les coefficîents de cette équation sont réels, car A^, B^, C^ sont 
réels, et par suite leurs dérivées le sont aussi. 

Fonctions de Lamé. — Parmi les solutions de l'équation de 
Laplace, il y a des polynômes de la forme RMN : on les appelle 
fonctions de Lamé. 

Remarquons d'abord que nous avons les identités évidentes : 

^p^d.^— /) = o. 
Il en résulte que l'équation (2) est équivalente à la suivante 



l 
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Supposons que l'on ait V= RMN et que R satisfasse h l'équa- 
tion 

M et N satisferont à des équations analogues, et la fonction V 
satisfera à l'équation de Laplace. Cette condition suffisante est 
aussi nécessaire, c'est-a-dîre que si V est de la forme RMN et 
satisfait a l'équation de Laplace, R satisfait à l'équation (3), et par 
suite M et N satisfont aux équations analogues obtenues en rem- 
plaçant p par ijL et v, et A par B et C. 

En effet si on fait ul=v l'équation doit être encore vérifiée ; 
R se réduit a M en ce cas, et l'équation se réduit à 

B' -j-j- + B -3 3 h HuL^M + KM = o. 

Existe-t-il de telles fonctions ? Considérons la fonction p [u) 
définie par l'équation différentielle 

p'u = 2 ^{pu—e,){pu — e,) [pu — e^), 



ou 



Posons 



on en conclut 



^1 +^2 + ^3= o- 
p2 =pu-{-h, 

c' = ^3 H- A ; 



3 



On a donc 



= 2Ap. 

On a d'ailleur§ 

, du 



donc 
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' A = 



du 



et par suite 



-i 



, dR dR dp dR 



L'équation 



devient 



rfp dp du du 



^4'[^f-]+("^^^ + ^)^ = ^ 



^'?- + (II?^+K)R = o. 



du' 

Peut-on choisir H et K de manière que R soit un polynôme 
en p, ou un tel polynôme multiplié soit par y/p^ — a^, soit par 

y/p2 — ^2^ g^Ji- pj^j. y/p2 — c2p 

Ce polynôme, si on y remplace p^ par p (w), sera une fonction 
doublement périodique admettant les deux périodes de p[u) 
2(0^ et atOg. Elle admettra une double infinité de pôles qui seront 
les pôles de p (?/), c'est-à-dire u = 2/^0)^ + 2 nd),^^ m et n étant des 
entiers positifs, négatifs ou nuls. Ces pôles qui sont des pôles 
doubles de/? (m) seront pour le polynôme des pôles d'ordre 2/2. 

De plus y/p^ — a^ ^ \J pu — e^, est une fonction doublement 
périodique admettant les périodes 2C0j et 4w^ et ayant comme pôles 
simples les pôles de p (?/);[donc si on ajoute a u la valeur 20)^, 2C0j 
ou 20)3, la fonction R garde sa valeur ou change de signe sui- 
vant qu'elle contient o ou 2 radicaux, ou bien qu'elle est le pro- 
duit d'un polynôme par i ou 3 radicaux. 

Donc R est une fonction paire, ou une fonction impaire, mais 
dans un cas comme dans l'autre on doit avoir le développement 

R=^^+_i^H "h L. 

u^ ^ u^-' ^ u'^-' ^ ' 



si on porte cette valeur dans l'équation, en se rappelant que 



^ ^ u\ 
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on obtient l'équation : 

Pour que cette équation soit vérifiée il faut d'abord que le 
terme en ^ „ disparaisse, ce qui n'est possible que si Ton a 

II = — n [n -\- i). 

On a d'ailleurs si q est le degré de R, et qu'il faille lui ajouter p 
radicaux, [p = i, 2 ou 3) : 

Formation des fonctions de Lamé. — Ainsi donc les poly- 
f j, nomesj de Lamé, doivent s'ils existent satisfaire aux conditions 
suivantes. 

Regardés comme fonctions de ^, y, z ce sont des produits de 
facteurs quadratiques multipliés par un ou plusieurs des facteurs .r , 
?/, z. Si on remplace jF5 2^,js par leur valeur en fonction de p^,[jL^, v^,on 
obtient une fonction de la forme RMN, où R désigne, soit un poly- 
nôme en 0^, soit un tel polynôme multiplié par 1,2 ou 3 des 
radicaux 

v/p=^ \/f^^\ \/7^^^ 

et M, N désignent les mêmes fonctions ou p est remplacé par [jl 
et V. 

On a pour Q huit formes possibles, à chacune desquelles cor- 
(. respond une forme pour R. Si n est le degré de V le degré de R 
en p^ sera suivant le cas 



n n — I n — 2 n — 3 

_ 9 ~ 



, ou 



2 



suivant que R contiendra o, i, 2 ou 3 radicaux. 

Enfin si dans R on remplace p^par pu^ on obtient une fonction 
qui, lorsqu'on remplace u par a + 20)^, change son signe ou reste 
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invariable, suivant qu'il y a un nombre impair ou pair des deux 

radicaux y/p^ — ^^ety/p"^ — c", car \^pn — e^ ne change pas de 
signe pour w + 20)^. 

De même quand on augmente u de 2(ù.^,\/pu — e^et\/pn — e^ 
changent de signe et par suite R changera de signe ou non, sui- 
vant qu'il contiendra un nombre impair ou pair de ces radi- 
caux. 

Ces résultats sont résumés dans le tableau suivant, oii la pre- 
mière colonne contient la forme de Q, la seconde la forme corres- 
pondante de R, la troisième le degré du polynôme qui entre 
dans R en fonction du degré de Q, la quatrième et la cinquième 
indiquant si la fonction change désigne pour l'adjonction à udes 
périodes 2(0^ et aco^; Q^ désigne uniquement un produit de facteurs 
quadratiques. 




R ne contient pas do radical 
Rconticnt y P"^ — a* 

V^p- — c^ Y p-^ — a* 
V/(p^-a^)(p^-6^)(p^-c^) 



» 

M 
» 
» 
») 




+ 




±s 



I ' 



i) 

+ 



\ " + ï; 



2 ; 

n — I 
•1 



Détermination de K. — Ceci posé, cherchons à déterminer K. 
Pour les fonctions du premier genre, R est un polynôme 

d'ordre contenant 1- i coefficients homogènes. 

Substituons dans l'équation différentielle, nous avons un poly- 

n 



nome de degré 



I, qui devra être identiquement nul. Le 



coefficient du terme de degré le plus élevé disparaît de lui-même 
h cause de l'hypothèse II = — n (/i+ i). Il reste alors — 



+ I 



coefficients à faire disparaître. 



Les 1- I coefficients du polynôme doivent donc satisfaire 
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il 1- I équations linéaires et homogènes, dont les coefficients 

dépendent linéairement de K; l'élimination des coefficients du 
polynôme donnera donc pour K une équation de degré h i . 

Yl -1- I 

Si R est du second type, on aura h déterminer coeffi- 

cients d'un polynôme de degré •• Le degré de l'équation 

en K sera donc, d'après un raisonnement analogue au précédent, 

n-\- \ 



Pour ceux du troisième, l'équation en K sera de degré — . 

n — 
Pour ceux du quatrième, l'équation en K sera de degré - 



1 

Ces résultats figurent dans la sixième colonne du tableau pré- 
cédent. 

Il résulte de ceci qu'à chacune des racines des équations en K 
ainsi obtenues correspond un polynôme et un seul. 

Si n est pair, il n'y a de polynômes que de première et de 
troisième espèce ; leur nombre est 

n ^ n 

h I + J = 2/1 + I . 

2 2 

Si n est impair il n'y a de polynôme que de seconde et quatrième 
espèce ; leur nombre total 

' O ' :=2/l+ I. 



Donc il y a au plus in-\-\ polynômes de Lamé. Je dis au plus, 
car les équations en K pourraient avoir des racines multiples ; 
nous verrons qu'il n'en est rien. 

Analogie avec les fonctions sphériques. — Supposons que 
.r, t/y z soient infinis du premier ordre ; [ji et v restent essen- 
tiellement finis, mais p deviendra infini du premier ordre aussi. 

L'ensemble homogène des termes de degré le plus élevé dans 
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la fonction de Lamé, satisfait évidemment à l'équation de Laplace. 
C'est donc un polynôme sphérique que l'on peut écrire \\{x,y^z). 
Le terme de degré le plus élevé dans RMN en p sera ap^MN. 
Passons aux coordonnées polaires 

X = r sin 9 cos îp, 
y = r sin 9 sin cp, 
z = /*cos9, 

si p est très grand il diffère peu de r; u, v dépendront au contraire 
de cp et de 9. 

En effet l'équation qui détermine [jl peut s'écrire 

^^ v' ^' 

^ ' -I=:0. 



Si Xyy^z sont très grands, l'unité est négligeable; il reste une 
équation du second degré; [x et vies deux racines non infiniment 

grandes ne dépendent que du rapport-^ et —c'est-a-dire de cp 

et de 9. 

On peut alors écrire 

/'"Il (sin 9 cos cp, sin 9 sin (p, cos 9) = ap"MN, 

et comme /• = p, on a 

n (sin 9 cos cp, sin 9 sin (p, cos 9) = aMN ; 

le produit MN apparaît donc comme une fonction sphérique, 
etllseradécomposable en un produit de facteurs quadratiques (*). 
M. Moutard a démontré que c'était le seul cas où un polynôme 
sphérique était décomposable en un produit de facteurs quadra- 
tiques et de facteurs linéaires. 

MN, considérée comme fonction de 9 et de cp, est donc une fonc- 
tion sphérique : il en résulte que toute fonction de [jl et v peut se 
mettre sous la forme d'une somme 



s 



aMN, 



(♦) Voir Journal de Liouville, t. XI, 1846, p. 278. 
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puisque toute fonction peut se mettre sous la forme d'une somme 
de fonctions sphériques; mais l'analogie des fonctions sphériques 
et des fonctions de Lamé, va apparaître davantage quand on étu- 
diera le cas où les surfaces homofocales sont de révolution. 

Cas des ellipsoïdes de révolution. — Cela peut se présenter 
de deux manières différentes suivant que Ton a a = i, ou i = c. 

1° Si « = 6, on a des ellipsoïdes de révolution " aplati fi quand 
A" est supérieur à a^ et des hyperboloïdes h deux nappes si 
a'^<).^<6"^. Les hyperboloïdes h une nappe se réduisent à des 
plans passant par Taxe Oz, 

2° Si i = c, on a des ellipsoïdes de révolutions -allofrgés-autour 
de Ox^ pour ")} > cû. Si ar > A^ > i* on a des hyperboloïdes à une 
nappe, les hyperboloïdes à deux nappes se réduisant à des 
plans passant par Taxe Ox, 

Dans ces deux cas, soit [ji, soit v ne peut servir de paramètre, on 
est obligé de prendre comme paramètre l'angle du plan avec un 
plan fixe, le plan xOz par exemple dans les deux cas. Nous 
avons : 

x ^ {?' — a') {a^ — }JL^) (a^ — v' j ^ ù' — c' ' 
Si b tend vers a, u. tend aussi vers a, et l'on voit que l'on a : 



2 l.i 



— y- 



X ^ a^ — yr 

Si cp est l'angle du plan passant par le point [x^ y^ z) et Taxe Oz 
avec le plan xOz, on voit que [jl peut être considéré comme fonc- 
tion de cp, et l'on a 

]j} = a^ cos^cp + i^ sîn^9. 

On voit que l'on a 



limV ^2_{^2 =sincp, 






, — ces cp. 
a- — b' 
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Dans les mêmes conditions Téquatlon quî détermine v peut 
s'écrire en supprimant l'unité, négligeable devant x, f/ elz, * 

x^ + if z^ 

' ^= o. 



v^ — a^ v^ — c^ 



Si on passe aux coordonnées polaires, on a 

+ -"2 :t- = o ; 



/•" sin^O 



et ensuite 



sm 



cosO 



/•* 


cos^Q 




V 


c' 




J 


n:" 


v^ 


\ 


a' 


c' 


J 


V 


a' 



V 



Dans ce cas on peut immédiatement former les fonctions de 
Lamé. 

Nous savons que MN est une fonction sphérique ; M ne dépend 
que de îp, N ne dépend que de ; MN est donc Tune des fonctions 
sphériques fondamentales considérées plus haut (p. 4fe). 

On peut alors écrire M = cos jrf ou sin /^cp ; et pour Non 
prendra 



en posant ; 



fj:(cosO) = fj:[v/-^] 






On voit que les zéros de N se déterminent en calculant les 
zéros de Fj, abstraction faite des zéros ±: i. Ces zéros sont tous 
réels et distincts ; car, comme on le voit par le théorème de 
Rolle, les zéros de FJJ séparent les zéros de F^~^ On a 






Donc sî p est pair, Fj ne changera pas de signe avec sin 0; au 
contraire sip est impair, F^„^ changera de signe avec sin 0. 

Si on remplace :r, y, ^par leurs valeurs en coordonnées polaires, 



i'i6 
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Q changera de signe s'il contient Tiin des facteurs x ow y \ 
mais s'il n'en contient aucun, ou s'il en contient deux, il ne chan- 
gera pas de signe; p sera donc impair dans le premier cas, et 
pair dans le second. 

On voit aussi que y et sin p o changent de signe en même temps; 
donc si Q contient le facteur ?/, M devra être pris égal à sin/?3), 
sinon, M devra être pris égal a cosp'f. 

De même si i^ = c^, on prendra comme coordonnées polaires 

X = rcos 9, 

y z=z r sin 8 cos ï), 

;:; = /• sinO sin rs. 

La fonction N sera sin/;'^ ou posp'f, la fonction M sera 

Fi; (cos 0) . 

On prendra y; impair si Q contient l'un des facteurs y ou z; 
p sera pair au contraire si Q ne contient aucun de ces facteurs ou 
s'il les contient tous deux. On prendra pour N, sin y^ ç ou cosp'f 
suivant que Q contient ou ne contient pas le facteur z. 

Résumé, — Le résultat de cette discussion est résumé dans le 
tableau suivant : la première colonne contient les 8 formes possi- 



1 




a- : 

o 

I 
I 



I 
I 
o 
o 


- h'' 

Cos 
Cos 
Sin 
Cos 
Sin 
Cos 
Sin 
Sin 


^2 

o 
o 

I 
I 
o 
I 
I 
o 


-C2 

Cos 

Cos 

Cos 

Sin 

Sin 

Sin 

Cos 

Sin 


p 

p XJO. . . 

Pxj. . 

Px =. . . 
Px^j . . 

VXzx . . 
P X ^J . . 

Px^j-. . 


R 


RV^p=* a' 


RV^p2 b' 


Rj/p2_c2 


R\/(?'-^'J(p--0 


R^(p2_c2)(p2_tt2^ 


RV^(p*-a*;(p2-6^) 


R ^ (p2 _ a2j (p^ — 62j (p^ — C^) 
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bles de Q, la seconde contient les formes correspondantes de R, h\ 
troisième et la quatrième contiennent les cas limites pour a^ = ù^; 
dans la troisième on a mis le signe o si y^ est pair, i si p est 
impair ; la quatrième indique si M se réduit à cos p'f ou sin p'f ; 
la cinquième et la sixième contiennent les indications corres- 
pondantes dans le cas où ù^ = c^. 

Les valeurs de K sont réelles, — Nous allons maintenant 
démontrer que les équations en K ont toutes leurs racines réelles. 

Supposons d'abord a^ = P, On sait alors trouver les valeurs 
de K. Nous pouvons poser [x" == a^ cos" o + b^ sin* cp, puisque ^ 
est compris entre « et i et lorsque a = b^ ceci est encore la 
valeur de M ; l'équation qui définit \k 

peut s'écrire, en prenant comme variable (f 

dm . .., „, . du 



■(v-'-c')-j-r + 'r-a^)^in. 



COS 



"i 



d'^^ ' ^ ^ * 'do 



et si on fait a^ = b^ on trouve l'équation 

. [a' — c') -^ —[n{n + i)a' — Khl=o. 

Cette équation doit avoir comme solutions cos/;cp et sin/^'^, on 
a donc 

K — AZ(/1+I) 



— -> — =p^ 



L'équation en K a donc n -j- i racines réelles, correspondant 
aux valeurs de /?, 

o, I5 2, . . . (/i — i), n. 

On voit que les équations en K ont toutes leurs racines réelles ; 
les valeurs paires de/;, séparent les valeurs impaires. 
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Lorsque b variera de a vers c les racines des équations corres- 
pondant aux deux groupes de fonctions de Lamé, contenant 
zéro et deux radicaux d'une part, sépareront les racines corres- 
pondant aux deux autres groupes contenant un ou trois radicaux ; 
cela ne cessera d'être, que si une des racines d'un groupe d'équa- 
tions, devenait égale h une racine \ d'un autre groupe. 

Soient, K cetle racine, R et S les fonctions de Lamé corres- 
pondantes, appartenant par hypothèse a deux groupes dille- 
rents : on a par hypothèse 

R/'4-(IIp2-+-K)R=o, 

S''H-(Hp^ + K)X=o; 

S 



on en conclut 



R*S*<^S''R=o. 



On peut intégrer immédiatement et écrire 

nfs* — S'R = a. 

Si on suppose R et S exprimés en fonction elliptique du para- 
mètre //, on sait que l'on pourra augmenter u de 2 (^, de 2 co^ ou 
de 2 CJ3 de façon que R par exemple ne change pas de valeur, et 
que S change de signe, et de même pour R' et S^ ; mais alors le 
produit RS' — SR' aurait changé de signe ; on a donc a = o ; 

on en conclut que -^ est une constante, ce qui est évidemment 

faux. 

Deux équations en K correspondant à une même valeur de n 
ne peuvent donc admettre de racine commune. 

Remarquons que dans l'équation en K le coefficient du terme 
de degré le plus élevé se réduit à l'unité; il est donc impossible 
que le degré de l'équation s'abaisse. 

Les fonctions de Lamé sont linéairement indépendantes. 
— Je vais démontrer maintenant qu'entre les p fonctions de 
Lamé, d'une même forme, et pour une valeur de n donnée il 

n'existe aucune relation linéaire à coefficients constants ; (/^ est 

V j u ^ i ^ — i n-\- I n\ 

1 un des nombres \- i. , ■ , — ) . 

2 2 22/ 
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En effet : supposons que Ton puisse avoir une relation de la 
forme 

en portant R^+i dans Téquation du second ordre, on aura à écrire 
que les p coefficients d'un polynôme seront nuls ; on aura pour 
déterminer a^ a^... a^. p -\-i équations du premier degré, dont la 
première sera satisfaite identiquement. Il restera donc p équa- 
tions homogènes à q inconnues. L'élimination de ces inconnues 
donnera un certain nombre d'équations en K de degré au plus 
égal à gr et ces équations devront être satisfaites pour les p va- 
leurs de K ; ces équations se réduisent alors à des identités. On 
peut donc déterminer les coefficients a^, ag... a^ quel que soit 
K, et on aurait alors un polynôme comme solution de l'équation 
de Lamé, et cela quelque soit K; c'est une contradiction. 

Les polynômes de Lamé ont leurs racines réelles. — Je 
vais maintenant démontrer qu'un polynôme R quelconque a 
toutes ses racines réelles et comprises entre a^ et c^. 

On peut partager l'intervalle de — 00 à + 00 en quatre autres 
par les nombres a^^ h^ et 6*^. 

Je vais d'abord démontrer que le nombre de racines compris 
dans chaque intervalle est constant lorsque i* varie de a^ à c^. En 
effet, il ne peut en être autrement que si, pour une valeur de i, 
deux racines réelles devenaient imaginaires. Mais alors, pour 
cette valeur de i, ce polynôme aurait une racine double : R' serait 
nul. On a 

R'' = FR = o, 

F étant (IPp + K) ; on a 

R'' = FR'+F'R, 
R^v ^ FIV'+ 2F'R'+ r^R'^ etc. ; 

toutes les dérivées seraient nulles, ce qui est impossible. 

Il pourrait encore arriver que Tune des racines devienne égale 
à a^^ h^ ou 6"^. Supposons par exemple qu'une racine devienne 
égale à a^ : comme on a posé 

pn—e^=f—a\ 
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la valeur de u correspondante serait co^, or on ay/w^ =0, o), est un 
"^ zéro de p'^^ ; donc ce serait un zéro de R' et l'on revient au cas 

précédent. 

Ceci posé soit s^ le nombre des racines imaginaires de l'ensem- 
ble des polynômes correspondant à une valeur de n ; 

£2, le nombre des racines réelles comprises entre «' et i* ; 

Eg, » » réelles comprises entre i* et c* ; 

e^^, » » réelles non comprises entre a^ et c^, 

. Le nombre total des racines est égal à la somme des degrés des 
polynômes 



n [ n \ , 

2 \ 2 / 



^4- 



Lorsque i* est égal à a', toutes les racines s^ se réduisent à a^ 
et e'3 se réduisent à a* ou à c^ ; lorsque h^ est égal à c^, toutes les 
racines Sg se réduisent à c^ et e^ se réduisent h a^ ou à c^. 

Mais si a^ = b^ les polynômes sont les dérivées successives 
de (v^ — i)" et ils ont toutes leurs racines réelles ; il y en a un qui 

a o racine, un qui en a une, un qui en a deux, ... un qui en a — 
abstraction faite des racines :± i ; on a donc 



de même on a 



n / n \ , 



Nous avons donc l'égalité 



^1 + Ê2.+ £3 + £4 =-7- (-^ + ^)+ ^1 +^4 + 2^2+ S'3; 



d'où l'on conclut : 



et par suite 



ei+Sii + e'2 + e's = <>- 



e,= e;=£'g = £\=o. 



Le théorème est donc démontré. 
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Pour savoir maintenant combien il y a de racines dans chacun 
des intervalles a^... b^, i^...c%on voit que dans le cas général il y 
en a autant que dans le cas particulier, où a^ = b^. 

Et on voit qu'il y a : 

un polynôme ayant o racine entre a^ et i^, et — entre b^ et c^. 

n 

» I » — -*-. I M 

2 

n 

)) — » o » 

2 



Développement d'une fonction en une somme de fonctions 
de Lamé, — Soit un ellipsoïde de la famille, déterminé par le 
paramètre p ; si je pose 



v/(p^-i-"/(?^-vV 

j'ai 

rZMN.M,Nycr = o. 

MN désignant un des produits déjà désignés, M^N^ en désignant 
un autre, et Tintégrale étant étendue à toute la surface de Tellip- 
soïde . 
Soit 

V = RMN, V, = R,M,N,. 
Puisque Ton a AV = o, la formule de Green donne 




V-^— V — W<r = o 

dn ^ dn / 



l'intégrale étant étendue à une surface fermée quelconque ; sup- 
posons que cette surface soit l'ellipsoïde ; R et R^ sont alors des 

dV d\ 
constantes. Calculons -^— et —j-^ : quand on se déplace sur la 

an an 
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normale à Tellipsoïde, ix et v demeurent constants, p augmente 
de dz ; et on a 



mais nous avons 



rfV 
dn 


ÔV rf? _ dV du 
dp dn du rfp 

ÔV du 
du dn ' 

dV -,„ dR 


d? 
dn 



du du ' 

l'intégrale devient donc 




MN. M,N, (RR', — R,R') ^ rfa = o ; 



mais RR'^ — R^R' ne dépend que de p et comme c'est une cons- 
tante, on peut la faire sortir du radical. On a d'ailleurs 

il vient donc en changeant le signe 

— MNM.N, ~dfj=o. 

Tout revient à calculer 

du du do 




dn dp dn ' 

nous avons vu que l'on avait, en posant p^ =pu'\- h 

d'où 

do , 

du - ' 

lorsque p tend vers l'infini, pu tend vers oo , u tend vers o, donc 
décroit ; on doit prendre le signe moins. 



I 
'il 
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On a d'ailleurs, puisque rélément linéaire est 

dn^ = ol' df + pWjjL^ + f rfv^ 

et que p seul varie sur la normale : 

d^ ^j Av/(iJL-^ — >/) 

dn OL Q 

On a donc, en posant 

du 



1 = 



dn ' 



+ I 



ce qui est l'égalité annoncée. L'intégrale, d'ailleurs, ne dépend 
pas de p. Car nous avons 

rfo" = Pyd^dv ; 
et par suite 

w _ Q^ y/— I rf[JL^v 

\^~BCv/(p'^-^o'^-v^) 

/ ,._ ([^^ — ^^) V^— ' ^!^^^ 
t^^'~ BC 

On a démontré que toute fonction déterminée sur l'ellip*^ 
solde peut se mettre sous la forme d'une somme de fonctions 
de Lamé; 

* =^ A4 Ma.N,. 

Le théorème précédent permet facilement de trouver les coef- 
ficients. On a 

fl<^Mk^kdcr^fYAJ^hNM^id<T; 
le second membre se réduit a A^t / M^N^rfo" : ce (][ui détermine A,,. 
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Les fonctions S. — Nous allons maintenant considérer une 
autre fonction de p satisfaisant aussi à Téquation de Lamé; nous 
la définirons ainsi qu'il suit. 

Soit R une fonction de Lamé de degré n définie par la relation 

R"+FR = o; 

la fonction S, associée à la fonction R, sera définie par la relation 

S" + FS = o 
où F a la même valeur : on en conclut que 

S''R — R"S = o; 

en intégrant on trouve 

S'R — RfS^=C% 

et nous prenons arbitrairement cette constante égale à an -j- i ; 
on a donc 

S^R — R^S^ _ 2/1+1 

R^ — R2 • 

Donc 



S / 2/1+1 



du ; 



on peut choisir arbitrairement la limite inférieure d'intégration 
et on a 



S I 2/1+ I , 





Pour p^ très grand 

R = Ap», 

A étant une constante que nous pouvons prendre égale à Tunité ; 
Il est alors très petit : puisque 



LES FoycriONS s lis 



p^ est approximativement égal à — j 

I 

p » » — 



On a donc 

R = z/-« , 



— -V- 



S=«-" / T. A.-»*'; 





s est donc sensiblement égal à z^** + * = ^^^ ; le produit SMN 



satisfait aussi à Téquation de Laplace. 
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Problème de Diricblet sur r ellipsoïde. — Abordons main- 
tenant le problème de Dirichlet pour un ellipsoïde, dans lequel 
on se donne, sur un ellipsoïde déterminé par le paramètre p^, 
les valeurs d'une fonction V, en proposant de trouver la fonction 
V satisfaisant h l'équation de Laplace et prenant sur rellipsoïde 
la série des valeurs données. 

Soit sur Tellipsoïde 

V-=«'=2a,M,N,; 

on peut poser 

A,=«,RÎS?, 

RJ étant la valeur pour p = p^ de la fonction R associée à M et N, 
et Sj » » S, » R,. 

On a alors 



^=Va,R?S?M,N,. 



Ceci posé, considérons la fonction définie par l'expression 



y .s. 



RMN 
à l'intérieur de rellipsoïde, et 

aR«SMN 



!■ 



a l'extérieur du même ellipsoïde. 
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Ces deux séries sont convergentes dans les espaces indiqués, 
la fonction est donc bien déterminée. 
Calculons les quantités 

(")V ÙV 

et 



on a 

= ;y^a,R?S?M,N,. 

On a de môme 

ùV 



0/2 



- = — ;V a,S?R'?M,N,. 



V est une fonction continue quand on traverse la surface, mais 
ses dérivées ne le sont pas. V est donc le potentiel d'une masse 
superficielle, dont la densité S est déterminée par la relation 

ÔV ùV _ g 



Au lieu de supposer la densité variable, on peut supposer une 
distribution de densité égale à i, mais d'une épaisseur î^, telle 
que Ton ait 

— — — — _/-r 

équation qui détermine Ç, car on a 

471?; =-}- ?^a,(R?S^? - R'?S?)M,N, 



= zVa(2«+i)MN. 



Si l'épaisseur de la couche est déterminée par une fonction de 
la forme 



2 /?MN , 
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on en déduira réciproquement le potentiel sur la surface par la 
formule 



V-y — 4— MN; 



d'où Ton déduit simplement le potentiel à l'intérieur et h l'exté- 
rieur de la surface. 

Application» — Voyons maintenant les plus simples des fonc- 
tions de Lamé. 

I® pour n==o on a 

S= li. 
2^ pour w = I on a trois fonctions R 



On a 



y = /t,R,M,N„ 



en posant 



A. = 



' V/(o^ — Z.=')(a*— f')' 



h,= 



\J{b-'—c'){b'— à') ' 



K-=-^ ' 



\/(c^ — a') (6-«— b^) ' 

3° Pour /i = 2, on a : 
d'abord deux polynômes du second degré 

p2 — A^ 
et 

p^ - A'^ 
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que Ton peut calculer facilement ; on sait que Ton a 

a^>h'>b'>h">c\ 
et ensuite les trois fonctions 

On a 

R^ = R2R3, 

R, = RjR2; 

le produit R^^M^N^^ est à un facteur constant près i/z. 
» R^MjNj » » zx. 



» 



RjMjNj » )) xy, 



Remarquons encore la signification des fonctions R^, R^, R3, 
qui sont les demi-axes de l'ellipsoïde correspondant à p. Son 
volume est 

T = -77- TrR^RjRg = — -TiR^R^. 

Supposons qu'à la surface de l'ellipsoïde, soit répandue une 
couche de densité constante, et d'épaisseur 

le potentiel à la surface sera 

les indices inférieurs indiquant le numéro d'ordre de la fonction 
de Lamé, l'indice supérieur indiquant qu'on y fait p = p^. A l'in- 
térieur, le potentiel sera 

le potentiel sera constant à l'intérieur ; on peut donc consi- 
dérer la distribution de matière comme la distribution d'une 
couche électrique en équilibre à la surface de Tellipsoïde sup- 
posé conducteur. 
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Le potentiel sera a Texténeur 

Les surfaces équipoteiitlelles seront donc les surfaces u = C^*, 
c'est h dire les surfaces p == constante ; ce seront des ellipsoïdes 
homofocaux. 

La couche correspondante est celle comprise entre deux ellip- 
soïdes homofocaux comme le montre la valeur de Ç. 

Attraction d'un ellipsoïde homogène, — Considérons (fig. 17) 
un ellipsoïde homogène dont le potentiel soit V. Déplaçons l'el- 
lipsoïde parallèlement à Taxe Ox^ d'une 
quantité c; le potentiel en un point [jCy y, z) 
P' va devenir 




y (^ — e, î^, z)=V {x,y,z)—t 



dV 

^x 



et l'ellipsoïde prendra la position E'. La 
différence de potentiel provient de ce qu'il 
y a en plus la région ombrée h droite, et 
Fig. 17. en moins la région ombrée à gauche. On 

peut supposer que c'est le potentiel d'une 

épaisseur Ç répandue a la surface de l'ellipsoïde, Ç étant positif 

à droite et négatif à gauche. 

Calculons !^ : si PP'' est la normale commune aux deux ellip- 

soides, Ç = PP"; PP' = e, et le triangle PP'P'' est rectangle 

en P'', on a donc 



PP''=£cos(PF,PP'0. 



Comme 



FF' cos (PP, PP'Q =1 projection de PP'/ sur 0.r, 
et, puisque PP'' est normale h l'ellipsoïde p, on a 

J^ cos P = -^ do 
- dp ^ 

^x dp , 
— ^ dn 



dp dn 

^x dp ,. 
dp dn 
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de ces égalités on déduit 

^ dx dp àx du , dx 

cos P = -T r- = -T 7— == — l 



Op dn du dn du 



Nous savons aussi que : 

dx 



du ' 
en posant 



=/ji;m,n„ 



^, ^ DR^ dp 

* dp du 







et comme 

A = -f V{?'-^n{?'-à'){?'-c^) = ^R.R.R, = -^ R.R> , 

r r r 



on en conclut 



et l'on a 



R'. = R„ 



^■^ =/»R,M,N,. 



du 



Il vient donc 

cosP=yr;m,n,, 

On peut donc maintenant déterminer le potentiel V relatif à 
la coucha d'épaisseur 2^, en appliquant la formule connue : on a à 
la surface 

V'=4-4^ £/^jsîr;m,n,. 

o 
A l'intérieur de l'ellipsoïde on a 

V'=-f-i^£R»S;R.M.N, 
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et, puisque R^M^N^ =x^ 



on a donc 



dV 

do; 



Rî 



A Textérieur on a 



|1 = i^ R!R?S,MA 
~~ R. ' 



dV 



A rintérieur -r— est donc proportionnel à a;, k l'extérieur la 

relation est plus compliquée; -r— représente la composante de 
l'attraction parallèle a Ox, 

Théorème d'Yvory. — Considérons deux ellipsoïdes homo- 
focaux, correspondant, le premier E^ au paramètre p^, le second Ej 
au paramètre p^; Pj > Po (fig* i8). Soient P^ et P^, deux points cor- 
respondants, c'est à dire ayant 
respectivement pour coordon- 
nées elliptiques (p„ a, v) (p„ 
[JL, v) . Les coordonnées rectan- 
gulaires de deux points corres- 
pondants sont proportionnelles, 
c'est-à-dire que Ton a 




X. 



=px 



1» 



yo = qvv 



Fig. i8. 



rz 



i> 



/?, y, /• étant seulement fonctions de p^ et de p^. 

Le théorème d'Yvory consiste à comparer l'attraction de l'ellip- 
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solde E^ sur P^ a Tattraction de Tellipsoïde E^ sur P^. La pre- 
mière a pour composantes suivant Ox 

la seconde a pour composantes suivant Ox 

^47rARîR;S}M,N,, 

le quotient est donc 

RI ~^ ' 

on trouverait de même que les rapports des attractions sui- 
vant Oy et 0^ sont respectivement 

R° R° 

^' et ^' 



Ri RJ 

Ellipsoïde de Mac Laurin. — Considérons un ellipsoïde 
homogène, que nous faisons tourner autour de l'axe Ox, d'une 
rotation infiniment petite o) ; le point de coordonnées x, y^ z, 
aura pour nouvelles coordonnées 

^, y + wr, z — ioy; 

son potentiel sera 

N Ar/ N . / ^V dV\ 

V(^, y + -w, z—ijio)=Y{x, y, -J+^Y-^ ^'d7/' 

V — V est le potentiel dû à la couche superficielle comprise 
entre les deux ellipsoïdes. 

Si {x^ y, z) représente un point intérieur, l'action en ce point 
a pour composantes 



T:rS; TySl TcS 



r; ' R^î ' R 

elle dérive d'un potentiel 



"> 



^•=^.-i(-'-l+^*-|+''-|)' 
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en appelant V^ le potentiel au centre ; et cette formule est encore 
vraie à la surface : donc 

dV dV 

et 



sont proportionnelles h y et h z, et on a 

V — V = A(ûyz 

= A(oR,M,N,. 

En appliquant la formule démontrée précédemment (p. 137), 
Tépaîsseur de la couche en un point est 



47r Sî 

Pour que l'ellipsoïde puisse être une figure d'équilibre, il faut 
qu'à la surface ou ait 



0)2 



VH {y'-^z') = C\ 

Nous allons traiter un problème plus général ; de savoir si un 
ellipsoïde à trois axes peut être une figure d'équilibre, sachant 
qu'il est soumis, outre les attractions mutuelles de ses parties, à 
une force dérivant d'un potentiel 



a, [3, y étant des constantes. 

Pour que l'équilibre ait lieu, il faut que 

soit constant à la surface de rellipsoïde dont Téquation est 

x^ 7/^ z^ 

1 — 1 I = o 

Rî ^ R^ ^ R^ 



Il faut donc que Ton ait 
(I) «RJ-TS.R, = ?R^-TSA=TRi-TS3R3. 
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Dans le cas particulier qui nous occupe, c'est-à-dire quand 
Tellipsoïde est en rotation, on a, a= o, ^ = y = to^, et nos 
formules deviennent 

^ TS^R, = (o^R^ — TS,R, = iù^Ki — TS3R3 ; 
on a donc 
(II) ^^ =^\ ^^383 — RjSt ^V R^Sa — R,S, 

1 \ J^8 A'-J 



et d'ailleurs 



T= -^ RjR^Rg; 



T est une donnée du problème où les inconnues sont p^ — a^, 
^—b\ f — c\ 

On a une première condition, qui est que a>^ doit être positif; 
donc il faut que 

R,S, — R,S,>o, 

R3S3— RiSi>o. 

Ces inégalités sont toujours vérifiées. Démontrons par exemple 
que la première Test, c'est-a-dire que Ton a 

R2S2— RA>o. 

Remarquons que R^, R^, S^ et Sg sont des quantités positives, 
puisque p^ > a^ ; il faut donc que Ton ait 



Si Rj ' 



c'est-à-dire 
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et on est ainsi ramené a un théorème d'arithmétique connu. Le 
premier rapport est compris entre la plus grande et la plus 
petite valeur du rapport des quantités placées sous le signe J": la 
plus grande de ces valeurs est i correspondant à m = o, p = oo, 

p2 — a^ 
la plus petite est ^ j^ ; le théorème est donc démontré. On 

en donnera plus loin une démonstration plus générale. 

On voit de suite que Taxe de rotation ne peut pas être Taxe 
Oy, car il faudrait que RgS.^ soit la plus grande des trois quantités 
analogues. 

La relation (II) a une solution évidente, c'est R2 = R3, qui con- 
duit à un ellipsoïde de révolution autour de Ox : c'est l'ellip- 
soïde de Mac Laurin. 

Ellipsoïde de Jacobi, — Les relations (II) peuvent aussi être 
satisfaites pour un autre ellipsoïde découvert par Jacobi, mais 
nous devons, avant de le rechercher faire la remarque suivante : 

Si nous nous donnons un ellipsoïde homogène quelconque, 
nous pouvons toujours choisir les constantes a, p, y, de manière 
à satisfaire aux relations (1) c'est-à-dire de façon que cet ellip- 
soïde soit en équilibre sous Faction de la gravitation et d'une 
force supplémentaire dont les composantes sont aj7, ^y^ yz ; nous 
allons supposer a = o. 

La condition d'équilibre d'un ellipsoïde homogène E soumis 
h la gravitation et a une force perturbatrice dérivant d'un poten- 
tiel — (Py^ -}- Y^^) est que, à la surface 

soit une constante. 

Supposons que l'ellipsoïde se déforme mais peu, de façon que 
la nouvelle surface S ne diffère pas beaucoup de l'ellipsoïde, et 
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cherchons à quelle condition 2 sera une figure d'équilibre, en 
continuant à supposer que le fluide reste soumis aux mêmes 
forces, c'est-a-dire à la gravitation et à la force a^, ^i/, yz. 

La condition d'équilibre est que, en désignant par V -|- f le 
potentiel dû à la figure 2 ; 

soit une constante à la surface. 

Soit ^ la longueur de la normale comptée de E à 2 ; on pourra 
développer ^ en une somme de fonctions sphériques : 

et Ton aura à la surface 

4irP 



=s 



R^S^MN. 



2n-\- i 
D'autre part, on aura, a des quantités près de Tordre de J^', 

et, en appelant g' l'intensité de la pesanteur à la surface, il vient 

dU 



.<r 



d'où 

La condition d'équilibre est que U soit une constante et comme 






u==v+^+— (?z/'+T-')=c^ 



la condition d'équilibre est que 
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uni une constante à la surface c'est-à-dire que 



(i, 



Y, ^êr ^*^*^'^' -S^'?^»^' = ^' 



Je dis que gl est une constante sur Tellipsoîde ; en effet, 
déplaçons Tellipsoïde parallèlement à l'axe 0.r, Téquilibre ne 
sera pas troublé, car par suite de la condition a := o, le travail 
de la force (o^r, ^y, y^) ^^^ i*ul- Nous avons vu que l'on avait 
pour t^ la valeur 

K étant une constante ; l'équation précédente se réduit en ce 
cas à 

ceci ne sera constant que si 



L'équation (i) devient alors 



KiS^-^^--^ 



et Téqualion d'équilibre devient, en supprimant Tindice zéro, 

R,S, R.S. 



2/« + I 3 



o. 



Si la masse est en équilibre sous l'action de l'attraction et de 
la force centrifuge, elle le sera encore si on lui donne un mouve- 
ment de rotation. Cette rotation est une des transformations que 
l'on peut faire subir à un ellipsoïde de Jacobi, sans troubler 
ré((uilibre. 

Nous avons vu que dans ce cas, on avait à un facteur constant 
près (p. i44) 

!;=M,N,; 
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on en conclut que rellipsoïde de Jacobi doit satisfaire à la relation 

R,S, R,S, 



6 

Réciproquement, supposons qu'un ellipsoïde soit en équilibre 
sous l'action de l'attraction et d'une force dérivant du potentiel 

je dis que si on peut le faire tourner d'un angle a> sans troubler 
l'équilibre, c'est que P=Y. 

En effet, a la surface on doit avoir, en appelant V le potentiel 
newtonien 

V + -^(P^/^4-v^^) = V, = C"; 

au sommet du petit axe V = Vj; si on fait tourner l'ellipsoïde 
d*un an^le o) autour de 0^, le sommet du petit axe restera sur la 
surface de l'ellipsoïde, la constante V, n'aura pas changé, V n'aura 
pas changé non plus. 
Or 

il faudra donc que sur la surface 
et en tenant compte des relations 

ot/ = (OZ 
oz == — toy, 

on devra avoir sur la surface 

(? — T) wî/c = o; 
donc 
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Discussion de l'équation des ellipsoïdes de Jacohi, — Nous 
avons donc à discuter Téquation des ellipsoïdes de Jacobi 

Cette équation ne convient pas aux ellipsoïdes de révolution, 
cela provient de la façon dont on a établi cette équation : dans le 
cas d'un ellipsoïde de révolution, on aurait en effet Ç=o. 

Au lieu de discuter l'équation (i), nous allons discuter l'équa- 
tion plus générale 



2W+I 2/Z-t-I 

oii l'on a p" > d^ ; m étant l'ordre de la fonction R^., n étant 

l'ordre de la fonction R,. 

F 
Fet-rp- sont de même signe, de plus R;^ ne peut pas s'annuler 

puisque les racines des polynômes R sont comprises entre a~ 
et c^ ; donc l'équation 

p .._ F _ S, S, R? 

'- RJ (2«i+i)R, (2n+i)R, RJ. 

a les mêmes racines que F=o. 

F 
Les racines de l'équation ^.^ = o sont séparées par celles 

de la dérivée : 

^F, I d S,, I R? d S, 



du im-\-i du R,^. iin-\-i R* du Rj 

I S, rf / R? 



2/2 -f- I R; du 



m 



On a par définition 



d Sj^ 2/72+1 

'dIi:'Ri^~ Rl~' 
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et par suite 

d?, 1 S, R, R;R;,— R,RI 



du 2/Z+ 1 R, ' R/. R 

aS, r;r,-r,ri. . 



» 



2/2+1 R 

mais S,- ne peut s'annuler : donc les racines de , ^ sont les 

racines de RÎRfc — RiR/t. Nous ne nous occupons évidemment 
que des racines supérieures h a^. 

Cherchons la dérivée de cette fonction par rapport à «, c'est : 

R;'R^ - R.Rl' ; 

et, d'après la définition de R, et R;>., 

R;/ = (IIp^ + K) R,, 
on a 

R^R^ *[«(«+,)_ ,„ [ni + 1 )] p' + K, — Kj.) . 

Les deux premiers facteurs n'ont pas de racines en p^ supé- 
rieures à a^; le troisième ne peut évidemment en avoir qu'une 

seule : donc ^ ne peut en avoir plus de deux, F^ ou F ne peu- 
vent en avoir plus de trois. 

R. 
Supposons que -^ soit constamment croissant, sa dérivée ne 

Ra 
s'annulera pas ; donc F^ sera constamment décroissant et ne 

pourra avoir qu'une seule racine. Au contraire, si — - est toujours 

décroissant, Fi sera toujours croissant. D'ailleurs pour u = o, 

p =00 , -^êst à peu près égal à— ; donc F part de zéro et est 

toujours croissant ou toujours décroissant, donc n'a pas de 
racines, sauf la racine p = QO qui n'est pas admissible. 

Dans le cas le plus général, F peut avoir, nous l'avons vu, 
trois racines : mais parmi celles-là se trouve la racine p =qo , 
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,R, 



F a donc zéro, une ou deux racines supérieures à a*; si-ijî- est 

toujours croissant ou toujours décroissant, elle n'en a aucune. 
Revenons alors à Téquation 

O 2/1+1 

(R\ ^ 
-q-* ) sera un polynôme 

dont toutes les racines sont réelles et comprises entre a^ et 
c^ ; donc la dérivée ne peut s'annuler pour une valeur de p^ supé- 
rieure à a^. L'équation F n'a pas de racine supérieure à a*. Sup- 
posons maintenant que R, ne soit pas divisible par yç^ — a% et 
que n soit plus grand que i, je dis que l'équation a au moins une 
racine. En effet, substituons a^, le premier terme s'annule, le 
second est négatif : on a donc une somme négative ; remplaçons 
p par 00 ; le premier terme de l'équation est a peu près égal à 



p V 3 2/1+1 / 



le résultat est positif; il y a donc un nombre impair de racines 
entre a^ et + oo , il ne peut y en avoir plus de deux, il y en a donc 
une et une seule. 

Ce raisonnement n'est pas valable si /i= i, c'est-a-dire si 

ou bien 

dans ces deux cas, la quantité 

^_l 

et les quantités analogues sont toujours croissantes lorsque p^ 
varie de oo à a^. Donc l'équation correspondante n'a pas de 
racines : c'est un théorème déjà énoncé (p. i45). 



,2 
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Appliquons ce prîncJpe h rdlipsoïde de Jacobi. Il faut consi- 
dérer l'équation 

\ i 1 3L. = o • 

3 5 ' 

R^ n'est pas divisible par R^ : donc il y a une racine, et une seule. 
Donc parmi les ellipsoïdes homofocaux, il y a un ellipsoïde 
et un seul, qui soit une figure d'équilibre ; c'est l'ellipsoïde de 
Jacobi. 

Figures dérivées de l'ellipsoïde de Mac Laurin. — Suppo- 
sons que b^ = c^ et cherchons s'il y a des figures peu différentes 
de l'ellipsoïde de Mac Laurin. Pour qu'il existe une figure d'équi- 
libre 2 peu différente de cet ellipsoïde, il faut que l'on ait la 
condition 

R,Sj RfS, 



(0 



3 2/l-|- 1 

Nous avons vu qu'en posant 



on avait 

N,== cos/?cp, 



M 



-'^(v^:-^)- 



h étant une constante. On a posé 

F (i)= (i — i^)i D»» + ^' (i — i7\ 

La condition nécessaire et suffisante pour que l'équation (i) 
ait une racine est que R, ne soit pas divisible par y p'^ — a^'^ donc 
Mf ne doit pas être divisible par y a^ — [jl% donc n-\- p doit être 
pair, n etp étant déterminés, une figure d'équilibre voisine d'un 
ellipsoïde de révolution est déterminée par la relation : 

Ç = £ IMi cos po 

ou, d'une façon plus générale, 

J^ == IMi (e cos p'f T-}- e' sin p'f) . 
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On peut supposer £' = o : cela revient à faire tourner les axes 
de coordonnées d'un angle convenable autour de Ox, La figure 
que l'on obtient ainsi est une figure admettant l'axe Ox comme 
axe de symétrie d'ordre/^. F (i) ne change pas si Ton change Ç en 
— Ç, donc le plan des j/z est un plan de symétrie. Si/> = o, Ç ne 
dépend pas de îp, la figure est de révolution. 

On ne doit pas conserver parmi les équations obtenues celle 
où n = Oy car le volume ne serait pas conservé; de même on doit 

exclure n= i, car alors R,- serait divisible par y p^ — a^ : donc 
la valeur minimum de n est 2; à cette valeur de n correspondent 
deux valeurs dep : ^ = o, et /? = 2. 

Quelles sont celles de ces figures qui sont ellipsoïdales : con- 
sidérons un ellipsoïde E^ et un autre E, qui en soit peu différent 
mais non homofocal. Le potentiel de chacun d'eux k son inté- 
rieur est une fonction du second degré des coordonnées carté- 
siennes ; la différence des potentiels est donc aussi un polynôme 
du second degré : c'est d'ailleurs le potentiel dû a la couche com- 
prise entre les deux ellipsoïdes. On petit par suite développer ce 
potentiel sous la forme 

YsRMN, 

R, M, N étant des fonctions ^ph é rique ^ d'ordre o, i ou 2. 
A la surface on a 



-s 



sR.MN. 



et à l'extérieur on a par conséquent 



^=S ' X ^^*^' 



l'épaisseur de la couche produisant ce potentiel est donc 



=E" 



IMN. 
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Il n'y a que des termes d'ordre o, i ou 2, mais nous avons dit 
que M ne pouvait être d'ordre o ou 1 ; il reste donc n==2 et 
p=o ou/^ = 2. Siy^ = o, la surface considérée est de révolution, 
c'est un ellipsoïde peu différent du premier; sip = 2, on a un ellip- 
soïde de Jacobi. 

Variations de w^ avec l' aplatissement, — Si on fait varier 
p- de a^ à + 00 , l'aplatissement des ellipsoïdes considérés varie 
de o à I . 

Lorsque l'aplatissement est nul to^ = o ; l'aplatissement crois- 
sant, O)'- va aussi en croissant; mais il peut avoir des maximums et 
des minimums; lorsque l'aplatissement augmente, le moment 
d'inertie principal augmente aussi. 

On a vu que 



0)^ 



TU, = H=W + — J 

2 

ne pouvait pas dépasser une certaine limite. On a de plus 



0)2 



dll = rfW + ayJdoy H dJ: 

2 

dans une position d'équilibre on a 

rfW+— rfJ=o; 
2 

il en résulte que pour une position d'équilibre on a aussi 

rfH = (o Jdtù, 

Lorsque l'aplatissement augmente, J augmente indéfiniment ; 
donc (O doit diminuer et tendre vers zéro ; w^ part de zéro, revient 
h zéro, donc a un maximum dans l'intervalle. 

Remarque. — Faisons varier p^ de 00 a a^. On obtient des 
ellipsoïdes de plus en plus aplatis : on en rencontrera pour les- 
quels on aura 

R,S, R,S, 

3 2n-\-i ^ 
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je dis que la première figure d'équilibre que Ton rencontrera 

R. 
est celle d'un ellipsoïde de Jacobi. En effet, si on suppose que -^ 

est toujours (Croissant dans l'intervalle considéré, on voit que 

— Jr-^ est touiours positii. 

5 2/e+i ^ ^ 

Pour p^= » on a 

R,S, R,S, R/S, 



>— ^^> 



ô 2/1+1 

si p décroit, on finira par avoir 

RjSt R4S4 R,s, 

3 ~ 5 ^ 2AI+I ' 

p continuant à décroître, on aura alors 

R^S, Ri Si 



3 2/2+ I ' 

mais on voit que l'on rencontre d'abord un ellipsoïde de Jacobi. 

Figures dérivées de l'ellipsoïde de Jacobi. — 11 reste à étu- 
dier les figures d'équilibre qui différent peu de l'ellipsoïde de 
Jacobi. Pour cet ellipsoïde on a 

R.S. _ R,S, 
3 ~ 5 ' 

pour une figure peu différente, on aura 

R.S. ^ R.Sj 

3 2rt + I ' 

la hauteur de la surface au-dessus de l'ellipsoïde étant 

On a donc à considérer les deux équations 

R.S, ^ R,S,, ^ R.S, 
3 5 2« + 1 ' 
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On a trois inconnues p- — a^, p- — i'^, z^ — c^ ; et on a pour le 
volume qui est donne 

V=-^V/(?^-a^)(?^-i')(?*-0- 

Pour que l'équation 

R.S, ^ R,S, 

ait une racine, il ne faut pas que R,- soit divisible parRj==y p^ — a^; 
pour que l'équation 

R,S, ^ R.S, 

5 2/1 + 1 

R . 

ait une racine il faudra que T=r-i- ne soit pas constamment croissant, 

R, a donc comme seules formes possibles les formes suivantes, où 
P désigne un polynôme en p'^. 

p, 

Pxy/ r — ^S 

P X y/ p^ — 6% 

Voyons si ces quatre formes sont possibles et rappelons que 

l'on a R, = v/(?'-^')(?' — ^•'). 

i" Type, Rj est un polynôme ; si a est la plus grande racine 
de R, on pourra poser 

R, = (?^-«) n,.. 

a est compris entre a^ et c^ et on a 

K, V p-^ - z.^ • V 7^31^ " 

Oj- croit constamment, 

i / wi_ _ croit constamment, 

V P" 6'" 

/■—> 

\/ -! croît si a est inférieur à Ir. décroît dans le cas con- 

V p^— i^ ' 
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traire ; îl ne faut pas que \ / -^ soit décroissant : donc %<b' 

V S" — b' 

et toutes les racines de R,- sont par suite comprises entre b' et c". 
Parmi les fonctions de Lamé d'ordre /i, il y en a une et une 
seule ayant toutes ses racines comprises entre b' et c' ; c'est celle- 
là que Ton doit prendre pour R,. 

2* Tt/pe, Prenons R, = y z' — b' P et appelons a la plus 
grande racine du polynôme P ; on a 






Les trois facteurs sont constamment croissants, donc toutes les 
fonctions de ce type sont à rejeter. 
5* Type. Posons 

n étant un polynôme dont les racines sont inférieures à a; nous 
avons 



Ri /-; r A I 2^ — a „ 



Pour que Téquation puisse avoir une solution, il est néces- 
saire que a< b' . Donc toutes les racines de Psont comprises entre 
b"- et 6•^ 

4® Type, On volt de même que les fonctions du quatrième type 

sont telles que -~ soit un polynôme sans racines supérieures à 

a^, donc II n'y a pas à considérer de fonctions du 4** type. 

Donc pour qu'une figure peu différente d'un ellipsoïde de Jacobi 
soit une figure d'équilibre, Il est nécessaire que R, soit une des 
fonctions que nous avons appelée Ro.„. Je dis que cette condition 
nécessaire est aussi suffisante : en effet, supposons que b^ varie 
de c^ à o}\ à chaque valeur de i" correspond une valeur de p'^ telle 
que 

R,S, R.S, 

1 1 » \ , 

3 ~ 5 ' 
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considérons alors la fonction, 

F== ^*^* ^»^» . 

5 2/1 -f- 1 ' 

R. V' 

pour i"^ = c^, -^ est constamment croissant, donc F est positif ; 

si P=a^^ la valeur correspondante de p^ tend vers a^, alors 

Rv=o, 

R . ^ o puisque R, ne contient ni yp'-^ — a^ ni y P~^» ^^ 4"® toutes 
ses racines sont comprises entre b' et c^ : donc 

F<o 

et il y a certainement un système de racines pour les équations 
considérées. 

Nous voyons donc qu'il y a des figures d'équilibre voisines des 
ellipsoïdes de Jacobi. 

Remarque. — Quelle est la première figure d'équilibre que 
Ton rencontre ? c'est-à-dire, pour quelle valeur de n la valeur 
correspondante dep est-elle maximum? 

Considérons deux fonctions Ro,„i et Ry „, 

je dis que 



Ko.n Ri 



est constamment croissant; nous allons montrer que la dérivée, 
prise par rapport a w, est constamment positive ; cette dérivée 
est à un facteur positif près 

R'„R.-RaR'., 

si R, n'est pas divisible par y/p^ — a% R/ l'est; et réciproquement. 

Mais ici R^ et R,^. ne sont divisibles ni pary p^ — a^ ni par y p^ — b^; 
il en résulte que 

r;r,-r,r', 

s'annule pour p^ = a^. 
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SI on exprime R, et R;^ en fonction de Targument w, on a 

R, = F{pu) , 
R:.=;y« F [/,(«)], 

et on voit que la dérivée n'a pas de racine supérieure à a^. En 
effet la dérivée de cette fonction (prise par rapport à u) est 

(ri;'r,-r,r;') = (?*-?.)R.R* 

comme on Ta déjà vu ; R, et R;^ ne peuvent s'annuler, ©,• — «pj^ ne 
peut s'annuler deux fois, puisque c'est un polynôme du premier 
degré en p^. La dérivée ne peut donc s'annuler plus de deux fois 
pour a^ ^ p^, c'est-à-dire pour o < m ^^(k\. 

La dérivée qui s'annule pour m = o, u==(i^^^ ne peut donc s'an- 

nuler dans l'intervalle; ~- est donc constamment croissant lorsque 

Rt 

p^ varie de a^ a -f- oo . Il en résulte que l'on a aussi 



in-\- 1 2/71+1 ' 

lorsque p^ varie de oo à a^^ on rencontre donc d'abord une 
valeur satisfaisant a l'équation 



5 2/1 + I ' 

puis ensuite la valeur de p^ satisfaisant h l'équation 

R,S, ^ R,S, 

5 2/1 + 1 

La première des valeurs de n que l'on trouve est donc /i = i , 
puis /i 1= 2 : on a déjà examiné ces valeurs ; il faut maintenant 
examiner la surface voisine de l'ellipsoïde de Jacobi et corres- 
pondant a Rç 3 ; nous allons étudier d'une façon générale la 
surface correspondant à R^.n- 

L'épaisseur de la couche sera 
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cette épaisseur sera nulle aux points où Ton a 



i6i 



ou bien 



M,.N,=o, 



R,M,N,=o. 



Ceci est décomposable en facteurs quadratiques de la forme 



X 



a — a' 



y 



a?—b' 



_2 



a 



c 






I 

\ 



V^N 





t^'ff- 19- 



les a sont d'ailleurs compris entre P et c^ ; les surfaces corres- 
pondantes sont donc des hyperboloïdes à deux nappes coupant 
Fellipsoïde suivant des 
lignes de courbures 
d'un seul système. 

D'un côté d^une de 
ces lignes Ç sera positif; 
et de l'autre côté ^ sera' 
négatif. On peut, avec 
ces données construire 
la figure dans l'hypo- 
thèse de Al = 3. On re- 
présente (fig. 19) en 
ponctué le contour apparent de l'ellipsoïde, en traits pleins l'in- 
tersection de la surface avec l'ellipsoïde ainsi que le contour 
apparent extérieur à l'ellipsoïde, en pointillé le contour apparent 
de la surface intérieure à l'ellipsoïde. L'une des lignes de cour- 
bure est d'ailleurs l'ellipse principale, comme on le voit en pe- 
n) arquant que 

On construira avec les mômes conventions la figure d'équilibre 
correspondant h R^, ,^ (fig. 20). 

Je remarque qu'on ne peut prendre pour les fonctions R, de 
fonction du second ordre, car il faudrait prendre pour la fonction 
Rq2 un polynôme de la forme p" — /r, ce qui donnerait un autre 
ellipsoïde et nous avons vu qu'il ne pouvait y avoir qu'un ellip- 
soïde de Jacobi pour une vitesse de rotation donnée. 

PoiNCARÉ. — Figures d'équilibre. 11 
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En résumé, si on fait varier b^ de c^ a a*, on obtient des ellip- 
soïdes a chacun desquels correspond une valeur cj et une seule ; 
iù ne peut avoir de maximum ou de minimum, car s'il y en avait 

un Wj, à une valeur voî- 





'"-----.- 





Fi g. ao. 



sine de (o^ correspon- 
draient deux figures 
ellipsoïdales ce qui n'est 
pas. Donc <*> varie tou- 
jours dans le même 
sens. 

Supposons que 6* et 
c^ tendent vers zéro, la 
forme de l'ellipsoïde 
tend alors vers une ai- 
guille allongée; (o a d'ailleurs la même valeur pour deux ellip- 
soïdes semblables. Il suffit donc de voir ce qui se passe dans le 
cas où é et c restant fixes, a augmenterait indéfiniment. On 
aurait alors un cylindre elliptique, la force résultante de l'at- 
traction du cylindre sur lui-même et de la force centrifuge 
doit être normale au cylindre, c'est-à-dire parallèle au plan yOz. 
Considérons un point d'une génératrice passant par l'axe de 
révolution : par raison de symétrie, l'attraction doit être parallèle 
au plan yOz^ il en résulte que la force centrifuge qui est per- 
pendiculaire a Oz ne peut être que nulle, et par suite <*> tend 
vers zéro lorsque a augmente indéfiniment. 

On peut d'ailleurs remarquer que le cylindre ne peut être 
que de révolution pour que l'équilibre soit possible. 



Stabilité des figures trouvées 



Représentation graphique des résultats précédents, — 

Prenons deux axes de coordonnées (fig. 21) et portons en abs- 

R? R^ 

cisses les valeurs de -:^|- et en ordonnées -57-; ces quantités sont 

les rapports de l'axe moyen et du grand axe au petit. 

On a un premier point A correspondant au cas où ces rapports 
sont égaux à i, la figure correspondante est une sphère. La partie 



^/ 


i II' /^ 
\ \ Vf 

1 / ^^"""""^--t^^c 




/ • 





JC 



Fig. 21. 
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de la bissectrice de l'angle xOy partant du point A correspond 
aux ellipsoïdes de révolution. Les ellipsoïdes de Jacobi sont 
représentés par une courbe CBC symétrique par rapport à cette 
bissectrice et la coupant en 
un point B. 

Pour les ellipsoïdes de 
Mac Laurin, w* va en crois- 
sant quand on part en crois- 
sant du point A, jusqu'à ce 
que l'on arrive en un point D, 
puis il décroit quand on s'é- 
loigne du point D jusqu'à 
l'infini où il arrive à la 
valeur zéro. Sur la courbe 
CBC deux points symétri- 
ques par rapport à AH re- 
présentent le même ellip- 
soïde ayant tourné de 90°; w^ décroît quand on part de B 
pour aller vers C ou vers C; les points C et C s'éloignant indé- 
finiment correspondent au cas du cylindre de révolution, la 
courbe CBC est asymptote aux droites x =. i , y = i. 

Nous marquons sur AH les points E, F, correspondant aux 
figures peu difierentes de l'ellipsoïde de révolution, de même sur 
CBC, nous remarquons les points GKLM, G'K'L'M', correspon- 
dant aux valeurs Rq„ donnant les figures peu différentes d'un 
ellipsoïde de Jacobi, que nous avons rencontrées précédemment. 

11 faut remarquer que les paramètres que nous avons pris ne 
correspondent pas aux vraies données du problème. Les para- 
mètres sont le volume T, la rotation 'o et le moment principal 
d'inertie J. Si on pose [jl=(o^J, jx et T sont des données du 
problème, puisque (o^J doit être constant; [jl est du cinquième 
ordre par rapport aux longueurs, T est du troisième ; si on 
pose 

iTX 5 



i» 



M sera une donnée du problème indépendante de l'unité de lon- 
gueur choisie. 



y /M^JL 



■'? 
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Lîouville a cherché à quelles valeurs de M correspondaient les 
divers points des lignes AH et CBC. Il a montré que lorsqu'on 
se déplace sur AH, M croît lorsqu'on va de A à D et décroît 
si on va de D vers H. Sur la ligne CBC, M décroît lorsqu'on 
s'éloigne du point B soit vers C, soit vers C. 

Donc si tù est très grand, il n'y a pas de figure d'équilibre 
possible; si (o est compris entre (Og et (o^, il y a deux figures pos- 
sibles qui sont des ellipsoïdes de révolution. 

Si (0 est inférieur à Wg, il y a quatre figures d'équilibre pos- 
sible, dont deux sont des ellipsoïdes de Mac Laurin, et deux des 
ellipsoïdes de Jacobi ; ces deux ellipsoïdes de Jacobi sont d'ailleurs 
égaux comme on l'a déjà remarqué. 

Courbes d'équilibre. — Considérons un système dépen- 
dant de n variables x^x., ... Xn, soumis à un système de forces 
dépendant d'un potentiel F [x^x^.., Xny a); la condition néces- 
saire et suffisante de l'équilibrej^est que F soit maximum par 
rapport aux x, X étant donné ; on a ici 

F=W--^J = W- ^' 



2J 

La condition d'équilibre est que 

ÔF ÔF ÔF 



o; 



pour que l'équilibre soit stable, il faut de plus que la forme 

soit une forme définie négative ; on peut la mettre sous forme 
d'une somme de n carrés : 



<E) 



= V ^k'4- 



Ces coefficients olj^ sont appelés coefficients de stabilité ; si l'un 
de ces coefficients s'annule, le discriminant de la forme 4> s'an- 
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noie, puisque cette forme est alors réductible à une somme de 
moins de n carrés ; ce déterminant est, en posapt 



A = 



hXj^X^ 


"!>. 


«1,1 


«1,2 • • • 


• «l,n 


«2,1 

• 
• 


«2.2 


«2,n 


• 
• 

«n,î 


«n,2 


^n,n 



Si deux coefficients aj^ sont nuls, tous les mineurs du premier 
ordre de A sont nuls ; nous supposerons qu'il n'en est pas ainsi. 
Si A est différent de zéro, on peut résoudre les équations 






ôF 



ôF 



' 0^2 ' ' ^Xn 



et on aura 



^i = W)' 



ôF 



Différentions la relation -^ — = o par rapport à X, on aura 



s 






^x^^Xi d\ 



ôT 



Ù^Tjd). 



ce qui peut s'écrire 



I 



a 



1. 



dXi 



d^F 



d^TjdX * 



dXi 
Si A est différent de zéro, on peut trouver les valeurs de y^ , 

puisque l'on a n équations linéaires à n inconnues et les x seront 
des fonctions uniformes de )v, dans le voisinage des valeurs con- 
sidérées. 

Supposons A nul, sans que tous les mineurs du premier ordre 
le soient, supposons que l'on ait par exemple 



ÔA 



ôa 



^ o. 



1.» 
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Nous pouvons alors résoudre les n — i dernîères équations 
linéaires par rapport à 

dx^ dx^ dx^ 



ce qui montre que x^^ x^^ ... x^ seront fonctions uniformes de x^ 
et de X. 

Substituons dans la fonction F les valeurs de x^^x^^ ... a:,,, en 
fonction de x^^ et \ tirées des équations 

ôF ôF 



o,...,-7^ — = o : 



on aura une fonction ^ {x^^ X), et on a 

^__ôF^ dF_j^ ôF dxn 

dx^ dx^ djTj dx^ "' ^Xn dx^ 



mais par hypothèse 



ôF dF 



=0,...,-^ — = o 



d^2 ^Xn 



quand on y remplace x^x^,.,Xn en fonction de x^ et de X ; on a 
donc 

d^ _ dF 



dx^ ôjc^ 



La condition d'équilibre est alors que Ton ait 

4 _ 
dx^ 

C'est une relation entre x^ et X qui définit une courbe, x^ sera 

une fonction uniforme de X si — r-^ est différent de zéro ; la 

dx\ 

. d% 
condition A = o équivaut donc à \ ' =o. 

ttXy^ 

Supposons que cette condition soit remplie sans que Ton ait 



même temps --7 — ^^- = 0, la courbe a alors une tangente ver- 
ticale (fig. 22). 
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7 
3^ 



Quand X traverse la valeur de l'abscisse de cette tangente, il y a 
pour x^ deux valeurs réelles qui viennent se confondre, puis 
deviennent imaginaires; de même ^* 
pour x^j ... x^. 

01 — ^ — ^Y" = o, le point corres- 
pondant de la courbe est un point 
double, à tangentes distinctes ou 
confondues ; Tune des tangentes 
peut être verticale ou non, on a 
deux dispositions principales : 
celle qui correspond k deux tan- 
gentes distinctes et non verticales (fig. 23) et celle qui 
correspond à une tangente verticale et une non verticale 
(fig. 24). 




> 



Fig. aa. 



y 




JC, 




Fig. a3 



Fig. 24 



On a deux cas à distinguer : dans le premier cas, deux valeurs 
réelles de x se confondent et redeviennent réelles quand on tra- 
verse la valeur correspondant au point double. 

Dans le second cas, deux groupes de valeurs réelles se confon- 
dent et deviennent ensuite imaginaires lorque X traverse la valeur 
critique dans un sens convenable, une autre valeur de x est 
réelle avant et après la valeur critique. Les autres dispositions de 
figures peuvent se ramener à celles-là, du moins au point de vue 
qui nous occupe. 

Échange des stabilités. — H y a équilibre stable, avons- 
nous dit, lorsque F est maximum, x^x^ ... x„ étant variables. Il 
faut a fortiori que F soit maximum si x^ étant donné on fait 
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varier x^^ jCg, ...,^n d'une façon arbitraire, il faudra donc que 

ÔF ÔF ôF 



ÔJTj ÔJTg '" ^Xn 



O. 



On tirera x^^x^^ -"i^n de ces équations dont le déterminant 
fonctionnel n'est pas nul par hypothèse, on portera ces valeurs 
dans F {x^x,^ .•.^n)^ ^t on aura une fonction qui devra être maxi- 
mum pour qu'il y ait équilibre. Donc ^ (x^, X) doit être maximum ; 

il y aura donc équilibre si , = o, mais il n'y aura équilibre 
stable que si l'on a en même temps > a < o. 

Un arc de la courbe -^7-^ = partage la région du plan où 



, 4 



, rf-î- 



il se trouve en deux parties, l'une où y^ > o, l'autre où ,^ <o ; 



SI c est un arc sans sin 



dx. 



dx. 



ingularité, , \ ne sera pas nul et sera 

négatif si la région où ^ est positif est située au-dessous de 

l'arc de courbe; et au contraire sera positif, si est positif au- 

ax^^ 

dessus de l'arc de courbe considéré. Le premier cas correspond à 
une position d'équilibre instable. 

En couvrant de hachures la région où .^ est positif, on voit 

que la figure aS correspond au cas d'un équilibre stable et la 
figure 26 au cas d'un équilibre instable. 



J/ 








y\ 




oa 



Fig. 25. 



Fig. a6. 



' Dans le cas où l'arc BAC a une tangente verticale en A, deux 
cas peuvent se présenter, suivant que c'est la partie convexe 
ou la partie concave qui est hachurée. 



ÉCHANGE DES STABILITÉS ■ 



T6.J 



Dans la figure a^ : 

l'arc BA correspond à des positions d'équilibre instable, 
l'arc AC » » stable; 

dans la figure 28 : 

l'arc BA correspond à des positions d'équilibre stable 
l'arc AC » » instable. 





Fig:. : 



Ce qu'il y a d'essentiel, c'est qu'en suivant l'arc et en passant 
par le point A, les deux systèmes de figures d'équilibre échangent 
leur stabilité. 

Dans le cas d'un point double sans tangente verticale, deux 
cas peuvent aussi se présenter; dans les figures 29 et 3o, les 




Fiy. ag. 



Fig. 3o. 



arcs BA et AC' correspondent à des figures d'équilibre instable, 
les arcs B'A et AC correspondent ii des figures d'équilibre stable. 
Remarquons encore qu'en suivant l'arc BAC ou l'arc B'AC les 
figures d'équilibre, de stables deviennent instables ou récipro- 
quement. 

Enfin dans le cas d'un point double avec une tangente verti- 
cale, on arriverait à la même conclusion, les figures d'équilibre 
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obtient en suivant un arc quelconque échangeant leur 
stabilité quand on passe par un point où ■j-~= o (fig, 3i et 32). 




Fig. 3a. 



Stabilité de l'équilibre des Ûg^ree trouvées. — Appliquons 
ceci à notre problème et considérons la figure 21; ce que nous 
avons appelé [i représentera ici ce que nous avons appelé X et 
<» correspondra à x^. 

La fonction F, qu'il faudra rendre maximum est 



W- 



aJ 



w" augmente lorsqu'on se déplace de A sur D, et décroît de 
D à H; sur AH, M va en croissant, et croît aussi de B à C ou 
à C. 

Il s'agit maintenant de traiter la question de stabilité. 

Partons de A, où b) == o et où on a une sphère, nous savons 
que dans ce cas l'équilibre est stable : donc tous les coefficients 
de stabilité sont négatifs. Ils cesseront d'être tous négatifs quand 
on rencontrera une figure de bifurcation. 

La première bifurcation que l'on rencontrera en suivant AH 
sera au point B, où l'ellipsoïde de Mac Laurin diffère peu d'un 
ellipsoïde de Jacobi. Le coefficient de stabilité s'annulera donc 
en même temps que 

R.S. R,S, 
3 5 * 

Après cela il ne s'annulera plus ; quand on arrivera en E, un 
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autre coefBcient s'annulera en changeant de signe, mais le 
premier coefficient restera positif. Donc de A a B nous avons un 
équilibre stable, au delà de B il est instable. 

Sur la branche BC, on a au contraire des ellipsoïdes qui sont 
stables jusqu'en G ; en G Tun des coefficients s'annule et au delà 
reste positif, donc à partir de G sur la ligne BC on aura des 
positions d'équilibre instable, mais à partir de G commenceront 
de nouvelles figures d'équilibre et ce sont celles dont il faudra 
étudier la stabilité ; il y a une difficulté, car, pour ces figures, 
M peut croître ou décroître. 

Il est aisé de voir que l'on a en tous cas un maximum ou un 
minimum ; on a en effet pour ces figures, si elles correspondent 
à des formes peu différentes de celle que l'on a en G 

en G même, e = o. Si on change e de signe, la figure a une 
forme symétrique de la première, M garde donc sa valeur; il en 
résulte qu'en G la quantité M passe par un maximum ou un 
minimum. Si M<Mo on a après une seule figure d'équilibre; si, 
au contraire, M > M^ on aurait trois figures d'équilibre dont une est 
stable. On n'a pas encore fait le calcul complet, bien qu'il semble 
plus probable que l'on a un minimum. 

Conclusion, — Quoi qu'il en soit voyons les deux hypothèses : 
supposons une masse fluide homogène animée d'un mouvement 
de rotation et soumise à un refroidissement assez lent pour 
qu'il ait le temps de se propager dans la masse ; supposons de 
plus que le mouvement de rotation soit le même dans toute la 
masse, et cela en supposant le frottement suffisamment grand. 

La vitesse du mouvement de rotation n'est pas constante, ce 
qui est constant c'est u = a>'^J, en vertu du théorème des aires ; le 
volume T doit diminuer puisque la masse se refroidit, donc M 
augmente. On a des figures d'équilibre stable jusqu'en B ; après 
cela, ce qui devient stable ce sont les ellipsoïdes de Jacobi et cela 
jusqu'en G ; le refroidissement continuant, on ne peut plus avoir 
d'ellipsoïdes de Jacobi, puisqu'ils sont instables. Si, pour l'autre 
série, M^ est un maximum, cette seconde série ne donnera aucune 
figure, et le corps se dissipera dans l'espace; si, au contraire, 
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Mo est un minimum, on aura une sorte d'ovale étranglée que 
l'étranglement divise en deux parties inégales. 

On n'a pas fait le calcul jusqu'au bout pour voir ce qui 
arriverait ensuite : peut-être l'étranglement s'accentue ra-t-il et 
la masse pourra se diviser en deux parties. Remarquons que ceci 
ne peut s'appliquer directement a l'hypothèse de Laplace, car 
la nébuleuse considérée par Laplace n'est pas homogène, elle 
doit au contraire être fortement condensée au centre. 

Quoi qu'il en soit, on ne sait pas jusqu'à présent si en G corres- 
pond pour M un maximum ou un minimum. 
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ANNEAU DE SATURNE 



Etude des ûgures annulaires d'équilibre. — On peut faire 
trois hypothèses sur l'anneau de Saturne : on peut admettre qu'il 
est solide, qu'il est fluide ou qu'il est composé d'une masse de 
solides indépendants ; la troisième hypothèse connue sous le nom 
d'hypothèse de Cassini parait être l'hypothèse véritable. 

HYPOTHISSE DE l'aXNEAU SOLIDE 

Historique, — Laplace est le premier qui ait examiné cette 
hypothèse : et il a remarqué que si l'anneau était homogène son 
mouvement serait instable ; pour peu que son centre s'écartât 
du centre de gravité de la planète, Tanneau se précipiterait sur 
celle-ci ; l'attraction de Saturne ne tendrait qu'à accentuer l'ano- 
malie. Si donc l'anneau est solide il doit présenter des irrégula- 
rités ; Laplace ne chercha pas à déterminer l'ordre de grandeur 
de ces irrégularités ; Maxwell a fait le calcul et a reconnu 
qu'elles doivent être considérables ; ces irrégularités sont inad- 
missibles. 

Equations du mouvement. — Parmi les forces qui sont appli- 
quées a l'anneau il n'y a pas a s'occuper de l'attraction de l'anneau 
sur lui-même parce que l'anneau étant solide ses diverses parties 
n'éprouvent pas de déplacement relatif; il n'y a donc à s'occuper 
que de l'attraction de Saturne et de celle des autres corps 
célestes ; nous négligerons cette dernière action et nous allons 
nous occuper seulement de l'attraction de Saturne. 

Nous allons supposer que l'anneau a son centre de figure au 
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centre O de la planète : l'observation montre qu'il en est ainsi; 
le centre de gravité G de l'anneau est supposé différent de O. 
Un point M de l'anneau, puisque l'on ne parle pas de l'attrac- 
tion de l'anneau sur lui-même, est soumis à l'action de l'attraction 
de Saturne ; l'ensemble des forces d'attraction fait donc équi- 
libre à la force d'inertie. Un point M de l'anneau est soumis à 
l'attraction de Saturne et à la force centrifuge : si l'anneau con- 
serve son centre de figure au même point cette force centrifuge 
et cette attraction sont constantes : M tournera alors avec une 
vitesse égale à celle d'un satellite placé à la même distance ; 
c'est la condition nécessaire et suffisante de l'équilibre. Si 
l''anneau est homogène, cette condition n'est pas nécessaire, car 
la vitesse de rotation peut être quelconque : chacun des points de 
l'anneau serait soumis à la même force, et l'équilibre aurait lieu 
quelle que soit la vitesse de rotation. Dans le cas général le 
centre de gravité G décrira donc d'un mouvement uniforme une 
circonférence autour de Saturne. 

Prenons comme unité de longueur le rayon de l'anneau, 
comme unité de masse la masse de l'anneau ; son moment d'inertie 
par rapport au centre de figure sera égal à i . Il reste à déter- 
miner Tunité de temps : on peut la prendre de façon que, en 
désignant par M la masse de Saturne et par /*le coefficient de 
la gravitation, on ait /'M= i. 

Etudions le mouvement : soient x e\ y les coordonnées du cen- 
tre de gravité, nous avons 



coss.rfs, 




en appelant p la densité de l'élément linéaire de l'arc ds de 
Tanneau supposé étroit. 



est la masse de l'anneau. 
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Les équations du mouvement du centre de gravité sont 




ods = 



p cos s ds 




d'où 



de 

d^y 
di' 



==07, 



y- 



Les intégrales de ces équations sont, en prenant une origine, 
des temps convenable, 

r ^ = a cos ^, 
/ y = P sin t. 

Le point G décrit donc une ellipse dans le cas général ; si 
c'est une circonférence, elle est décrite avec une vitesse angulaire 
égale à Tunité : c'est la vitesse d'un satellite placé à l'unité de 
distance et décrivant un cercle et dans ce cas, a = p. 

Stabilité du mouvement, — Il reste a voir si ce mouvement 
est stable. Supposant négligeable, par rapport h la masse de 
Saturne, la masse de 
l'anneau, nous pou- 
vons supposer Sa- 
turne fixe. 

Par son centre O, 
faisons passer deux 
axes fixes OX et OY 
(fig. 33) ; par le cen- 
tre de figure de l'an- 
neau, faisons passer 
deux axes ox^ oy in- 
variablement liés à l'anneau, leur position, et par suite celle 
de l'anneau, sera déterminée, si l'on se donne les coordonnées 
a; et y de par rapport à xoy^ et l'angle de l'axe ox avec OX. 

Soit G le centre de gravité de l'anneau, qu'on peut supposer 




Fig. 33. 
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sur ojT ; on a ôG = a, c'est une donnée du problème ; les coor- 
données X et Y du centre de gravité sont 

X = -4-y sinO + (a — jr) cos 9, 
Y :=: — y cos + (a — x) sin 9. 

La force vive de Tanneau est égale à la forcé vive de trans- 

lation (X'^ + Y'^) augmentée de la force vive de rotation W^. 

^ . . . ^" 

Le moment d'inertie par rapport au point o est donné par la 

relation I = J + a^. La force vive est donc une fonction de x, y, 9 

et de leurs dérivées et on a 

X/2 _^ Y'^ ^^/2 _^ y, 2_^ Q/2 [-^2 4_ (^ _ ^)2] _ ^0' [x' y + [a —x)y']. 

L'énergie potentielle est 

U =— /MV, 

V désignant le potentiel de l'attraction de l'anneau au centre de 
gravité, et comme on a 

on a 

U=— V. 

Supposons que la densité de l'anneau soit une fonction déve- 
loppable en série de Fourier. Si di^ est l'angle sous lequel on 
voit du centre de l'anneau un élément linéaire, on aura, en appe- 
lant p la densité au point déterminé par l'angle ^, 

r . , 2a I 28 ... 1 

P==?o| ' + ^aGos^-\ — ^ cos 24^ + -^ sin 2y + ... I; 

si \ et 7| sont les coordonnées d'un point de l'anneau par rap- 
port aux axes mobiles, on a 
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prjrfs=o, 



p(~$^ + ri^)<fc = 




a 



On trouve encore 



2ic 




p^Tj rfs = I psîn2^rf^ = -^. 



Le potentiel en G de l'attraction de Saturne sur l'anneau 
est, comme on le voit en effectuant le calcul, 



/pds 



Les équations de Lagrange donnent 

rf /dT\ dT ÔU 



\àf/) 



G 



O. 



dt \ àf/ ] ^q àq 

mais comme U ne dépend pas de q' on a 

d d(T — U) d(T — U) 
dt ^q' àq 

On peut écrire 

T-u= Ae--+^BQ'+c +v, 

où l'on a posé pour abréger 

C = x" + t/"y 

PoiNCARÉ. — Figures d'équilibre. i'^ 
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remarquons que T — U ne dépend pas de 9 ; la première des 
équations de Lagrange donne donc : 



r a(T-u) i 

L ô9' J ' 



_£r d(T — U) 
dt 



et, comme U ne dépend pas de 6, il reste 

dH _ 

p étant une constante. Il est d'ailleurs facile de calculer cette 
constante, étant données les conditions initiales ; si on déve- 
loppe l'équation on a 

A9'+B=/;. 

Supposons j;,:=:o, y, = o, et par suite x\ et y', très petits, on 
trouve 

B, = o, 

et par suite p = ^'^\ donc c'est la vitesse de l'anneau, supposé en 
équilibre : on a, comme on l'a démontré plus haut, ^'q= i ; 
donc p=i : l'équation obtenue est Téquation des aires. 

Transformons maintenant les équations de Lagrange : posons 

Les équations de Lagrange conserveront la même forme, si on 
y remplace T — U par H ; car, si q est l'une des variables, on a 

ÔH d(T — U) ôT W W 



dq dq * dO' àq ' dq ' 

et, comme on a 

_ ôT . 



on conclut : 



ÔH _ d(T— U) 
d^^ dq 
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De même on voit que Ton a 

ÔII _ ()(T — U) 
^(f ~ àr/ 

On peut faire disparaître 6' de ces équations, on n'a en effet qu'à 
remplacer 8' par 

A ' 

et comme, ni T ni U ne contiennent 9, II ne sera plus fonction 
que des deux variables x et y, et de leurs dérivées. On trouve 



o A o I 



2A 

Supposons maintenant que le centre de l'anneau soit très près 
de Saturne, c'est-à-dire supposons :r et y petits, et voyons si dans 
le mouvement ultérieur a: et y resteront petits. 

Les termes de degré zéro qui entrent dans II peuvent être 
négligés, car ils n'interviennent pas dans les équations de 
Lagrange. Les termes du premier degré en x et y doivent dispa- 
raître d'eux-mêmes : en effet, dans les équations de Lagrange ils 
donneraient des termes constants^ et ces équations doivent 
admettre la solution x = o^ !/ = o. Quant aux termes du premier 
degré en a/ et y^ ils disparaissent dans les équations de Lagrange. 
Il ne reste donc que les termes du second degré, que nous 
appellerons W, le reste étant négligeable. 

Développons H, en ne conservant que les termes du second 
degré ; on trouve pour W le développement suivant : 

Ex* H- aCxy + Fy " + — + ' ~^ y'^ ^,r'y +{i — sa') xy' 



en posant 



E=_-2a+-, 

4^4 
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Les équations de Lagrange, comme on l'a remarqué, peuvent 
se réduire à 

d ôW _ ôW 
dt àx' dx ' 

rf dW ôW 
dt dy' ôy ' 

c'est-à-dire 

x" — 2 (i — a^)y=2E^ + 2 Cy, 
(, _«2)y/+ 2 (i — a')y= 2C^+2Fy. 

Ce sont deux équations linéaires a coefficients constants ; pour 
intégrer on pose 

On aura quatre valeurs possibles de w, et l'équation générale 
du mouvement sera 

les l et les /ti étant des constantes quelconques déterminées parles 
conditions initiales. Pour que le mouvement soit stable, quelles 
que soient ces conditions, il faut que les valeurs de w soient à 
partie réelle négative ou nulle. 

Pour déterminer w, on remplace a; et ^ dans les équations 
différentielles et on a 

— lo^jc — 2(1 — a^)ia>2/^= 2 Ea7,+ 2 C^, 
— (i — a^)a>^y^+2(i — a^)((i)jr = 2Cjr+2F^, 

on peut éliminer ^ et y entre ces deux équations, et on trouve 

L'équation en o)^ doit avoir ses racines réelles : car si elles 
étaient imaginaires, les racines carrées de ces deux racines ne 



POTENTIEL D'UNE CIRCONFERENCE HOMOGENE i8r 

seraient pas des imaginaires pures, et par suite, comme elles 
sont égales et de signes contraires, Tune serait a partie réelle 
positive. 

Or, pour un anneau homogène, oii 

a = a= p = o, 
l'équation se réduit à 

(»>* — (o^H — ^ = o; 



les racines en sont imaginaires : si a, a, ^ sont voisins de 
zéro, elles seront encore imaginaires. Il faut donc que a, a et ^ 
soient grands, il y aurait donc de grandes irrégularités dans 
Tanneau : et sa régularité apparente serait inexplicable. 

Maxwell a fait le calcul dans l'hypothèse où l'anneau serait 
homogène, sauf en un point, où il y aurait une masse supplé- 
mentaire : et il a trouvé que cette niasse devrait être au moins 

les -^ de la masse totale de Tanneau, c'est-à-dire que l'on aurait 

a>058. 

M. Radau a examiné le cas où l'épaisseur de l'anneau serait 
variable d'un point à un autre, et il a trouvé que la densité 
devrait varier de 2,7 à o,o4, ce qui ne paraît pas probable. 

De plus, la section de l'anneau étant mince, il lui faudrait 
une grande résistance pour résister à l'attraction des satellites : 
Hirn a calculé qu'il lui faudrait une résistance mille fois supé- 
rieure a celle de l'acier : donc Vanneau de Saturne n^est pas 
solide, ' . ' 
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Potentiel d'une circonférence homogène. — Cherchons le 
potentiel d'une circonférence homogène, en un point P de l'es- 
pace. Soient PA = a, PB = b les distances minima et maxima 
du point P à la circonférence. La circonférence étant donnée 
de rayon R, aet i déterminent complètement le potentiel. 



i8a 
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Soit 10 l'angle AOM (fig. 34) , on a 



PM' = a' cos^ — + i* sin* — 

2 2 




le potentiel 



Vp = 



= M 




Fig. 34. 



211 XPM 



Comme le potentiel ne dépend que de a, de b et de la masse 
de la circonférence, on peut ramener 
•le calcul de V au potentiel d'une cir- 
conférence, en un point de son plan : 
on décrit une circonférence de dia- 
mètre (<2 -|- i), et on calcule le po- 
tentiel de cette circonférence en un 
point situé à une distance a d'une des 
extrémités du diamètre (fig. 35). 

Le potentiel en F est le même, 
on a 




Fig. 35. 



V=M 
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\/a 



^COS^- ' '2. ,-.2 



2 2 



mais on a aussi 



V=M 




4 



27iv/oZ 



cos 



2tli 



Q étant le point de la circonférence dont la projection sur le 
diamètre AB est en F, ^ étant l'angle AOM ; on a donc 

? {a, b) = 'f (OA, QP) 



=?(-^'v/^)- 



POTENTIEL D'UNE CIRCONFÉRENCE HOMOGENE iSJ 

Si donc on forme une série de quantités 

a b 



telles que Ton ait 

^„ _ 1 + i„ _ 1 



«n 



on aura 

V = Mf(a„,è„), 

et cela quel que soit l'entier n. On démontre facilement que a„ 
et bj^ tendent vers une même limite qu'on appelle la moyenne 
arithmético-géométrique de a et de b. 
Soit m cette limite, on a 




d^ 



271 1/ m'^ cos^ — -|- m'^ sin 



•2 r 



= M 




îp (<2, i) est une fonction homogène de a et de i, de 
degré — i , c'est-à-dire que Ton a 

'^ [la, Ib) = > ~ * '^ (a, i) ; 

on peut donc écrire 

et, remplaçant cp (a, i) par cp (a), on a 



?(«.i)=rf(y) 



Que devient la formule si P est rapproché de la circonférence? 
a est alors une quantité très petite. On sait que Ton a, en dési 
gnant par L le logarithme népérien 

r 
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p désignant la distance du point P au point P^ le plus voisin du 
cercle, e étant une quantité tendant vers zéro avec p, jjl^ et K étant 
des constantes ne dépendant pas de la position du point P par 
rapport à P^ non plus que de la position du point P^ sur le cercle. 
En reprenant les notations précédentes on a 

, K 





' ^J*« '^ a ' 


et comme 






M _ 2irjA,R, 


on a 






M kR o 



b diffère peu de 2R, à l'ordre d'approximation adoptée ; on a 
d'ailleurs 

Cette relation permet de calculer la constante K : on a en efiet 

,K*,=,(v/S:îi±^) = ^,(^); 

si on fait i = i on a 

s/a 



?W = TT^?(TfT)- 



11 vient donc dans le cas qui nous occupe, en négligeant a 
devant l'unité 



mais on a aussi 



donc 



ç (a) = L — ; 



\ « / \ i\/a / 

LK — La = 2 LK — 2 log 2 — La ; 



d'où on conclut 



MASSE ANNULAIRE DE REVOLUTION i8 = 

LK = 2 La, 

K = 4, 



et finalement 



M 



=^■-(4) 



avec Tapproxiination donnée. 

On peut chercher une forme plus exacte pour f [a). Posons 

?(«) = -^ log— ["i + A.a + AX-l- --J 
en se servant de la relation 



on trouve facilement que Ton a 

f A, = o, B„=B, = o, 

A =-L R=— i 

4 

de façon que Ton peut prendre pour f (a) la fonction 

4 



l'H 



B,= ^La, 



7 >»8 (i) 



et cela en négligeant seulement les termes ea — r^ . 

Masse annulaire de révolution. — Soit G le centre de gra- 
vité de la section méridienne, menons Gy parallèle à Taxe de 
révolution et Gj7 dans le plan méridien, perpendiculaire a G/y 
et tel que le sens positif soit dirigé de Taxe de révolution vers 
le point G, soit l la distance du point G à cet axe de révolution. 
Un élément de surface rfd' de la section méridienne, de densité p' 
et de coordonnées x' et y' engendre un anneau dont la niasse est 
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Cherchons le potentiel de cet anneau en un point P du plan 
méridien, dont les coordonnées soient a: et y ; il faut d*abord 
calculer les quantités a et b. Si on suppose le point P du même 
côté de Taxe de révolution que l'élément rfo-, on a 

a=v /(x-xO'Hr(y-;y)S 

Supposons l'anneau mince et le point P situé près de G, de 
façon que ^, y, x\ y' soient négligeables devant l; a est une 
quantité du même ordre ; on peut prendre pour valeur approchée 
de b 

car [y — y'y est négligeable ; il faut maintenant calculer îp (a, b). 
On a 

?K*)=t?(t) 

= — !— L^*-. 

il -\-x-\-x' a 

Le potentiel en P est donc 

b a 

on a 

l + a:' l-\-x' 



il-\-x-\-x' 



comme par hypothèse x et ^' sont petits, on peut écrire approxi- 
mativement 



l + x' 



=^[■+■^1 



Il vient ensuite 



a a * 21 
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ce qai peut s'écrire, avec le de^é d'approximation donné 

on peut alors écrire aussi 

Le potentiel de Tanneau sera donc 

les autres termes étant négligeables, et les intégrales étant 
étendues à la surface de la section. Nous allons supposer Tanneau 
homogène et nous écrirons p' = i. 

Les deux premiers termes sont des potentiels logarithmiques 
qui jouent, dans Tattraction des masses planes^ le même rôle que 
le potentiel newtonien dans l'attraction des corps célestes ; les 
propriétés du potentiel logarithmique dans le plan sont anulo* 
gués à celles du potentiel newtonien dans Tespace : on démontre, 
par exemple, que le potentiel logarithmique d'un cercle homo- 
gène est le même que si toute sa masse était concentrée au centre. 

Si le point P est à la surface d'un tore, à une distance a^ du 
centre du cercle générateur, il suffit de considérer le potentiel 
logarithmique en P du cercle de rayon a^ ; sa surface est tî^iJ, 
Son potentiel est donc 



8^ 

0' 



TZûa Li 1 



«0 



la première intégrale est par suite 

8^ 



HTzal L 

^0 



La seconde est la somme de deux autres : 




a fia 
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dont la seconde se calcule immédiatement, on trouve pour sa 

valeur 

X 2 T 8^ 
-y- TtaJ L ; 

la première va se calculer plus tard. 
On tro.uve de même 





d3'=0 



puisque le centre de gravité est à l'origine des coordonnées. 
Il reste à calculer 



/x' , SI , 



nous allons employer un procédé de calcul dont nous nous 
sommes déjà servi dans une autre question : supposons que la 
densité du cercle soit 



al x"+y'^ 



le potentiel sera 



U = 




la densité ne dépendant que de la distance au centre, le potentiel 
logarithmique du cercle est le même que si toute la masse est 
concentrée au centre. Cette masse est facile a calculer, c'est : 
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«aj 



4 

Le potentiel est donc, en un point A, situé à une distance a^ du 
centre du cercle a^ > a^, 



L 



4 «1 

Supposons que le cercle C se déplace d'une quantité e parallè- 
lement à l'axe des x, de façon à 
occuper la position C^; le poten- 
tiel en A sera (fig. 36) 

u ^" 




la différence de potentiel sera donc 

OU 

ÔJ; Fig. 36. 

Nous allons étudier d'où provient cette différence : en un 
point de l'intérieur du cercle la densité était p', elle est devenue 

— Ê . 

Dans la région du cercle C extérieure au cercle C, la densité 
était nulle, elle est devenue égale à p' qui est petit puisque la 
densité p' est nulle sur la circonférence du cercle C, et l'aire de 
la région considérée est de l'ordre de e. Il y a encore à consi- 
dérer l'aire égale à la précédente, intérieure au cercle C et exté- 
rieure au cercle C, la variation de densité en sens contraire de 
la précédente est 'aussi de Tordre de e au plus, puisqu'on un 
point de cette aire la densité p', petite sur le bord du cercle, est 
devenue nulle. 

La variation de potentiel due à ces deux aires est donc négli- 
geable, et on a 

£ — = — A -^ L — da'. 
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Or 

il vient donc, en supprimant le facteur e, 



,T 8/ ,, dU 




a dx 



L'intégrale que nous voulons calculer est donc, au facteur Z près, 

^ dU dU 1 1» • . 

égale a -^ — ; or — — est la composante de 1 attraction suivant 

• I • X 

Ox ; cette attraction est , sa composante suivant Ox est — 5- . 



a^ a; 



On a donc 



x' 81 j Tza^x 

L a(T = 



2 ' 



la ^l a^ 



et, si on suppose le point situé sur la surface du tore, c'est-à-dire 
si l'on a a^ = a^^ il vient pour la valeur de l'intégrale 

Tzai 



4r 

On a donc pour le potentiel en un point de la surface du tore 
= 7.Tzal L r ira; L [- -77- izal, 

a} 
en négligeant les quantités de l'ordre de —^ • 



Section méridienne de l'anneau. — Si Ton appelle M la 
masse de Saturne, to la vitesse de rotation de l'anneau, une 
condition nécessaire de l'équilibre est que, à la surface du 

^2l>2 

tore, la fonction U -j soit une constante, U représentant 

le potentiel de l'anneau sur lui-même augmenté du potentiel de 
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Saturne, ce dernier peut être représenté, avec l'approximation 

M , 

employée, par— y— et Ton doit avoir à la surface de l'anneau : 

En posant 

W=l^ + llX', 

la distance d'un point de la surface au centre du tore est 



ce qui, avec 


l'approximation que l'on cherche, 


peut s'écrire 




p^-(/ + x)% 






d'où 


p_/+^; 




• 


enfin on a 










I I I 


X 






p l-^ X l 




La condition d'équilibre peut s'écrire 








, _ 8Z X 2 ï 8^ 1 
«0 ^ % 


5ar 

4^ 


■Kal 






Ix — 


= 0"; 



il faut donc que le coefficient de x soit nul, ce qui donne 

5Tzal iraS , 8/ M 



^ l,— 11 ^ ^^2; ^ ^^ 



41 lai 

Celte équation donne w'^, aux quantités du deuxième ordre 

près. 

^2 
Comme -j- est aussi petit, on a 

s M 



19a 
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La vitesse de rotation d'un tore engendré par un cercle, de rayon 
suffîsamment petit par rapport à sa distance de Saturne, est 
donc la même que celle d'un satellite placé à la même distance. 
S'il n'y a pas de corps central, M = o : on voit alors qu'en 
l'absence de corps central il peut y avoir une figure annulaire 

d'équilibre : a^ est alors très petit, mais le rapport tend vers 


/ 

l'infini quand a^ tend vers zéro, car L — augmente alors 

indéfiniment. 

Si on voulait pousser plus loin l'approximation, il faudrait 
calculer les termes de la forme de 



r 



Q/ 

^ a ' 



et 



x""' y" Tx""' tj"" 



d^. 



Supposons maintenant que le tore ait une section elliptique ; 
nous devons alors calculer le potentiel logarithmique d'une 
ellipse. Or le potentiel new^tonien d'une droite indéfinie en un 
point, est le potentiel logarithmique en ce point de la projection 
du point sur la droite ; il en résulte que le potentiel logarith- 
mique d'une ellipse en un point de son plan est le potentiel 
new^tonien d'un cylindre long par rapport a sa section droite, 
i en un point situé à une petite distance de sa surface. On est 
donc amené à calculer le potentiel d'un ellipsoïde très allongé 
pour lequel on supposerait c^ infini. 

Le potentiel intérieur est donc une fonction du second degré 

Y==Ax^ + By' 

dont les coefficients sont faciles à calculer (page i4o)- 
Pour le potentiel extérieur on a 

V = A,+ A,x + A,x' + Ay + .. . 
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On peut calculer A^, A^, A^ comme on Ta fait précédemment, 
mais les intégrales devront être étendues à la surface d'une 
ellipse dont l'équation est 

> •> 

a b- 

Lcs quantités du troisième ordre sont négligeables. La fonction 
de force peut de même se développer en série de la forme 

La condition d'équilibre est que l'on ait 

(A, + C.) X + (A,+C,),r^ + (A, + C,)t/ ^ C- 

à la surface de l'ellipse ; il faut donc que l'on ait 

A, + C, = o 

(A,+ C,)a^ = (A3+C,)i^ 

Laplace avait fait le calcul, et a discuté l'équation mais sans 
tenir compte des termes en A^ et A3. 

M'"® de Kowalewski a été plus loin et a fait voir que la section 
méridienne du tore n'était pas symétrique par rapport à Taxe 
des y. 

Nous étudierons plus tard la question de la stabilité ; nous 
allons d'abord discuter Thypothèse de Cassini. 
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Vraisemblance de rbypotbèse- — Cette hypothèse est la seule 
compatible avec les travaux de Maxwell sur l'anneau de Saturne. 
L'observation confirme cette hypothèse, car l'anneau est trans- 
parent et la lumière le traverse sans trace de réfraction ; il faut 
donc supposer qu'il est composé de corpuscules solides (ou 
liquides), isolées les unes des autres. Les observations spectros- 
copiques montrent d'ailleurs que la vitesse d'une molécule de 
l'anneau n'est pas la même sur le bord interne, ou sur le bord 
externe. 

PoixcARÉ. — Figures d'équilibre. i3 
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Equations du mouvement — Noas allons commencer par 
^tudi'rr un vu^ trêâ simple : nous allons considérer p satellites 
Je m^me masse ;x, équidîstants sur un cercle : la distance 

de deux «afellîtes sur ce cercle sera — ^ = \th; le rapport -7- 

n'est pas infiniment grand. 

l'n mouvement possible du svst^me est celui où chaque satel- 
lite parcourrait le cercle d'un mouvement uniforme, avec une 
vitesse u»^ déterminée par l'attraction de la planrte, augmentée 
d*un U;T%u(i où il faudrait tenir compte de la masse des autres 
satellites; et nous allons chercher s'il va un mouvement peu 
différent de celui-là. 

Soit M4 le satellite de rang /•, son rayon vecteur sera i -J- S;, 
son angle polaire sera tôt -\- '^hh -f- 7/ : dans le mouvement non 
troublé on aurait 

7f, = o ; 

si le mouvement est stable Zf^ doit rester petit. 

I.a force vive T, et la fonction de forces U, sont fonctions de Zj. 
et de or^; les équations de Lagrange donnent 






dt 

d / ôT \ 4]iï m_ 



= : 



fri posant 



T — IJ= H, 



on îi 



d «MI «MI 

y^ o. 



dt (>pV ^? 

Nous pouvons supposer que Saturne est fixe, c'est-à-dire que 
nous négligerons les perturbations apportées dans le mouvement 
de Salurne par celui des satellites ; la mise en équation du pro- 
blème est alors Indépendante du mouvement de Saturne^ 
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La force vive du système est 

et, en négligeant les infiniment petits du troisième ordre : 
La fonction de forces est 

R représentant la partie complémentaire du potentiel due h 
l'attraction des satellites les uns sur les autres ; si l'on pose : 

on a 

Aj 2sin4; I 2 4 8 " ' 

les termes non écrits étant négligeables. On peut écrire aussi 

Les termes les plus importants de R sont évidemment ceux 
pour lesquels 4* est petit, sîn à étant un infiniment petit du pre- 
mier ordre, — . ^, est un infiniment errand du second, et . ", 

sm ^ ^ sm ^ 

est un infiniment grand du troisième. 

U et T comme on le voit ne dépendent pas de ty et sont du second 

degré par rapporta p,? ,^ et o* ; — — , , ' sont du premier 



#K:gr'; par rapport ii ces quantitr?^ ; enfin —^ — ( — r— )est dupre- 
mi^îr A*iyijh u\'f:e rl^5 dérivées du second ordre. Les équations de 
La^rang^ sont donc des équations ici coetïicients constants de la 
forrnf! 

*\ • I) ^« ' « ■^ • ^» •* «^ -~^ o, 

K B ' & ' 



• t ' t 



A étant finrr constante, et B étant une Tonction linéaire des ar^ru- 
ments; mais ces équations doivent admettre la solution p = a' = o; 
il faut donc que A disparaisse de lui-même; les équations sont 
alors des équations linéaires et homogènes, qui s'intègrent par 
les exponentielles, et on a. 



r=y 



ae" 



les A étant des racines d*équations algébriques, et les a étant des 
c.r»ristantes. 

Le nionvcnieiit sera stable si p ne croit pas au delà de toute 
limite ; en mettant en évidence la partie réelle et la partie imagi- 
naire de A, on a 

A -^ \ 4- ^^-i ; 

le mouvemeiil sera stable si A^ est négatif ou nul ; mais si nous 
démontrons <[ue ré([uation en A a ses racines deux à deux égales 
et (l(^ signes contraires, nous verrons que )v,,ne peut être que nul. 
Or si nous cliangoons l en — t les équations ne changent pas, 
puis(|ue les termes du second ordre seuls interviennent dans les 
é(|ualionH de Lagrange ; il en résulte que si l'équation en A 
admet une racine, elle admettra aussi la racine de signe con- 
traire, et par suite, la condition nécessaire de la stabilité de 
ré([uilibre, est que les racines de l'équation en X soient purement 
imaginaires. La condition est d'ailleurs évidemment suffisante 
puis([ue, dans ces conditions, p s'exprime par des fonctions trigo- 
nométri(pies de la variable réelle t. 

Vous avons y,/> é([uations du second ordre, donc ). satisfait à 
une équation algébrique d'ordre 4/^; c'est une équation en A*, de 
degré 'jp dont toutes les racines sont réelles et négatives. 
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Les premières des équations de Lagrange sont 

pi'— or ( I + p,) — 2(0^'/,+ I — 9.p, = ;j. (-^ j ; 

«n peut écrire, en mettant en évidence dans'R les termes de 
degrés o, I et 2 

R^.R, + R, + R,; 

et Ton a 

OR ÙR, OR, 

OR, OR. ^ ^ . T , . 

-^ eslune constante, est une lonction linéaire. 



Les secondes équations sont 



, OR, OR., 



Dans la première équation les termes constants doivent dispa- 
raître : on a donc 

OR, 

(0- — I = — ;^ -— -i- 

4 ^ sinA ^ 



^ prend les valeurs Q, 26, 30, ,.., (/; — i) 0. Dans le cas ou [ji=o 
on aurait (o^= i. 

Les équations sont donc 

p;; — 2<o< — ((0-^ + 2) p,. == [^ -T-^, 

\ A 



■ÎA aW?/. = .'A 



' 0,, • 



Changement de variables. — Faisons le changement de 
variables, 



^ i,.->.- , 



Y — ^ t^'i* »• ^P' 



\ 
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On a par exemple 

Posons de même 

Les Ç., d'une part et les r^^ de l'autre sont deux à deux imagi- 
naires conjuguées; d'une façon plus précise je dis que E. et i^_., 
sont imaginaires conjuguées. En effet, on peut écrire 

^k — / Çp_.,e / .' — / ^p-.e 



Cette seconde valeur de p/^ se déduit de la première en chan- 
geant i en — i et Ç.. en Ep— .: il faut donc que Çp_y se déduise de E., 
par le changement de / en — i; donc Ç^ et $p_.. sont imaginaires 
conjugués. On démontrerait de môme que r,., et T,p_., sont imagi- 
naires conjugués. Si p est pair, Ç^ est réel, car il est à lui-même 



• • 



son imaginaire conjugue. 

(j \ 2 
' ^ j sera un polynôme du second degré en ;^., qui ne ren- 
fermera d'ailleurs que les termes de la forme ^'.^ ip-v. Kn effet, T 
et R sont symétriques par rapport aux pf^ et aux o*^; si Ton per- 
mute circulairement les p^ d'une part et les o-/,. d'autre part, c^est- 
à-dire si l'on remplace ^.. par ^.. e-'':\ ^^ est multiplié par le même 
facteur. 

De même pour yi^ et r/.^ . Par ce changement ni T ni R ne doivent 
être modifiés : il faut donc qu'un terme ^'^ Ç'.^/ se transforme en 
lui-même ; donc il faut que le changement de variable multiplie 
ce facteur par e-''*^' ; il faut donc que 

21 (y -f- v') == 2ih/j ; - 

donc 
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En appliquant toujours le môme raisonnement, on voit qu'en 
définitive ï et II ne doivent contenir que des termes de la forme 

c' c' 7' y' 

\ ^P- Y ^-l 'iP— i 

Ç; %p - Y Ç; '/ip — y 

*9 Y ^^P — Y '^'V ''*P ~ Y 

Si maintenant on regarde de plus près la formation de la fonc- 
tion II, on voit que Ton peut poser 

II^II,+ IÏ,+ ... + IÏ^ + ...+ II,, 

en écrivant IL, à la place de 

\'.. l'p _ ,+9.(0 (;, y/^, _ ,+ Ip _ ,r/,) -f- o>^;, Ip _ ,+7/, 7/^, _ ,+ a^^Sp _ , 
+ [x j^_L.^.^^^_.+ M, (;,7,p^,+7,,;^,_,) + ^^r,,7,p_, j. 

L, M, N sont donnés par les formules 



I 



VI / sin-y^ cos^'!^ cos-y} \ 

' jiij\ 4si"'^ s>- sin '{^ /' 



^'f'i cos ij; 



\j VI sin ayo co 



^, , sin-y^]; cos" tjy sîn-y^ 



s(- 



2 sin *A 4 siu i 

Les équations de Lagvange seront alors 

d /DII,\ DH 



/ OU , \ OH., 



dt 

et ou aura en supprimant les indices 

i" — l (to- + 2) — 2(.)r,'= [x [14 + M/i], 

r/'+20)=' =ui[M$ + Nr,]. 

L'équation en "k est 

().^ _ 0,^ _ 2 — 1-ui) (À-^ + Nu) — ([A-M- — 4o>-).-) == o. 
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Discussion des équations. — Quel que soit y, les racines de 
cette équation en X^ doivent être réelles et négatives ; on voit 
que cette condition sera remplie si [x est suffisamment petit ; 

2,3 

Maxwell a démontré qu'il fallait que Ton eût [jl < — ^ — . Dans 

ces conditions [jl/j tend vers zéro quand/; augmente indéfiniment; 
c'est la une difficulté, car la masse de l'anneau précisément égale 
à ^p^ serait négligeable ; mais cette difficulté disparait d'elle- 
même si on suppose les satellites répartis non sur une circon- 
férence, mais dans un certain volume. 

Si les racines de l'équation en )."^ sont négatives, on pourra 
poser 

n étant une quantité réelle, et en prenant 

B=B;»'>', 

on aura 

Ç^= A,e"''^'>, 

On aura une solution particulière si à Ç„ r,.. on ajoute les solu- 
tions, 

•/io = o. 






et on aura alors une solution particulière pour les p^^ en posant 

à cette solution on peut ajouter la solution imaginaire conjuguée; 
on obtiendra ainsi des solutions de la forme 

pA:=AoCOs(/l^ + cp+2A'vQ), ' 
(7A:= Bq COS (lit + 'f + 2kyh), 

Si on suppose un système d'axes entraînés constamment avec le 
satellite M;^ dans son mouvement non troublé, et tel que l'un d'eux 
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coïncide avec le rayon OM^, ^j. et or^ seront les coordonnées du 
satellite, par rapport à ces axes, et on voit que le satellite décrit 
par rapport à eux une petite ellipse. 

On passe de la trajectoire d'un satellite à celle du suivant en 
changeant k en A: -}- i ; le second satellite est donc en retard sur 
le premier, et ainsi de suite; il y a ce que Ton appelle un déca- 
lage égal à yO. Si Ayô est un multiple de 9.7: le satellite de rang 
A H- I sera dans la même position que le premier satellite; celui 
du rang A* + 2 dans la même position que le second et ainsi de 
suite. 

Le mouvement sera comme s'il se transmettait d'un satellite 11 
l'autre avec un certain retard; on peut aussi dire que le mouve- 
ment se transmet comme une onde, dont la longueur serait . 

î 

Les ondes qui sont les plus dangereuses pour la stabilité 

sont les ondes courtes; les ondes qui compromettent le plus la sta- 
bilité sont donc celles qui correspondent aux grandes valeurs 
de y. Considérons en effet une valeur de n, suffisamment grande; 
pour une grande valeur de y> l'onde serait très courte, deux satel- 
lites voisins pourraient se rapprocher d'une manière sensible, et 
leur action mutuelle ne serait plus négligeable par rapport à 
l'action de Saturne. 

La masse totale [xp doit être inférieure à -^; si on suppose 

p= 100, on voit que la masse de l'anneau devra être inférieure à 

2 3 

^ , c'est un inconvénient de la théorie de Maxwell ; mais 

10 000 

c'est un inconvénient artificiel, car Thypothèse que nous avons 
faite est trop simple et il faut supposer une distribution de satel- 
lites occupant un certain volume de l'espace. 

La perturbation apportée au mouvement d'un satellite par l'ac- 

tion d'un autre est de l'ordre de i..^ , en appelant la distance 

des satellites. 

Si les satellites sont répartis sur une circonférence de 

longueur /, on a 0= — : le mouvement sera stable si la per- 
turbation a une valeur moindre qu'un nombre donné A. Dans 
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rhypolhèse de satellites répartis sur une circonférence, il laut 
donc que l'on ait ' 



a/;' 



ou bien 



<A 



Af 

i ^ ri 
P 

Si les satellites sont répartis sur une couronne comprise entre 
deux cercles de distance comparable h la longueur de l'un deux, 
le nombre des satellites sera comparable à jr. S'ils sont dans un 
anneau de volume fini, leur nombre sera comparable à y/, la 
masse de l'anneau sera égale à ^p^ et la condition de stabilité 
sera que Ton ait 

la masse de l'anneau devra donc être inférieure à un nombre qui 
ne dépend pas du nombre des satellites. 

Stabilité de l'anneau supposé liquide. — Nous allons main- 
tenant, d'après Maxwell, examiner de plus près la stabilité du 
mouvement, dans le cas où l'on suppose l'anneau, soit fluide, 
soit composé de corpuscules indépendants. Dans cette partie du 
travail de Maxwell, les raisonnements ne sont pas tout à fait 
rigoureux, ni même très clairs, mais on doit en accepter la con- 
clusion; il est probable qu'un calcul rigoureux conduirait à des 
résultats peu différents de ceux auxquels nous arriverons. 

Soient .r, y, z les coordonnées d'un point dans le mouvement 
normal; //, r, (v, les composantes de sa vitesse ; 



ù 



l^ dr, . r/r 



"+rfr' ^'+^' "+7r' 

les mêmes éléments dans le mouvement troublé, 
l.a force vive du svstème est 



'"i"' 




r- / «+-^ + •'+ 
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rénergle potentielle — U est égale h 






eu appelant V et p le potentiel et la densité dans le mouvement 
normal, V + V, f + p', ce que deviennent ces quantités dans le 
mouvement troublé. V est une somme de quatre intégrales : 

La première intégrale est une constante ne dépendant pas 
de ^, Y,, ^; la seconde et la troisième sont égales, comme on Ta 
vu au commencement du cours. Elles représentent la moitié de 
Ténergie potentielle due à l'attraction de la matière répartie, 
comme dans le mouvement normal, sur l'excès de matière qu'il y 
a en cluujue point dans le mouvement troublé. 

On peut poser en négligeant le cube de ;, y», Î^, 



v+v»,v+î-^+,4^+"^^ 



(\r dtj i)z 

La somme do la seconde et de la troisième intégrale peut donc 



s cerne 






od'z. 



di/dz I ' 

Il faut maintenant calculer la dernière intégrale; je peux sup- 
poser que le déplacement ^, y,, i^ est la somme de deux déplace- 
ments, l'un dans le sens de la propagation de ce mouvement, 
l'autre transversalement \\ ce sens : on posera ainsi : 

>- _.../• __i /. 

'fi - '^.i4-'^,, 

_ _ V \ Y 






aoi AyyEAl' DK SATl/tyiL 

Ou sait (railleurs, d'après les formules connues des vibra- 
tions longitudinales et transversales que 

est une difFérentielle exacte et que Ton a 



(l.r (fj/ (fz 

On a d'ailleurs la relation 

, a SC; il 07. ' fi Z^ 

' d.v ajf (fz 

Dans ces conditions on peut écrire 






YoWt-::— / ov/t r.+y.-. + rV- 



Prenons maintenant des axes mobiles tournant avec une 
vitesse co autour de Taxe de Tanneau, Texpression de V ne 
change pas. ï devient dans ces conditions 




— <<//+ -0 |'+[''+7J/' +'•)(,.+;} 



+(«■+§)'; "^' 



La condition d'équilibie est que 



rr—vuii 



soit maximum; si Ton tient compte de la relation 
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on trouve les éqiuitions 



0; ' ù.f ' 



r. 



./' + .,o4'=__ + 



•>o *V ' 



_ .1P dO 



P est un polynôme du second degré en (;, '/i, î^}, 8 est une 
variable auxiliaire liée aux trois premières par la relation. 

O.r 0// Oc 

Conclusion. — Maxwell est arrivé au résultat suivant : pour 

(fue l'équilibre soit stable, on doit avoir : ?<—, — . 
* ^ ^ ooo 

Ce résultat est vrai pour un amas de poussières cosmiques 
comme pour un anneau liquide ; mais il y a une condition de plus 
pour un anneau liquide : la pression à la surface doit être dirigée 
vers l'intérieur ; or, Laplace a trouvé pour l'anneau à section 

elliptique que la densité doit être supérieure à --, pour l'anneau 

extérieur et à 2 pour l'anneau intérieur. 

On peut alors se demander s'il y a d'autres figures d'équilibre. 
Rappelons que l'on a trouvé qu'il fallait que 



o 



to" < 27:0. 



\ 



Dans cette formule, co doit être pris égal à la vitesse angulaire 
d'un satellite dont l'orbite coïnciderait avec l'anneau. On a donc, 
en désignant par M la masse de Saturne et par L le rayon de 
l'anneau. 

.. M 
Si p, est la densité de Saturne, R son rayon, on a 
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On trouve donc 

4 H' 






Le second membre est supérieur à — ;- : si on rapproche ce ré- 

ib 

sultat du résultat précédent, on voit que ranneau de Saturne ne 
peut être liquide. On a d'ailleurs vu plus haut que les observa- 
tions montraient qu'il était probablement composé d'une masse 
de petits corpuscules. 
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ÉVREUX, IMPRIMERIE DE CHARLES HERISSEY 



ERRATA 



Page 3o, ligne i5 en remonlant. — Au lieu de — w'^o;, il faut — to'gS. 
Page 37, ligne 10 en remontant. — Au lieu de -^-jz — , il faut . 

Page io3, ligne 3 en remontant. — Au lieu de — r-, il faut — 



2r* ar* 

Page 104, lignes 6 et 10. — 

Au lieu rfe — 6 V^ CxYp'd-z', il faut 3 V -^ fxyp'dx. 
Page io5. — La dernière formule de la page doit être r mr^ ( i + 



e,Y 



3 

Page 107. — Les dernières formules de la page renferment quelques erreurs 
de signe faciles à apercevoir, d'où il suit pour la formule finale 



4 = '-^+ 

00 ^ 



t(-^: 



Page 124» ligne 10 en remontant. — Le numérateur de la fraction sous le 
second radical doit être a^ — c^. 

Page 128, ligne 12, — Lire S"-|-(Hp2-|- K) S = 0. La formule suivante 
doit être R"S— S"R =z o, d'où en intégrant R'S — S'R = a. 

Page i34, lignes 12 et i5. — Lire S'R — R'S. 
Page i37, ligne 7. — Mettre le signe — devant /. 

Page 187, ligne 7. — Dans le dernier terme de la ligne, il faut x ■\- x' au 
lieu de x — x\ 
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